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PREFACE TO RUSSIAN EDITION 


THE theory of electron-atom collisions is given in many monographs, text- 
books and reviews. The “‘Theory of Atomic Collisions” by Mott and Massey 
is the most substantial among these publications (the second edition was 
issued in 1950). Massey’s articles ‘Theory of Atomic Collisions’? and 
“Excitation and Jonization of Atoms by Electron Impact”’ in Handbuch der 
Physik (1956) are also concerned with this subject. Yu. N. Demkov’s 
‘“Variatsionnye Printsipy v Teorn Stolknoveni’ (1958) is devoted to 
variational methods. 

However, the theory of atomic collisions has recently developed very 
rapidly in directions which are either not considered in these books or are 
only briefly mentioned. For example, substantial progress has been made 
in the application of dispersion relations to atomic collisions. Important 
results have been obtained on the extremal properties of the scattering length 
derived by variational methods. Computers have facilitated extensive 
calculations of scattering cross-sections with allowance for such factors as 
the polarization of the target atom by incident electrons, the coupling 
between different channels and so on. Calculations of this kind are being 
carried out at a very intensive rate both in the USSR and abroad. It has 
been found that the factors just mentioned are in many cases very important. 
The semi-classical theory of electron—atom collisions, first developed by Niels 
Bohr in connection with the interaction of fast particles with matter, has 
been developed still further. 

The present book ts largely devoted to recent developments in the theory 
of electron—atom collisions. It expounds both the general theory and its 
applications to specific problems. Chapters | and 2 are devoted to central 
force problems. The general theory of electron—atom collisions is discussed 
in Chapter 3 in the simplest case, 1.e. the case of the hydrogen atom. 
Equations describing the collision of an electron with an arbitrary atom or 
ion are considered in Chapter 8. The remaining chapters are concerned with 
concrete calculations. 

The general theory given in Chapters 3 and 8 is valid at any energy of the 
incident electrons, and covers both excitation and ionization processes. 
However, in concrete applications of the general theory we shall confine our 
attention to the interaction of slow electrons with atoms, 1.e. elastic scattering 
at energies of a few electron volts, and inelastic scattering near the excitation 
threshold. It is precisely this energy range which is of the greatest interest in 
connection with the study of plasma processes and many other phenomena. 
On the other hand, approximate methods based on perturbation theory are 
not valid in this energy range, and a more complete analysis becomes neces- 
sary. For obvious reasons, most attention will be given to the interaction of 
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electrons with the simplest systems, i.e. the hydrogen atom, the helium atom 
and the helium ion. 

The account which we shall give of concrete calculations of collisions 
between electrons and the simple systems mentioned above is not intended 
to be exhaustive. We shall be mainly concerned with recent calculations in 
which the various complicating factors are taken into account. The results 
of some earlier papers will also be summarized. Ionization processes will 
not be considered. 

I wish to thank Academician V. A. Fock, Yu. N. Demkov and L. D. 
Faddeev, whose comments have contributed to the improvement of this 
book. 


Leningrad 1962 G. DRUKAREV 


PREFACE TO ENGLISH EDITION 


As the study of atomic physics progresses, the scope of textbooks and mono- 
graphs appears to become more limited, but correspondingly their depth 
increases. Editorial failure to understand and keep pace with this involution 
process would lead to a paucity of books which describe physical behaviour 
in sufficient detail or depth. 

With depth of treatment in mind Professor Drukarev seeks only to 
discuss the quantum theory of electron—atom collisions, omitting detailed 
reference to other atomic and molecular collision processes. Since the 1949 
edition of Mott and Masscy’s standard work “Theory of Atomic Collisions”’ 
few attempts have been made to cover the more general field, whilst all the time 
the subject has advanced, in principle, in magnitude and in refinement. 

Important contributions to these advances have come from the University 
of Leningrad, where, under the aegis of V. A. Fock, Drukarev and others have 
concentrated their studies on electron collisions, developing the brilliant 
team who have established the neighbouring school at Riga. 

“The Theory of Electron-Atom Collisions” should find a place not only 
with the professional quantum theorists but with those physicists who require 
merely to understand the effects produced by elastic and inelastic collisions, 
without personally manipulating the mathematics. No longer is it possible 
for the more classically inclined ionized gas physicist to aver that the quantum 
theory methods only nibble at the collision problem. Consider the achieve- 
ments as of 1964: the ability to predict elastic scattering and excitation 
functions for simple atoms, with pointers in the direction of success for 
many-electron atoms; the prediction of the existence of virtual and real 
energy levels, and the scattering effects they produce; the establishment of 
upper limits with confidence which is far greater than that usually applied 
to theoretical cross-sections; the prediction of scattering cross-sections for 
electrons with ions and with excited atoms, collisions which have hardly yet 
been attempted in experiment. Only in ionizing collisions is the theoretical 
position rather less satisfactory. 

The successes have come from such treatments as are described by 
Professor Drukarev: perturbation theory, distorted wave calculations, 
variational methods, allowances for exchange and for coupling, together 
with improved computational techniques. 

All this has made it well worthwhile to undertake the translation into 
English of “‘The Theory of Electron-Atom Collisions”, and a debt of 
gratitude is owing to Mr. Chomet for his capable efforts. 


November 1964 J. B. HASTED 
University College London 
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Chapter | 


THE SCATTERING OF A PARTICLE BY A CENTRE OF FORCE 


|. TIME-DEPENDENT THEORY OF SCATTERING 


THERE are two ways of approaching the scattering of a particle by a centre of 
force. Onc of these 1s to be found in ncarly all the existing textbooks of 
quantum mechanics; it involves the consideration of the continuous-spectrum 
wave function with a known asymptotic behaviour in a given field. This 
corresponds to the scattering of a stationary beam of particles with well- 
defined momentum at a large distance from the scattering centre, and 
represents an idealization of the real state of affairs. 

The other approach is used in the time-dependent theory of scattering 
{1-5] and involves the study of the development of the process with time. 
Here one investigates the collision of a single particle with the ccntre of force 
rather than the scattering of a stationary beam of particles. The state of the 
particle is represented by a wave packet with a certain momentum distribu- 
tion. This is closer to the conditions obtaining in real scattering experiments 
which nearly always involve the use of various collimating devices. 

Any scattcring experiment can be divided (somewhat schematically) into 
three parts: the formation of the state of the particle, its interaction with the 
centre of forcc, and finally, the determination of its state after the intcraction. 
When the mcasuring and forming deviccs are at a large distance from the 
region of the interaction between the particle and the centre of force, the 
momentum distribution in the wave packet is independent of the distance 
traversed by the particle. We shall simplify the real conditions by assuming 
that the statc of the particle is formed at time ¢ = —oo while the final state 
determination is carried out att = +00. A definite limiting momentum dis- 
tribution will be supposed to exist in the wave packet as tf + —oo and 
t— +00. The problem 1s then to establish the relation between the initial 
and final distributions. 

The expressions for the scattering amplitude in terms of the wave functions 
are, of course, the same on either of the two approaches mentioned above. 
However, the time-dependent theory is in some respects more satisfactory. 
It is particularly advantageous in the case of scattering of a particle by a 
complicatcd system, for example an clectron by an atom; and also in the 
analysis of the asymptotic behaviour of the wave function for the system 
along differcnt directions tn configuration space. The time-dependent treat- 
ment is also more convenicnt in the analysis of certain general problems. 

We now proceed to the mathematical formulation of the thcory, using the 
system of units in Which #4 = m = 1. In this system the dimensions of all the 
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parameters can be reduced to the dimensions of length by multiplying them 
by the appropriate powers of A and m. For example, the total energy - and 
potential energy V ofa particle have the dimensions of [L]~* when multiplied 
by m/f?. With a suitable choice of the scale of length it is possible to reduce all 
the parameters to a dimensionless form. The length characterizing the range 
of action of the field V(r) is a natural unit of length. In the ensuing analysis 
of the scattering of an electron by an atom it will be convenient to take for it 
the Bohr radius ay; the system of units is then known as the atomie system. 
Let P(r, t) be the wave function for the particle. In the atomic system of 
units it will satisfy the equation 
cD 
i 71 H® (1.1) 
where the energy operator H may be written as the sum of the energy of 
free motion Hy = —4A and the energy of interaction V between the particle 
and the centre of force: 


We shall assume that V ts spherically symmetric and falls off with distance 
more rapidly thanr~' as roo. The case V~r~* requires special considera- 
tion and is discussed tn § 2.7 

Let @ be the eigenfunctions of the operator H, in the form 


QO (E,v,r) = k*(2n)7# inne (1.3) 


The plane wave e'*"" describes the free motion of a particle with energy 
E = 3k? and momentum k = ky. The normalizing factor £4(2n)~? is chosen 
so that 


| oo(E’,v',r) 0 (E,v, r) dr = 6(E— E’) d(v—v’) 
We expand ® in terms of the eigenfunctions 95: 
M(r, t) = | e FW OE, v, 1) @o(E, v, r) dE dQ (1.4) 


We shall use the symbol dQ to indicate integration with respect to anglcs over 
all the directions of the unit vcctor ». This expansion represents the wave 
packet mentioned above. The description of a system with the aid of the 
function ‘’, is usually referrcd to as the interaction representation. On this 
particular representation the operator //, is diagonal. The ftnction 
VY o(E, vy, t) is the state vector on the interaction representation. It is physically 
reasonable to assume that when ¢ > —oo there exists a definite limit of this 
function, which we shall denote by ‘¥_ (the function ®(7, ¢) has no limit when 
t+ —o). When t > +0, the state vector ‘¥,(£, », 2) tends to a definite 
limit which we shall denote by W,. The existence of these limits reflects the 
fact that @(r, 7) satisfies the equation 7dM/ct = //)>@ for a free particle when 


Tt Some of the results given below are valid only when the requirements imposed on V 
are more stringent. This will be pointed out at the appropriate places in our analysis. 


1. THE SCATTERING OF A PARTICLE BY A CENTRE OF FORCE 3 


t— —oo. The relation between the final state ‘, and the initial state ‘V_ 
may be expressed in the form 
W .(E,v) = i S(E, v, E', v') _(E', v’) dE’ dQ’ (1.5) 


The matrix S which is defined by this relation is of fundamental importance 
in the theory of scattering. It follows from 


[\¥.P dE do = [ |¥_/?dEdo (1.6) 
that the S-matrix is unitary, that is, 
{ S*(E, v, E’, v’') S(E, v, E", v") SE dQ = S(E'—E") S(v'—v") (1.7) 


By considering the time-reversed motion, i.e. the situation when ¢ is replaced 
by —¢ and M(t) by ®*(—1), it 1s possible to show that the S-matrix 1s sym- 
metric with respect to time reversal; that 1s, 

S(E,v, E’,v’) = S(E’, —v’, E, —v) (1.8) 


where the right-hand side corresponds to the time-reversed motion. 

In order to obtain more detail about the structure of the S-matrix, let us 
consider the Heisenberg representation in addition to the interaction repre- 
sentation, and establish the relation between them. 

Let ~ be the eigenfunction of the operator H which satisfies the cquation 


Ho = Ep 
or in another form 
Ag+(k?—2V)o =0 


We shall consider those functions g which have the asymptotic form: 


I 
Pro ~90+0(~) 
and are normalized to 0(E— E’) d(y—y’), as is @. We shall see later that this 
does not define the function g unambiguously. 
Let us expand ® in terms of the functions g: 


(r,t) = | e~ W(E, y) o(E, v, r) dE dQ (1.9) 


where the nght-hand side includes summation over the discrete spectrum. 
The function 1s the state vector on the Heisenberg representation. 
The relation between p and ¥ is: 


Y o(E, v, t) ils w(E,v, E’,v') VCE’, v') dE’ dQ’ (1.10) 
Comparison of the expansions (1.4) and (1.9) shows that 
W(E, v, E’,v’) =| ob(E. VN) Q(E'7,v,rdr (1.11) 
It will be convenient to transform (1.10) and (1.11) to another,form. Let 
P=Potz (1.12) 


4 THE THEORY OF ELECTRON—ATOM COLLISIONS 


so that in view of (1.2) we have 
(GI—E)y =—V@p (1.13) 


In order to write the solution of (1.13) in an operator form let us introduce 
the resolvent 
G(E) = (H-—E)7! (1.14) 


The operator G is uniquely defined for complex E, since in this case the 
operator H—E has no zero eigenvalues. The ambiguity arises in the case 
of real E. If we regard G(E) as lim G(E+ie) then the result of the appli- 
e—> +0 

cation of G will depend on the sign of ie. It will be shown later that this 
sign, and therefore the resolvent itself (and thus the function @ also), are 
unambiguously defined by specifying the initial condition for the state 
vector ‘V,(E, v, t) as f + —oo. When expressed in terms of the resolvent G, 
the solution of (1.13) is of the form 


14=—G(E)V Qo (1.15) 
Substituting (1.15) into (1.12) we have 

~ = Po GV Ho (1.16) 
We can now employ the equation 

G=6,—-G6,VG (1.17) 


which can be obtained from simple operator algebra. In this expression Gy 
is defined, analogous to G, by the formula 


Go =(H)—E)7! (1.18) 
Equation (1.16) can now be written 
~ = Po—-GoVe (1.19) 


which is the integral equation for gy. Substituting (1.19) into (1.11) we 
have 


w(E, v, E'w’) = (E-E) S(v—v')— { p(E¥, 1G E' V(r) CE’, v0) de 


(1.20) 
Since 
| pC, Vodr = [ pV(Gopo)* dr 
and also (for complex E£’) 
Copal Es, 1) = ———— ol Ey 857) = = 0 GOED) 
Oo(E, v, r) = ——« Se he oe 
0ofo ) H—E’ Po : Ep? 
we can, using the substitution 
T(E,¥, Ev) =| o8(E,¥, NV) OE WW dr (1.21) 
rewrite equation (1.20) in the form 
FEAL 
inven yes O-r Gea ES (1.22) 


E-E 
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We note that this expression contains all the singularities of w in an explicit 
form. The matrix 7(E, v, E’, v’) has no singularities at E = E’. If V(r) 
falls off more rapidly than r~*? when r > oo, then T has no singularities 
when v = v’. It will be seen later that the matrix 7 plays an important 
role, since the effective scattering cross-section can be expressed in terms 
of this matrix. We now substitute (1.22) into (1.10), giving 

—i(E’-E)t 


W (E, v, t) = YE, v) — -| T(E, v, E’,v’) —— > W(E',v') dE’ dQ’ (1.23) 


In this expression £ is real and £’ complex. When E”’ becomes real, there 
is a Singularity in the integrand at the point E = E”’. In order to establish 
the behaviour near this singularity, let us substitute the initial condition for 


t = —oo. We shall require that when t = —oo the state vectors on the 
Heisenberg and on the interaction representations should be the same: 
Y o(E, v, — 00) = P(E, v) (1.24) 


In order to establish on this basis the contour of integration, let us transform 
(1.23) so that 


Y o(E, v, t) = P(E, v)— 
-| T(E, v, E’, v’) P(E’" v') — T(E, v, Ev) MEV) ce - Ege dol 


E-—E' 
—i(E'- E)t 
- | TE V, soe v') P(E, vy dQ) eseerae (1.25) 
C 
The expression 
0 = —_Irte, v, E’, v’) ¥(E’, v’)— T(E, v, E,v’) ¥(E, v’) | 


E— E’ 


which enters into the first integral on the right-hand side of (1.25) exhibits 
no singularity at E = E’ and is integrable in the interval 0 < E’ < o. 
Therefore, from the Riemann-Lebesgue theorem we have 


lim | Qe Et dE’ =0 


t~ to 56 


oe” HE’ a 
ie =e | 
E’- 
C 
tends to zero as tf + — 00 if the contour of integration in the plane of complex 


E’ is chosen so that the point E = E’ is excluded, as shown in Fig. 1. In 
fact, assuming that for t < 0 


—(E’—E)|t| =z 


The integral 
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the integral J may be transformed to read 
E|t| 


J= | Sa: 
z 


where the singularity at z = 0 is by-passed, as shown in Fig. 2. When 
It] + co we have 


J= | Saz=0 
Zz 


This method of excluding the singularity corresponds to the definition of 
the resolvents G and Gy as the limits 


G,(E) = lim G,(E + ie) 


and 

G(E) = lim G(E + ie) (1.26) 

e> +0 
MZ 
Ime 
E Re E' Rez 
Fic. 1. Contour of integration in the Fic. 2. Contour of integration in the 
E’-plane z-plane 


We can use the above limits, together with (1.14) and (1.18), to rewrite 
equation (1.16) in the form 
ge =— lim isG(E+ie)@, (1.27) 
e+ +0 
We now consider the expression (1.25) in the limit t— oo. The first 
integral on the nght-hand side vanishes as before in view of the Riemann- 
Lebesgue theorem. However, the integral 


eo (E’—E)t 
ia { pope 
C 


does not vanish as tf > oo. Substituting (E’— £)t = z when rt > 0, we have 


te @) 


e 
= dz 
vA 


— Et 
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When ¢ > o0 and the point z = 0 is excluded by deforming the contour 
of integration into a small circle above this point, we have J = 2zi. It 
follows that 


WY, = P_—2ni i T(E, v, E’, v') WE, v’) dQ’ (1.28) 
Comparison of (1.28) with (1.5) yields 
S(E, v, E’, v') = 6(E— E’)[ 6(v—v’)— 2ni T(E, v, E, v’) | (1.29) 


We have thus expressed the singularities of the S-matrix in an explicit form. 

Let us now establish the connection between the effective scattering 
cross-section and the 7-matrix [5]. We suppose that in the initial state the 
particle has a well-defined momentum near ky which is parallel to the z-axis 
and is localized in the volume U, so that fw = U~*4C exp (ikgz), where 
C = | and 


a, ay. a a a, a, 
—-—<x<—}], —~<y<—)], —~<z<- 
( 2 *) ( 2 4 ( 2 5) 


C = 0 when at least one of the coordinates lies outside this interval. 
In momentum space the wave function will be of the form 


(6) = (2)' yea sinjk,a, sin}k,a, sin 4(k,—ko)a, (1.30) 
TU Ky ky k, 
If we normalize along the energy scale, then 
dk\? 
Y_(E,v) = k(—]) W_(k 1.31 
Ey) = (FE) ¥-0) (1.31) 
Let 
1 sindaé 
d,(S) = z 
(1.32) 
2 \* sindaé 
NO) = | ae 
na 4 
so that 


lim 5,(¢) = 6(¢) 
sie (1.33) 


lim Az(¢) = 6(€) 
Since E = 4k? and U = a,a,a, we have 


, 4 
W_(E,v) = (. . 2754, (kx) 5a,(k,) da,(ke— ko) (1.34) 


erry 
Let j be the particle current density in the initial wave packet. The cross- 
section of an arbitrary packet is a,a,, therefore the total flux of particles 
per unit time 1s 


n= aya, J (1.35) 


2—2 
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Next, the probability of finding a particle moving within the solid angle 
dQ. in the final state in the energy interval dE is equal to 


|\Y .(E, v)|? dE dQ 


Therefore, the number of particles which have been scattered per unit time 
into the solid angle dQ within the energy interval dE is 


dn, = a,a, j|V,(E, v)|? dE dQ (1.36) 
In view of (1.15) and (1.34), 
dng = jl2n [ S(E,v, B’,v') 5,,(ka) 5,(k) Aa,(ki— ko) 
x(k’ dE’ dQ’|? dEdQ (1.37) 
Substituting (1.29) into this expression, we have, for v ¥ v’, 
dn, = (2n)*j! { T(E, v, E,v") 5g,(ki) 5g, (k,) Aa, (ki — Ko) 
x (E—E')/k' dE’ dQ’)? dEdQ (1.38) 
Consider the integral 
J =| T(E, v, E,v) 5,,() 5o,(k)) be, (ki— ko) (E—E')\/k' dE" dQ’ (1.39) 


It will be convenient to replace integration with respect to E’ and Q’ by 
integration with respect to ky, ki, and ki, and replace E’ by 4(ki° +k? + k??). 
The result of these substitutions 1s 


jx | T(E, ¥, E’, v') 5g, (Kt) 5g, (ki) Ae, (ki — ko) 
1 
x SE-AK2 +k +k) Ta dkdk dk; (1.40) 
7 


If it is supposed that a, and a, tend to infinity, we have (provided |k/| = k’) 


fo @) 


k’2\ dk’ 
ke (1.41) 


J=]| T(E,v,E’,v)A,(k —k E-— 
[re )A,,( 0 ( ain 


O 
and hence 


J = k-3T(E,v, E,v')A,.(k — ko) (1.42) 


On substituting (1.42) into (1.38) and taking the limit a, > oo in A,., we 
have, in view of (1.33), 


dn, = j(2x)*|T(E, v, E,v’)|? 0(k—kp)k~3 dE dO (1.43) 


When (1.43) is integrated over a small interval of energy in the final state, 
and the result denoted by dN, then 


dQ 
dN 4 = j(2n)*|T|? L2 (1.44) 
0 
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Hence, by definition of the effective cross-section do: 
dN 2 
do = ar ihe ey, 7 “ee v, E, v’)|? dQ (1.45) 


The total scattering eer ie can be obtained from this expression 
by integrating over all angles; the ae 1S 


(1.46) 


On substituting (1.29) into see and assuming that v” = v’ then: 
i[ T(E, v, E, v)—T*(E, v, E,v)] =2n | IT(E,v, E,v')|2dQ’ (1.47) 
OT 
—Im T() = x | |TP dQ (1.48) 


where 7(0) = 7(E, v, E,v). The result (1.48) is known as the optical 
theorem. Finally, in view of (1.46) and (1.48): 
4,3 


——ImT(0)=o (1.49) 
ee 


2. THE WAVE FUNCTION FOR A PARTICLE IN A CENTRAL 
FORCE FIELD 


Let us now consider in detail the eigenfunction of the energy operator H 
defined by equation (1.19) and the relation (1.26). It will be convenient to 
replace @ by another function w defined as 


WE, v,r) = (27)?k~* @(E, v, r) (2.1) 
In view of (1.3), (1.19) and (1.26) the integral equation for ¥ Is 
Ww = el" — lim Go(E + ie) Vw (2.2) 
e> +0 


On the coordinate representation the resolvent is an integral operator 
centred upon a Green’s function. The Green’s function corresponding to 
the resolvent Go(E+ie) will be denoted by G(r, r’, E+ie). In some cases 
we Shall omit E-+-/e and write simply G)(r, r’). This will be done whenever 
there is no danger of ambiguity. In this notation (2.2) may be rewritten in 
the form 
wir) = eh | Golr Vr) Wer’) dr’ (2.3) 
Let us determine the Green’s function G)(r,r’,) corresponding to the 
resolvent Go(E+ie). In the momentum representation: 
1 eik'rew’ e iki’ pte dk’ dQ' 
G,(r,r’) = lim ail SS ee 
" +0 (27)" $k’? — (4k? + ie) 


If we complete the integration with respect to Q’ and substitute 
Ir-r'|=p;  k’'p=z'; — (k? +2ie)p? = z7(e) (2.5) 


(2.4) 
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we obtain 
1 e!* dz r e' dz 1. « 
G, = lim —— ———— ——§— | = lim —e” 2.6 
olin al [aot [eas] 7 lta” 29 
and finally 
1 
Go = ——_ eit lr 2.7 
Y 2n|r—r'|° ee 


Therefore the integral equation (2. : becomes 
elklr- r’| 
gk Vw dr’ (2.8) 
|r-r| 
For large r the asymptotic behaviour of w is 
ikr 


wn~ ek + aly, n)—, where n = : (2.9) 
and the amplitude g is given by 
q(v,n) = — = clad Viv, vr’) de’ (2.10) 
Comparison of (2.10) with (1.21) yields 
jis — OO rE, n, E, v) (2.11) 


With the aid of (2.11) and (1.45) we can express the effective cross-section 
in terms of g: 


do = Iq|? dQ (2.12) 
The optical theorem (1.48) may now be written in the form 
k 
Im q(n, n) = al (n, vy)? dQ (2.13) 


The effective scattering cross-section given by (2.12) is obtained in the 
Stationary theory by considering the asymptotic behaviour of the wave 
function (2.9). This is done by finding the current density associated with the 
scattered particles, i.e. with the term g e’/r in (2.9). This current density 
turns out to be equal to klg|*.. The quantity do/dQ which is equivalent to 
(2.12) is obtained by dividing this current density by the current density 
of the incident particles, which are described by the plane wave exp (iky.r). 
The latter current density 1s, of course, equal to &. The time-dependent 
theory provides the justification for this procedure. The expression (2.13) 
can also be obtained directly by considering the asymptotic form of the 
wave function (2.9). This can be done by setting the current through a 
sphere of arbitrary radius equal to zero, i.e. 


[ (vy*—v*vy)ds =0 
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This condition represents the conservation of the number of particles on 
scattering. 
Consider the wave function wy at r = 0. On setting r = 0 in (2.8) we have 


i 
pie =f elt ss Vb dr’ (2.14) 
TC r 


The second term in this expression is a measure of the distortion of the 
wave function at r = 0 due to the field V(r). For large & the value of this 
integral is small owing to the rapidly oscillating factors exp (ikr’) and w. 
The maximum distortion occurs at kK = 0. In the latter case (2.14) 1s. of 
the form 


WO) = 1-2 | VO) WO’) dr’ (2.15) 
0 
Applying the mean value theorem to [/(0)—1] we obtain 
0)-1 i 
[¥O)~1] = 2 | |V|r dr’ (2.16) 
Wm 


where ,, is the mean of ||. Therefore the relative distortion of the wave 
@ 
function is characterized by the integral I Vr dr. 


O 
Let us now separate the angular variables in (2.8). To do this, we shall 
expand exp (ikv.r) and w into partial waves: 


re [4n(21+ 1)]#i! - un Yio(9) (2.17) 
where @ is the angle hauatns v and r and 
+ 
tr 
uy) = (=) Ji44(kr) (2.18) 
W =D Lanai nyt? y.(0 (2.19) 


By analogy with spectroscopic notation the functions /O(k,r), f,(k,r), 
filk,r)... are referred to respectively as the s, p,d... waves. It will also 
be convenient to expand the Green’s function Go defined by (2.7) in the 
form 


ae) a ogo 
Golr,r’) =} 2——— ¥,,(0, @) ¥1,(0'9') (2.20) 
Lu 
Y,, are the normalized spherical harmonics, I, is the one-dimensional 
Green’s function 
Pee —uj(kr)o(kr’) . r’ 
—ofkr)u(kr’) re<r 


+ This estimate is adequate only if yw does not change too rapidly over the range of r 
in which V is appreciable. 


(2.21) 
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where 1, 1s given by (2.18) and v, 1s defined by 
1 /nkr\* 
fees (=) H{D,(ker) (2.22) 
Substituting these expansions into (2.8) we obtain the integral equation 


for f(r): 


r \° 8) 


f(r) = u,(r)+ 20,(r) { u(r Vir) fir’) dr’ +2u; | V(r) Vir) fir )dr’ (2.23) 
0 r 
The asymptotic form of f; when r > 00 Is 
| i gost 
fir) ~ *sin Ca — 4 + ge (* 2) (2.24) 


where q, is the partial scattering amplitude given by 
a 
gq = —2[ u()VOfi(r)ar (2.25) 
0 


The amplitude g can be expressed as a sum of the partial amplitudes 
if we substitute the above expansions into (2.10): 


q = ), [4n(2!+ D}*ar¥,, (8) (2.26) 

The total cross-section o can also be expressed as a sum of the partial 
cross-sections ¢;: 

Oo = > 0; (2.27) 

o, = 4n(21+ 1)|q,]’ (2.28) 


The optical theorem can be used to establish the relation between the 
imaginary part and the modulus of q,. Substituting (2.26) into (2.13) we 
have 


, 41 
Since the amplitude g, may be written in the form 
qi = |qi| (2.30) 
there follows 
a l ; 
gq, = eft sin 6, = ae 11) (2.31) 
or 
= (2.32 
11 k cot 6, —ik ae 


We can now use (2.31) to transform the asymptotic expression given by 
(2.24) so that it includes the phase 6,. After some rearrangement we find: 


eft lnm 
fir ae kr — ae 5 (2.33) 
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It is clear that the relationships given by (2.30) and (2.33) define the phase 
only to within a multiple of z. This limitation can be removed by means of 
the natural condition 

lim 6, = 0 (2.34) 


Ea 


The partial cross-section a, is given in terms of 6, by 


o, = 4n re sin’ 6, (2.35) 
It is evident from this formula that 
6, < (21+ Wo (2.36) 
The function /, is the solution of the differential equation 
wf + (2-52 av) j,=0 (2.37) 


where 
FO) =0 (2.38) 


The function f(r) has the required asymptotic behaviour at large r. 
In addition, let us consider the functions which have the given behaviour 
for small r. We shall denote them by x,(r). In some respects they have ad- 
vantages over f(r). Since both functions are solutions of the same differential 
equation (2.37) and both vanish when r = 0, their ratio must be independ- 
ent of rand is a function of k only. We may therefore make the substitution: 


x(k, r) = L,(k) fi(k, r) (2.39) 


and require that the functions x, should behave as u,(kr) for small r. 
Substituting (2.39) into (2.33) we have, after some rearrangement 


r 


xr) = u(r) + [o(r)u,(r’) —u(r)o(r’)|2V (r')x,(r’) dr’ (2.40) 


0 


Lk) = 1-2 ] vp, Vx,dr (2.41) 
0 


We note that the function v, has a pole of the form r~' at r = 0. Moreover 
V has a Coulomb singularity at r = 0. It follows that v,V has a pole of the 
form r~'~* or r~'. However, since f; falls off as r'** when r — 0, it follows 
that the integrands in (2.40) and (2.41) remain finite at the lower limit. 
The function L, defined by (2.41) is identical with the so-called Jost 
function [6]. In view of the connection between the Bessel and Hankel 


functions, we have 
vp = k[(-1)!*fu_y- 4, —iny] (2.42) 
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We can use this relation to rewrite (2.40) and (2.41) in the form 


x, =u,t+(- 1k | Lar a(r’)= (De 20 xr) dr’ (2.43) 


L=14+(-1)k | eg VG dr-+ik | 4,2Vx, dr (2.44) 
0 0 


It follows from the properties of Bessel functions that at given r the expres- 
sions u,k~' and u_,;_,k'*! are analytic functions of k*. Therefore k~ ‘x is 
an analytic function of k? at fixed r. ‘Hence, L, may be written in the form 


L, = Ak) +ik?!* 'Bi(k) (2.45) 


where A, and 8B, are even functions of a real k, defined by 


A, — 1+(- I'k | wy, 2Vxydr 
. ° (2.46) 
B, = ko! [ 1 2Vxrdr 


0 


The partial scattering amplitude g, can be expressed in terms of L,. Using 
(2.39) and (2.46) we can rewrite (2.25) in the form 


k??B(k) 
— = > 2.47 
q1 L, ( ) 
Let 
Ly = |L,| el?! (2.48) 


The quantity g, is found to be closely related to the phase 6, introduced 
above. This can be shown by substituting (2.48) into (2.47) and utilizing 
(2.45): 


| 
q=- pe sin QD, (2.49) 


Comparison of (2.49) with (2.31) then shows that gy, =—6,+22N. The 
term 2x/N is unimportant and can be removed by a suitable choice of origin 
of ~, We shall therefore set 


The properties of L, when k 1s complex are of interest. They can be 
conveniently investigated with the aid of (2.41). By using the asymptotic 
expressions for wu, and v, and the expression for x, given by (2.40) we can 
verify that L, is analytic at least over the upper half-plane. Its properties 
in the lower half-plane depend on the behaviour of V(r) asr—7 oo. If V 
falls off as a power of r, that is V ~ r~", then in the lower & half-plane L 
is not analytic. If, on the other hand, V ~ e “/r, as is often the case, 
then L, is analytic in the lower half-plane within the band |Im k]| < 4a. 
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If the potential energy V is negative, then the existence of bound states 


with energies —|£,|, —|E 2]... is, in general, possible. Negative energies 
correspond to imaginary wave numbers k; = ik; = iV 2)E)). It may be shown 
that the function L,(k) vanishes at the points ix, ik2,.... Let us substitute 


Ak = ik into (2.40) and require that x,(ik, r) vanishes when r > oo, aS in a 
bound state. In view of the asymptotic behaviour of the Bessel functions 
of complex arguments, we find that as r > oo the coefficient of u,(r) must 
vanish, which in turn implies that L, = 0. It may be shown that all these 
roots are simple ones. 

When V(r) behaves at infinity as a Coulomb potential, as for example 
in the case of the scattering of an electron by an ion, the above equations 
must be modified in view of the particular asymptotic properties of the wave 
functions for the continuous spectrum of the Coulomb field. To make 
this modification, let 


V(r) = Wr) - ~ (2.51) 


considering the attractive interaction in which we shall be interested later. 
We introduce the Green’s function Gc for the Coulomb field: 


Z 
(44 + 4k* 4+ *) Go = —-O(r—-r’) (2.52) 
Instead of (2.3), we now have 
Y =Ve—| GoW par’ (2.53) 


Where the function We is the wave function for the continuous spectrum 
of the Coulomb field whose asymptotic behaviour for r > © Is 


We as eilkv.rty In (kv.r kr) Ge 1 eikrt+y In 2kr (2.54) 
r 
where 
Z 


and gc is the scattering amplitude in the Coulomb field and is given by 
y T(it+iy) 1 fy In 2 sin? © 
[os er eT € : 

2k TU—iy) . 0 

sin 7 


(2.56) 


The problem is to find the solution of (2.53) which has the asymptotic form 
given by (2.54) but with an amplitude corresponding to the potential energy 
(2.51). In order to separate the angular variables, let 


ur(r) 


r 


Wo =) [4n(21+ 1} Yio(9) (2.57) 
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where uy is the regular solution of 
d*ue Ii+1) 22 
~ + om fale ur = 0 (2.58) 
dr r r 


which, when r > oo, has the asymptotic form 


1 l pal he fre 
ur ~ sin (i — = + yin 2kr + eG ne ) (2.59) 


where qc,, is the partial scattering amplitude in the Coulomb field: 


ri+1+iy) 


= 2.60 
4et = T+ 1— iy) oo 
We shall take the Green’s function Ge in the form 
Pir, r’) es eee 
Gelr.r') = Y 2 ——— FiO, 9) Yin(0'9) (2.61) 
lp 
where 
ul (r)ur(r’) ae 
D(r,r’) = 2.62 
(rr) ee eee ( ) 


where wr is defined by (2.58) and (2.59), and cf is a non-regular solution of 
(2.58) whose asymptotic form is 


I 
i( kr +3 In 2hr 


ae 2 (2.63) 
We write the function w in the form of a series analogous to (2.57): 
all) » 


For the functions /,; we obtain 
io 0) 


flr) = uf(r)+2 v r(ryul(r Wor fit’) dr’ + 2 fu r(ryelr Wr’) f(r’) dr’(2.65) 


r 


For large r we have 


ri eet 2kr 


| l oh 
fir) ~ i sin (11 — = + yln 2kr) + (dcitde ( (2.66) 


where 


oO. 


je | uN WO f()) dr (2.67) 
0) 


It follows that the partial scattering amplitude q, in the field (2.51) ts 
given by 
N= Icei1t 4d (2.68) 
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If we now let 


| . 1 
= ar a I), = dea = ral 1) (2.69) 
and 
Op = Nt (2.70) 
we obtain 
j, = ES 7! (92! _ 1) (2.71) 
2ik 


When r — © we have, using (2.66) and (2.69), 
Lex l 
fir~ esin (ir — = + yin2kr + 5) (2.72) 


The integral equation (2.65) will now be transformed to a form analogous 
to equations (2.40) and (2.41). To accomplish this, equation (2.65) will be 
rewritten in the form 


fir) = uf(0) [ +2 f of KW) ir| i 


+2 LF UEO') =F UEOIWO IAC) dr’ (2.73) 


If we substitute 


yk.) = L(A f(k, r) 
we obtain 


oD 


vir) = uF) +2 | LoFOuF()= oF OME (dr (2.74) 


0 
IS(k) = 1-2 f voWy, dr’ (2.75) 
Consequently, , 
—2 fugwy, dr 
oleae eee (2.76) 
1-2 | of Wy, dr 
0 


3. SCATTERING AMPLITUDE AND PHASE 


It is evident from equations (2.47) and (2.50) that the properties of the 
amplitude g, and phase 6, are determined by the quantity L,. Using equations 
(2.48) and (2.45) we can express cot 6, in terms of 4, and B;: 


k7**) o6t6, = — tiie (3.1) 


18 THE THEORY OF ELECTRON-ATOM COLLISIONS 


Hence it follows that cot 6, is an odd function of k. Since 6(+00) = 0 we 
find that the phase 6, 1s itself an odd function of k, 1.e. 


6(k) = —6(—k) (3.2) 


On this definition 6, will in general experience a discontinuity at k = 0. 
Since L, vanishes on the imaginary axis at points corresponding to bound 
states, we can establish with the aid of (3.2) an important relationship 
between the number AN, of levels with given / and the quantity 6,(+0), 
i.e. the limit of 6,(k) for k — 0 on the side of positive k. This relationship is 


6 +0) = Nn (3.3) 
and is known in the literature as Levinson’s theorem. In order to prove 
this result we recall that, according to Cauchy’s theorem, the change in 
phase of L,(k) for integration over the contour shown in Fig. 3, 1.e. the 
quantity —[6(—0)—6(+0)]/2z, is equal to the number of zeros of L,(A) 
within the domain enclosed by this contour in the plane of complex k. 
All the zeros of L, lie along the imaginary axis and their number if Nj. 
Consequently, 

—[6(—0)—6(+0)] = 22N, (3.4) 
If we now combine equations (3.2) and (3.4) we obtain the result given 
by (3.3). 
Imk 


Rek 


Fic. 3. Levinson’s theorem. The points on the imaginary axis are the zeros of the 
function L 


A special case arises when L,(0) = 0. This situation occurs, for example, 
when one of the discrete energy levels merges with the boundary of the 
continuous spectrum, i.e. tends to zero as the depth of the potential well is 
reduced. Under these conditions L,(0) = 0 and the phase is given by 


§(+0) = Nx +> (3.5) 


where NV is the number of levels (not including the special level with zero 


energy). 
Let us consider the limiting value of g, as k ~ 0. We shall start with the 


expression given by equation (2.47). Substituting equation (2.45) into 


(2.47) we obtain 
Ak) \~? 
n= (ik ao | (3.6) 
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Suppose that A,(0) 4 0 (when 4,(0) = 0, it follows from equation (2.49) 
that L,(0) = 0, i.e. we have the special case mentioned above). It follows that 


lim gq, =0 when / #0 
k-0 
B,(0) (3.7) 


i.e. all the partial amplitudes except the zero-order amplitude vanish and 
the amplitude of the s-wave is, in general, the only one to have a non-zero 
value. The limiting value of this amplitude is usually denoted by —a, where 
a is the scattering length. Thus, 


B,(0 
A (0) 
According to equation (3.1) we have 
] 
lim k cot dg) =— - (3.9) 
k-0 a 


We can use the latter relation to show that the scattering length is equal, 
apart from the sign, to the derivative of the phase angle 6 with respect to 


the wave number kK, 1.e. 
a= (=) (3.10) 
~—— \dk/x=o 


If there is an s-level with zero energy, then the amplitude q (A) increases 
without limit as 1/k when k — 0 (resonance scattering). Under these 
conditions, equation (3.10) is no longer valid. If there is a p-level with 
zero energy, then 

A,(k) 

k*B,(k) 


lim 
k-0 
is finite and therefore 
lim g, #0 
k-+0 
The partial amplitude for /= 1 has a non-zero value when k = 0. 
Partial amplitudes for / > 1 vanish independently of whether a zero level 
with given / exists or not. 
Let us consider in greater detail the sign and magnitude of the scattering 
length. From the definition given by equation (3.8) and the expressions in 
(2.46) we have 


2 | Vx,rdr 
q = (3.11) 
142] Vx dr 
0 
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In a repulsive field V > Q. For zero energy, the function x9 satisfies the 
equation 


2 a r+2| (r—s)Vxo ds 
% 


and is therefore everywhere positive. It follows that the scattering length a 
is positive. In an attractive field V < 0, and if there are bound states the 
function x9 has no nodes (apart from the origin) and is everywhere positive. 
Therefore the integral in the numerator of (3.11) is negative. The denomina- 
tor remains positive as V passes from positive values through zero to negative 
values. The denominator in equation (3.11), which is equal to A,(0), 
vanishes (and the scattering length becomes infinite) only when a bound 
state with zero energy appears. In an attractive field the scattering length 
is negative in the absence of a bound state. In the presence of bound states, 
the scattering length can be either positive or negative, since the sign and 
magnitude of the integrals in equation (3.11) depend on the number and 
position of the nodes of the function xg. All the possible cases can be 
conveniently illustrated by introducing a parameter a, which characterizes 
the intensity of the field of force, and considering the dependence of a on «a. 
Let 
V(r) = aV(r) (3.12) 
and assume that V, > 0. On this definition a positive « corresponds to 
repulsion and a negative « to attraction. 
On the basis of the general properties of functions satisfying Volterra’s 
integral equations, we can state that for given A and r the function Yo will 
be an analytic function of « It follows that a will be a meromorphic 


function of «. 
Figure 4 shows the general character of the dependence of a on «. The 


‘e) 
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Fic. 4. Dependence of the scattering Iength and the energy levels £o, £,,... (broken 


curves) on the parameter a 
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points &, a,,... are the values at which the first, second, ..., levels with 
zero energy appear as the depth of the potential well increases. The broken 
curves represent the dependence of the energy of the first, second, and so on, 
levels on the depth of the potential well. As can be seen from Fig. 4, the 
scattering length, and with it the scattering cross-section at zero energy, 
become equal to zero (Ramsauer effect) for certain values of the parameter a. 
As an example, consider the dependence of the scattering length of an 
electron in a Thomas-Fermi field due to a nuclear charge Z [7] (see Fig. 5). 
Figure 5 is a mirror image of Fig. 4, since —Z plays the role of a. 


Fic. 5. Dependence of the scattering length for electrons (in units of 0-885/Z) on the 
nuclear charge Z in the Thomas-Fermi field 


Let us now consider the behaviour of the amplitude g, as a function of 

k for small k. If the potential energy falls off as r — oo more rapidly than 

any negative power of r, then L, will be an analytic function not only in the 

upper half-plane of complex & but also in a band in the lower half-plane. 

Consequently, A, and B, can be expanded into a series in powers of k 

(within the limits of convergence), and in view of the fact that A, and B, 

are even, these will in fact be power series in k* when k is real. We can 
then write 

A(k) Af+k7AM 4... 

B(k) BO +k?B{P +... 


If k is so small that B, vanishes in the interval between 0 and k, we can 
expand 1/(Bf+k? Bf") + ...) in powers of k’, giving 


(3.13) 


A,(k) A( ; A) Al BO) 4 
BK) ~ BO ** [po “HOE “s 
Let 
AT. AP APB? 
BO a,’ BO —_ ~ Bz = —4r, (3.15) 


E.A.C. 3 
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so that equation (3.14) may be rewritten in the form 


A, 1 
— = ——4rke+... (3.16) 
B, a, 
On substituting the latter series into equation (3.1) we have 
] 
k?**! cot 6, =— — + 4k7r, 4+... (3.17) 
ay 
We can now use equation (2.32) to obtain the following expression for q;: 
<i 
q, = (3.18) 


l 
—- + ther, tik?! 
ay 


In the special case / = O, the expansion given by (3.17) becomes identical 
with the well-known formula [8] 


1 
kcotds =— = +4k?r, (3.19) 
0 


originally obtained by a different method. The quantity r, in this expression 
is called the effective range of the forces. 

If for given / there exists an energy level such that x, < 1, it is possible 
to obtain a relation between a,, r; and «, by using equation (3.16) and 
assuming that A = ik. Since L,(ix) = 0, it follows from equation (2.45) 
that 


A,(ik) _ 


Bix) (—1)t3R2t! (3.20) 
Therefore we have 
(—1)/k7/*1 = -+ LK?r, (3.21) 
and in particular when / = 0 
- = K(1+4xAro) (3.22) 


0 


If the potential energy falls off as r~"” when r > oo, then L, is not analytic 
in the lower half-plane of complex k, and therefore A, and B, cannot be 
expanded in powers of A. It follows that the expressions given by equations 
(3.19) and (3.22) will no longer hold. 

As an example, consider the energy of interaction between a point 
charge and an induced dipole which at large distances is of the form 

V~— p (3.23) 

r 
This expression 1s encountered, for example, in the case of scattering of an 
electron by an atom when polarization of the atom by the incident electron 
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is taken into account. In this case, a is replaced by [9] 
k 
kcotdy =—- a: | +S Pn (AU a (3.24) 


When / ¥ 0 we have instead of equation (3.17) 


eee) Oe) a 
np 


Therefore as k > 0 the amplitude g, falls off as A for all /. 

Let us also consider the behaviour of the amplitude as k — 0 when the 
field is of the Coulomb type. According to equation (2.68) the analysis 
can be reduced to the investigation of g,, since the scattering amplitude 
gic 1n the Coulomb field is known. We shall confine our attention to the 
case / = 0. 

The dependence of q on k can be analysed with the aid of equation 
(2.76). However, it will be useful to express the functions vy and uy (with 
given asymptotic behaviour) in terms of a function having a given behaviour 
at the origin. Let x, and .x; be two linearly independent solutions of the 
equation 


d*x, 2Z 
Ete («2 + =) 610 (3.26) 
dr r 
where 
X,or 
ae (3.27) 
when r > 0. We then have 
u& = eCox,(r) (3.28) 
ft A 
ypose (= xX, +ikCox, ) (3.29) 
Co 
where 
2nZ/k \* 
oe (Saar) (3.30) 


When k < 27Z, the quantity Cy is approximately equal to (2Z/k)*. It 
follows from equations (3.28) and (3.27) that as k — oo the function u& 
increases as k~*. According to equation (2.74), the function u,(r) should 
also increase as kK-* when r— 0. It is evident from equation (3.25) that 
v© falls off as k* as k > 0, so that the product v°y, remains finite as r > 0. 
The quantity q, increases as 1/k when k — 0 (this 1s evident from equation 
(2.76)). 

So far, we have discussed the general properties of the phase 6 and 
amplitude g without considering their numerical values. Let us now proceed 
to the calculation of the phase 6. In principle, this requires the solution of 
equation (2.37). However, this equation has an analytical solution only 

32 
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for a certain special form of V. In all other cases it can be solved only 
numerically. Equation (2.37) is inconvenient for the evaluation of the 
phase, since the function /; is oscillatory and does not tend to any definite 
limit as r > 00. 

However, it is possible to reduce the problem to the integration of a 
first-order non-linear equation for the function ny,(r), whose limit as r > 00 
is equal to the phase 6, [10], and which exhibits monotonic variation for 
monotonic V. 

Consider the logarithmic derivative 


1 df, 
=-— 3.31 
Os As (3.31) 

It satisfies 

d idi+1 

nag go? =2V + A i (3.32) 

dr r 
and may be written in the form 

dt 


du_)_ 
cosni(r)+(—- 1)! ia ie siny(r) 
dr dr 


ee ee aa 
u,cosn(r)+(— L)'u_y-, sin 7,(r) (3.33) 


if we substitute equation (3.33) into (3.32) and bear in mind the equations 
satisfied by u,and u_,_,, and the fact that 


du_,_ du 
) aS Sey - =(-1)'*!'k (3.34) 
we obtain 
dn 2V , 
7, = —, Lin cos n(r)+(—1)'u_y;-, sinn(r)]’ (3.35) 


For large r the function n,(r) in equation (3.33) may be replaced by »,(00), 
so that by using the asymptotic form of the Bessel functions we have 


d l 
o~ anet sin (ir — 5) cos n,(0c0)+ 


ar 


! d ! 
+ cos (1 — >) sin n(eo)| ar In sin kr — - + neo) (3.36) 


On substituting equations (3.36) and (3.31) we find that the asymptotic 
expression for f/f; which corresponds to equation (3.36) is 


Ir 
fir) ~ const sin kr mes n(o)] (3.37) 
from which it follows that 


n,(00) = 0; (3.38) 
Consider the behaviour of y,(r) for r > 0. 
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Since f, ~ r'*! for small r, we have 
lim gr = 1+1 (3.39) 


r-Q 


On the other hand, it follows from equation (3.33) that for small r 
1—1(1+1)y,(r)(kr)77'~! 


=(1+1 4 
Or OT Endy kn re 
Hence in view of equation (3.39) we obtain 
lim y,(r)(kr)~ 7/7? =0 (3.41) 
r-0 


i.e. numerical integration of equation (3.35) should be carried out subject 
to the condition n,(0) = 0 for all /. 


4, DISPERSION RELATIONS 


Consider the matrix element T(E, v, E, v’). According to equations (1.21) 
and (1.16) 


T(E, y, E,v’) = | o*Voodr — lim otVG(E+ieVoodr (4.1) 
e~ +0 


If we recall the explicit form of the function @9 given by equation (1.3) and 
replace the resolvent G by the Green’s function G(r, r’, E+i0), we have 
in the case of forward scattering (v = v’) 


k 1 , 
T(E, v, E,v) = dl V dr — { evr) Vr)G(r, r')V(r') dr dr’ 


a2 
- (22)? 4n 
(4.2) 


If we multiply the latter expression by —(2z)’/k and let — { Vdr = 27qp, 
we have in view of equation (2.11) 


1 
Q-o= mal ek) Vr) G(r, vr’ )V(r') dr dr’ (4.3) 
TT 


Let 
G-dg=F (4.4) 


and consider the analytical properties of F in the plane of complex energy 
z = E+ie cut along the positive real axis; this corresponds to the upper 
half-plane of complex k. The properties of F are determined by the prop- 
erties of the Green’s function G(r, r’, z). For a sufficiently rapid decrease 
in V(r) at large distances (faster than r~°), the function G has the following 
properties in the z-plane [11]: 

1. G(r, r’ z) is analytic in the z-plane with a cut along the positive real 
axis and has, in general, simple poles on the negative real axis E = —|E;j. 
These poles correspond to discrete energy eigenvalues in the field V. Using 
the eigenfunction expansion of G 

*/ of 
G(r,r’,z)=) YiOir) (4.5) 
j zs—F. 


Jj 
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or in the presence of degeneracy, 


Yair’) 
Grr,z=) "eo 4. 
(r,r’, Z) 2 pas (4.6) 
we see that the residues at the poles z = E; are equal to 
2 YW Rr’) 
2. The limit lim G(r, r’, E+ie) exists for all E > 0. 
e—->+0 

3, G(r,r’, E+i0) = G*(r,r’', E—i0) (4.7) 


4. For large |z| 
IG-G,| < Cie ie 


where Gog is the Green’s function for free motion (2.7); Cg — 0 when E > oo. 


Fic. 6. Derivation of the dispersion formula. £,, Eo, ... represent the bound-state 
energy levels 


It follows from these properties that the integral on the right-hand side 
of equation (4.2) converges absolutely for all z and is an analytic function 
of z everywhere on the cut plane except for the poles of the Green's function 
G; moreover, 


F(E+i0) = F*(E— i0) (4.8) 
lim F(z) =0 (4.9) 
|z] +o 


We now derive the dispersion relation [I1, 12]. Consider the contour C 
in the z-plane (Fig. 6). From Cauchy's theorem 
F(z’) dz’ d, 


l 
Bd oa | “yor SB] 42 ee) 
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where 


1 ; 
u=gad {er Viadr | on™" “Vw idr’ (4.11) 
are the residues of F at the poles z’ = —|E,| (we recall that d; are real). As 
the radius of the circle tends to infinity, the contribution due to the integrand 


on the circle tends to zero in view of equation (4.9). There remains only the 
integral along the boundary of the cut: 


1 f¢ F(E'+i0) 1 f F(E’—i0 
F(z) =— | ———— dE' —- —— | ———— 4.12 
2) aa E’~z : a E'-z BS gis een) 
0) 0 


From equation (4.8) we have 


Im F(E’ ) AE 
F(z) = 4.13 
2) ic E’—z Vocases a 
0) 
Now let Im z > 0 so that z > E. Since 
] 
lim ———— = 7 —E 4.14 
eb hee BSE one) 
where ¥ represents the principal part, we have 
Im < 1 d. 
Re F(E) = to [ er E’ — : (4.15) 
~ JE +E 
Since 21g) = — [{ V dr is real, the expression given by equation (4.15) may 
be written in the form 
Im g(E’) dE’ 
Reg=q)+- Lo | meee DE (4.16) 
J 


We can now use the optical theorem (2.13) to express Imq in terms of the 
scattering cross-section: 


Img =— 
mq =7 9 


so that 


Re q(E) = qo + 73 io [EE > (4.17) 


ae 


It is interesting to consider the limiting case of equation (4.17) when E = 0. 
Since the amplitude g at zero energy is related to the scattering length by 
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q(0) = —a, it follows from equation (4.17) that when E = 0: 


1 | d, 
a=ay—- =>7|]odk+) — (4.18) 
“ane / a 


The importance of the dispersion relations in collision theory lies in that 
they can be used to derive rigorous formulae for quantities directly related 
to measured values, and can therefore serve as a criterion for the analysis 
of experimental data (see § 16). 


Chapter 2 


APPROXIMATE METHODS FOR THE DETERMINATION OF 


WAVE FUNCTIONS AND SCATTERING AMPLITUDES 


5. PERTURBATION THEORY 


IN § 2 we estimated the effect of the field of force on the wave function near 
the centre of force for k = 0, and saw that it could be characterized by the 


quantity { Vr dr. 


co O 


If i [Vir dr < 1, this effect is small for k = 0 and even more so when 


0 
k # 0, so that V can be regarded as a small perturbation. Let V = «Vy as 


in equation (3.12) and expand the solution of (2.8) into the series 
y= Dey 

Let us take a plane wave for pO 
yy) a gikvr 

The expression for y™ is then 


gikir— r’| 


yO =- Tory ov YOY) dr’ 


The scattering amplitude q will also be expanded into a power Series: 


q = yy gman! 
It follows from equation (2.10) that 


1 | 
gq” a es e@ ikver Vow” dr’ 
2n 
The first term in the expansion (5.4), namely 


gq” ss = | Vo eo aN dr 
2i 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


represents the well-known Born approximation. If Vo is independent of 
the angles, we can carry out integration with respect to the angles in 


equation (5.6). The result is 
oC 
gre? |, ee 2 dp 
Kr 


0 
29 


(5.7) 
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The quantity 
86 
K =2k sins (5.8) 


is the momentum (in atomic units) transferred to the centre of force as a 
result of scattering through an angle 0, where @ is the angle between v and w’. 
From the formula 
sin Kr (kr 
MOAT =F [4n(2t-+ ]Y, (0) 
r 


and using equation (2.26), we obtain the following expression for the partial 
amplitude on the Born approximation: 


(5.9) 


r2 


gq)? =—2| Vou? dr (5.10) 
0 


It is evident from the expression for q, given by equation (2.47), and from 
equations (2.45) and (2.46), that the partial amplitude q, 1s a meromorphic 
function of the interaction constant «. Consequently, the expansion given 
by (5.4) (and the equivalent expansion for the partial amplitudes) is only 
valid for values of a smaller than the radius of convergence (a detailed 
investigation of the convergence of the series for g in powers of & is given 
in [13]). 

However, the expression for the partial amplitude can be written in the 
form of a ratio of two series of powers of a [14], which is valid for all a. 
With this purpose in mind, let 


iS 2 atx () (5.11) 
in (2.40) and assume that 
x(0) = Ul, (5.12) 


Then, 
x) = 2k(-1)'7! fu (ur’)—u_ y(n u(r) ]x!"" Vo dr’ (5.13) 
h 


It follows from the general properties of Volterra’s integral equations that 
the series given by (5.11) converges for all a. 

If we substitute equation (5.11) into (2.44), we obtain the following 
series for the Jost function 


L, = y an Ey” (S: 14) 
where 
LO) = | (5.15) 
[9 = 2(- L)'k | Vouju_,-,dr42ik Vou dr (5.16) 
0 0 


UP = 2(-1)'k ui xfPu_,-, dr + 2ik } Vou xt? dr (5.17) 
; ; 


QO 
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Finally, substituting equation (5.14) into (2.47) we obtain 
LY ot Im Ly 
a OL 
The approximate expression for g, can be obtained by cutting off the series 


in the numerator and denominator. In particular, if we retain first-order 
terms only, we have, in view of equation (5.15), 


(5.18) 


a / Im Lh) 
pe Mere (5.19 
ss k ( 4 
or on substituting the explicit expression given by (5.16) 
~2] Vuj dr 
q, = (5.20) 


1+2k | Vuuju yy dr-+2ik | Vu? dr 
0 0 


Comparison of (5.20) and (5.10) shows that the Born approximation can 
be obtained from (5.20) if the integrals in the denominator are neglected 
in comparison with unity. 

As an example of the application of Born’s approximation consider 
the scattering of a charged particle by a distributed charge of density p 
where 


i pdr=Z (5.21) 
In this case, the amplitude g® can be directly expressed in terms of p. 


Suppose that Z, is the charge on the particle. Its potential energy V must 
then satisfy the equation 


AV =—42Zy)p (5.22) 
Moreover, we can write 
A eiker = — K2 eik.r (5.23) 


On multiplying equation (5.22) by exp (iK.r) and (5.23) by V, we have on 
subtracting and integrating over all space 


i [el AV —VA (e!K")] dr = K? { el Vdr—4nZy{e™* pdr (5.24) 


According to Green’s theorem, the integral on the left-hand side vanishes 
if /” has no singularities and falls off sufficiently rapidly as r > oo. If this 
iS SO, 


. 4 . 
1 Vdr=-— | e'™" odr (5.25) 
so that 


22 ; 
g™ La <A |e pdr (5.26) 


32 THE THEORY OF FLECTRON-ATOM COLLISIONS 


The quantity 


1 
F(K) = =| er odr (5.27) 
is known as the form factor. It can be used to rewrite equation (5.26) as: 
2252 
gi = Saar F(K) (5.28) 


If V(r) has a Coulomb singularity, this singularity must be considered 

separately. For example, when the scattering system is an atom of charge Z, 

we have 4 
Ligh 

Voip (5.29) 


where W is finite at r = 0. We then have 


2Z 
gi = ra [Z, + ZF(K)] (5.30) 
For a neutral atom, 
Z = —Z,;, ($.31) 
so that the formula given by equation (5.30) becomes 
=e 21 
gO = 2 [F(K)-1] (5.32) 


6. THE METHOD OF TRIAL POTENTIALS 


When the potential energy V cannot be regardcd as a small perturbation, 
one can try to choose an expression for VY, such that there is a rcasonably 
simple solution for scattering in this ficld, and at the same time V—V can 
be regarded as a small perturbation [15]. If we confine our attention to 
simple trial expressions for Vg which contain only one or two parameters, 
we cannot, of course, expect to obtain a good approximation to the exact 
wave function and the exact amplitude over a wide range of encrgies and 
angles 0. 

However, the situation is more favourable when &R <s 1, where R is 
the range of the ficld V and w is a slowly varying function of A and 0. Having 
chosen V, so as to achieve the best approximation to W at zcro energy, we 
can use the same expression for V, to obtain an adequate approximation 
in the entire intcrvalO < kR < 1. 

It 1s natural to choose V so that there are no effects of the first order 
in V—Vo. 

Consider the limiting case of zero energy. The wave functions w and 
Wo corresponding to V and Vg respectivcly satisfy the equations 


Ay —2Vy =0 (6.1) 
Ao —2Vo Wo = 9 (6.2) 
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and have an asymptotic behaviour of the form 


q 


w=1+ rs (6.3) 
jae 7 (6.4) 


Let us multiply equation (6.1) by Wo, (6.2) by w and then subtract and integrate 
over all space. In view of (6.3) and (6.4) we have 


2n(q—4o) =| (V—Vobow dr (6.5) 

We shall require that g = qo, so that from equation (6.5): 
| (V—V,)Wowdr =0 (6.6) 

On a first approximation we can replace w by Wo, so that 
| (V—V,)W2dr =0 (6.7) 


which is the condition for first-order effects to vanish in perturbation theory 
(see, for example, [17] Section 10). 

Consider now the condition which must be satisfied by the difference 
OV = V—V, for it to be looked upon as a small perturbation. Let 


by =Y—Wo (6.8) 
The equation for 6 follows from (6.1) and (6.2) and Is 
Ad —2V, dy = 26V (6.9) 
Furthermore, consider the irregular solution of the equation 
Ag—2V,9 =0 (6.10) 


which behaves as 1/r when r— 0. If we multiply equation (6.9) by @, 
(6.10) by é6W and then subtract and integrate over all space, we obtain, taking 
into account the singularity in g at r > 0, 


i 
5Y(0) = = { 5Vow dr (6.11) 


The order of magnitude will be retained if g is replaced by 1/r in the integrand. 
Hence, 


|5Y(0)| <(max |])2 | |5V|r dr (6.12) 
0 
The quantity |dy(0)|/max || may be looked upon as a measure of the 
relative departure of wy from Wo, within the sphere of action of the field of 
force. It is evident from (6.12) that the condition for this departure to be 
small is 


2|\sv\rdr< 1 (6.13) 
0 
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When V, = 0 this condition becomes identical with the condition for the 
applicability of perturbation theory for V (§ 5), as should be the case. 

If the trial expression for Vj contains one parameter, equation (6.7) can 
be used to define it. If equation (6.13) is then not satisfied, it is necessary 
to introduce additional parameters in order to reduce the magnitude of 


co.8) 


{ }OV|r dr. An example illustrating the possibilities of this method will be 


O 
considered in § 13. 


7. THE VARIATIONAL METHOD 


The problem involving the scattering of a particle by a centre of force 
may be formulated as a variational problem involving the determination 
of the stationary value or an extremum of a functional. This formulation 
is used to establish certain general relationships such as, for example, the 
virial theorem [16], and also in the development of the variational method 
for approximate calculations. 

In the variational method one chooses a suitable analytical expression 
for the continuous-spectrum wave function, i.e. the so-called trial function 
which contains a number of parameters. The parameters are determined 
by requiring that the functional which includes the wave function should 
exhibit either an extremum or a Stationary value. 

Apart from the latter condition, there are a number of integral relations 
which the wave function must satisfy. The exact function satisfies all the 
conditions whereas the trial function may satisfy only a few of them. The 
various modifications of the variational method differ from each other by 
the number of integral relations satisfied by the trial function. 

In this discussion we shall confine our attention to Kohn’s method and 
the interesting consequences recently established on this method in the case 
of zero energy. 

Consider the functional 


P= { o*(v>,r)(H—E) W(v,, r)dr (7.1) 


where w and ¢ have the asymptotic forms 
ikr 
*4: e 
Wy a gikver ae ee (7.2) 


~ikr 


e 


(7.3) 


(p ae pikver 4 q' 


and g and q’ are not necessarily the exact amplitudes. The functional 
vanishes for functions w which are exact solutions of the equation 


(H—E)y =0 (7.4) 


We shall vary the functions w, » so that in the asymptotic expression q, q’ 
change by 6g and dq’ respectively. The details of the calculations are given 
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in {17]. We shall not reproduce them here but merely state the final result: 
6Lf = 27 dq(v,, V2) (7.5) 
It follows that the expression 


Q=- = L£+q (7.6) 
2m 
is stationary with respect to variations in the wave function, and that its 
value is equal to the amplitude g when exact wave functions are employed. 
Suppose now that instead of the exact wave function W we use a trial 
function w, containing a few adjustable parameters and having the asymp- 
totic form 
ikr 


WV, oe evr ta, (7.7) 


Let the value of the functional (7.1) calculated with the aid of this function be 
denoted by Y,. 

The approximate expression for the scattering amplitude in Kohn’s 
method is taken to be the stationary value of (7.6) obtained using the function 


Wy, 1.€. 


an=st(= 5 ¥,+41) (7.8) 
Tl 
In particular, if Y, includes the parameters c,, C2, ..., c,, the conditions that 
SY, should be stationary with respect to c,, C2, ..., Cay G are [in view of (7.5)] 
OL 
—=0 (7.9) 
dc; 
OL 
—=2n (7.10) 
0”, 


In the above discussion, we have only used the fact that the functional QO 
is stationary. The problem as to whether (7.8) yields the upper or lower 
bounds remains unresolved. In order to determine the nature of the 
approximation which can be achieved by the variational method, consider 
[18] the phase of the s-wave. Let y = f(r)N 4x r. The function f(r) will be 
normalized so that it has the asymptotic form 


fr 7sin kr+dAcoskr (7.11) 
in order to avoid complex quantities. If yw is the exact wave function, 
i = (tan d)/k. If, on the other hand, w is the trial function, then according 
to Kohn 

1 i? @) 
ctané = 4-2] f(H—E)fdr (7.12) 
0 


where we must take the stationary value of the functional. 
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Let us take the trial function in the form 


f=9 + 2 Cid (7.13) 
where 
1 
O= re sinkr+Aq,(r)cos kr (7.14) 


and @,(r) is a function which vanishes at r = 0 and is equal to unity when 
r — oo, but is arbitrary in all other respects; y; are functions which fall off 
rapidly as r — oo and ¢; are the constants to be varied. We then have 


J 

ptand=A+ ) cie;Mij +2) CR +B (7.15) 
i,j i 
where 
Mi = My =2| x(E-H)y,dr (7.16) 
0) 
R, = 2 | x{E-H)g dr (7.17) 
O 

B=2[ 9(E-Hodr (7.18) 


0 
Variation of the constants c; yields 


» Mic; =—R, (7.19) 
1 


or in matrix form 
MC =-R (7.20) 


The matrix M~* may turn out to be singular when a determination is 
made of & and the other parameters entering into y;. The solution given 
by (7.20) will then not exist. Near these values the phase expressed as a 
function of A undergoes rapid changes and may be either much greater 
or much smaller than the true value. 

Specific calculations in a particular example [18] (see also § 16) show 
that as the number of terms in equation (7.13) is increased, there is an 
increase in the number of singularities in M~', but that the phase fluctuations 
near these points are gradually smoothed out. 

However the limit of zero energy is a special case. It turns out that 
Kohn’s expression for the amplitude [equation (7.8)] gives the upper bound 
for the scattering length. The case when there is no bound state of the particle 
in the central force field and the case when the bound state does exist must 
be considered separately. 
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Let w, be the trial function at zero energy which, in view of equation (7.7), 
has the asymptotic form 
Ay 


on ees (7.21) 


r 


and consider the difference between w, and the exact function w: 


g= W.-W (7.22) 
We shall use the following identity which can be verified directly 
1 I 
a=a,+ | dtthdr = | ota dr (7.23) 
2n 271 


where a is the exact scattering length. The quantity 
1 
apy =a,t+ | W,Hw,dr (7.24) 


is, according to equation (7.8), the Kohn scattering length if y, is taken to 
be the function giving the stationary value of ax. Consequently, 


I 
a =ag—5- | gig ar (7.25) 
Consider the integral 
j= | gHg dr R (7.26) 


In the absence of bound states, H is a positive definite operator in the 
space of square-integrable functions. However, f is not square-integrable. 
We therefore introduce the integral 


J(A) =| 9 e~""Hge *dr (7.27) 


which includes the square-integrable function g exp(—Ar); for 4 = 0 this 
becomes identical with (7.26). The integral J(A) is such that J(4) > 0. We 
shall show that J remains positive when 2 > 0. To prove this, it will be 
sufficient to show that J(/) is continuous at 2 = 0. 

Consider the difference 


J) —J(0) = | g(e7?*"—1)Hg dr +22 | g? eo 2** dr—2) i gVge72** dr (7.28) 


Since {g(r)| has an upper limit for all r and Vg > 0 when r > ©, it can be 
shown that each of the three terms on the right-hand side of (7.28) becomes 
as small as desired when 24-0; this proves the continuity of J(A). 
Therefore the integral given by (7.26) is positive, or in the limiting case zero. 
It then follows from (7.25) that 


agdg (7.29) 
i.e. the Kohn scattering length is the upper-limit approximation for the 
scattering length. 
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Consider now the case when there are bound states. We shall assume 
for simplicity that there is only one such state. Let 


Hy, =| viHy, dr i,k = 1.2 (7.30) 


be the matrix elements of the energy operator H calculated with the aid 
of the orthogonal functions Y, and w2, where w, is an adequate trial function 
for the single bound state, so that H,, <0. 


Let us reduce the matrix 
a | 
H,, H4, 


(; 0 

0 ES 

One of the eigenvalues, £,, will be negative and the other, £3, positive 
or zero, so that 


to the diagonal form 


E,E, <0 (7.31) 


Since the product £, F, of the eigenvalues is equal to the determinant of 
the matrix Ge ag it follows that 
Hy, #2 

. H,,H..—Hi, <0 (7.32) 
We note that since y, enters twice into each of the terms in the inequality, 
the normalizing factor of Y, is unimportant. The inequality (7.32) remains 
valid when w, is replaced by another function w, which is not orthogonal 
to w,. This may be verified by substituting 


Yo =Ui-W [Wide (7.33) 
The inequality (7.32) then becomes 
H,,H5.-H <0 (7.34) 


where the dashed matrix elements represent the replacement of w, by >. 
If we now take a suitable bound-state function for wy, and let 


= —E, (7.35) 


where E, is the value of #/,, calculated using w,: and if we replace wy, by 
g, we have instead of (7.34) 


5 


] ~ 
| otigar <~—( [na dr] (7.36) 


t 


Since Hg = H(y,—W) = Aw, (because JI = 0) we may write 


- | atigar <=( | vi ttidr) (7.37) 
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Then (7.25) becomes 


2 
a<a,+ =. (| Ww Hy, ar) (7.38) 


Therefore the right-hand side of (7.38) is the upper bound for the scattering 
length. 

It is important to note that (7.38) includes the trial function and the 
variational energy but it does not include the exact bound-state wave 
function and the exact value of the binding energy. 

The inequality given by (7.38) may be rewritten in the form 


a< ax (7.39) 
where ax is Kohn’s scattering length calculated using the function 
J 
Vim War | Wi Hbeae (7.40) 
{ 


The orthogonality condition [ w,W; dr = 0 is satisfied. It is evident from 
this that in Kohn’s method the trial function can reasonably be taken in 
a form containing the variational function of the bound state in addition 
to other expressions. 


8. THE METHOD OF TRIAL FUNCTIONS 


In the variational method the optimum parameters of the trial function 
are determined from the condition that the functional should be an extremum. 
However there are other methods for determining the optimum values of 
the trial-function parameters. One method is based on the use of certain 
relationships which are satisfied by the exact wave function [20, 21]. Depend- 
ing on which particular relationships are employed in the determination 
of the parameters, it 1s possible to modify the method in various ways. 
In the simplest case the trial function contains only one parameter. According 
to [20], the parameter may be determined from the integral equation for the 
scattering amplitude given by equation (2.25). It is convenient to choose 
the parameter so that it plays the role of the amplitude in the asymptotic 
expression for the trial function 


1 l i nem 
f~tsin(kr ~ 2) + qe z) (8.1) 


The equation for q, is then of the form 


q=—2 { u Vf, dr (8.2) 
0 
If the trial function F, depends on the m parameters c,, c,...., Cc, 1n addition 
to g, and has the asymptotic form 
l l i as 
fi~ ,, Sit (11 _ =| + q,e ( 2) (8.3) 
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then one of the possible methods of determining these parameters is as 
follows. Substitute /(c,, ¢2,...,¢C,,9, 7) Into the right-hand side of the 
integral equation (2.23). This yields the function 


GC 15 Cas 0063 Coe QP) = u(r) +2 { fo(nu(r)—ulryo(r) Vr f(r) dr’ (8.4) 
0 
Since f; is not an exact solution of (2.23), it follows that g, # ff. We expand 
g, and f; into series in powers of r and require that the first 7 terms be identical. 
This yields a system of equations for the coefficients c,. In order to determine 
the amplitude g,, we substitute /, into the integral relation for the ampli- 
tude (8.2). 

The application of one of the variants of the method of trial functions 
to the scattering of electrons by atoms has been described in a number of 
papers [21]. 

The disadvantage of the method of trial functions is the absence of a 
criterion of the type given by (7.29) which would determine the nature of the 
approximation in the limiting case of zero energy. Moreover, determination 
of the coefficients c, by means of the expansion of the trial function f(r) in 
powers of r becomes ineffective when the potential well is so deep that the 
wave function has a number of zeros inside the well. 


Chapter 3 


COLLISIONS OF ELECTRONS WITH HYDROGEN ATOMS: 
GENERAL THEORY 


9. BASIC RELATIONSHIPS 


THE collision of an electron with a hydrogen atom is one of the simplest 
processes, particularly when the energy of the incident electron is insufficient 
to excite the atom and the scattering 1s elastic. Nevertheless, even the 
description of elastic scattering of electrons by hydrogen atoms 1s theoretically 
quite a difficult problem. 

We now consider the main factors which characterize the scattering of 
a slow electron by an atom. Suppose for the sake of simplicity that the 
electron shell of the atom is “‘frozen”’ so that it does not exhibit any changes 
during the interaction. The effect of the atom on an electron can then be 
determined by considering an equivalent centre of force. The interaction 
between the electron and the atom is in fact found to be attractive in 
character. 

It sometimes happens that when the atom 1s regarded as a centre of 
force there exists a shallow zero energy level for the electron, i.e. a negative 
ion can be formed, or that there is a virtual level. As a result, low-energy 
scattering approaches resonance scattering. Even this single fact complicates 
quantitative discussion of the problem quite considerably (for example, 
it prevents us from using the perturbation theory). 

When the electron shel! of the atom is “‘unfrozen’’ a new factor arises, 
namely, polarization of the atom by the scattered electron. At large distances 
between the electron and the atom the influence of polarization may be 
described by an interaction of the form U = —f/r*, where f is the polariza- 
bility of the atom. This expression represents the interaction between the 
point charge and the induced dipole. However, at small distances, a single 
coefficient f is insufficient to describe the polarization effect. In contrast 
to the static polarization of an atom in a homogeneous field, we must now 
deal with the more complicated dynamic polarization. Moreover, the 
indistinguishability of the electrons has a major effect on the various detailed 
features of the scattering process. As a result, the so-called exchange- 
scattering, in which the incident electron changes places with the atomic 
electrons, becomes indistinguishable from simple scattering and both must 
be regarded as a single process. Mathematically, the indistinguishability 
of the electrons is reflected in the symmetry properties of the wave function. 
When these symmetry properties are taken into account, it is found that the 
effective interaction between the electron and the atom becomes non-local, 

41 


42 THE THEORY OF ELECTRON-ATOM COLLISIONS 


and the equations for the wave function F(r) of the scattered electron include 
an exchange operator of the form 


WF = i W(r',r) Flr’) dr’ 


in addition to the central-force term VF. 

When excitation is possible the theory becomes even more complicated. 
In elastic scattering calculations it is necessary to solve only one equation 
when polarization is neglected. When excitation is possible the corresponding 
calculation (on the same approximation) requires the solution of at least 
two equations even when there is only one inelastic scattering channel. 
When there are several channels it fs necessary to consider a complicated 
set of equations, since the cross-sections for processes in the different channels 
are related, in view of the conservation of the number of particles, which 
is expressed mathematically by the unitary property. It follows that the 
elastic scattering 1s also modified under the influence of inelastic scattering. 

The effect of all these factors on the collision between an electron and 
atom, their relative importance, and the connection between them, have 
not as yet been adequately explored, although it is possible to write down 
the equations which formally represent the exact theory including all the 
above factors. 

We now consider the theory of time-dependent states, which is analogous 
to the theory in § 1.t The energy operator is of the form 

H = Ho(ry)+Ho(r2)+ Volri) + Volr2) + VCr1, r2) (9.1) 

where H,(r,) and H)(r,) are kinetic energy operators for the first and 
second electrons respectively, V)(r,) and V,(r,) are the energies of inter- 
action between the electrons and the nucleus, and V(r,, r,) is the energy of 
interaction between the electrons: 


Vir)=-—,  Volr=-—,  Vur=——) 2) 
ry r2 lr = >| 

The wave function for the electron + atom system depends on the coordi- 
nates and the spin variables of the electrons. In the non-relativistic theory 
the total spin of the system is a constant of motion, so that we can separate 
out the spin variables and confine our attention to the coordinate part of 
the wave function. Since the total wave function is antisymmetric with 
respect to the interchange of coordinates and spins, there is a connection 
between the symmetry properties of the coordinate part of the wave function 
and the magnitude of the total spin of the system. 

In the case of the two-electron system which we are considering here, 
there are two possible spin states, namely, a triplet and a singlet state. In 
the triplet state the coordinate function is antisymmetric while in the singlet 
state it 1s Symmetric. 


+ The ensuing analysis is a variant of the so-called theory of rearrangement collisions 
for two particles in the field of a fixed centre of force. It will be assumed that the two 
particles do not form a bound state. The various problems in the theory of rearrangement 
collisions are discussed in [22-25]. 
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Let M(r,, 12, ¢) be the coordinate part of the wave function satisfying 

the equation 

oP 

is = H®D 
As has been pointed out above, the function ® should be symmetric or 
antisymmetric with respect to the interchange of r, and r,. However, it is 
in many ways more convenient to begin without assigning symmetry 
properties to ® and to symmetrize the result in the final stages of the 
analysis. 

Proceeding as in § 1, we shall express ® on the interaction representation, 
that is, we expand ® in terms of the eigenfunctions for a non-interacting 
electron-atom system. Two such systems can be formed. One 1s described 
by the operator 

H, = Ho(ry)+ Vo(r1)+ Hora) (9.3) 
and the eigenfunction 
QL (ry,r2) = k*(2n)7F el (r 1) (9.4) 
where w, 1s the atomic function. The subscript a represents the set of 
quantum numbers which characterize the state of the atomic electron, 1.e. 
the energy £,, the angular momentum, and the component of the angular 
momentum in a special direction (in the present case, the direction of 
motion of the incident electron). The plane wave exp (ékv.r,) describes 
the free motion of the electron with energy E = 4k? in the direction of v. 
The functions g{'? are normalized to 6(E~E’)d(v—v’)d,¢ and _ satisfy 
the equation 
Hy, py = &gy” (9.5) 
where 
S=4k*tE, (9.6) 
The eigenvalues & of the operator H, form the continuous spectrum 
beginning with the negative quantity —|Eo|, which is equal to the ground- 
state energy of the hydrogen atom. 
The other system is described by the operator 


H, = Ho(ry)+ Holr2)+ Volra) (9.7) 

and the eigenfunction 
QLO(ry, 2) = kF(2m) 77 eh (12) (9.8) 
whose significance 1s analogous to equation (9.4), except that r, has been 
replaced by r,. Either the functions given by (9.4) or those given by (9.8) 


can be used as the basic functions. Accordingly, we shall consider the 
following two equivalent expansions for D(7,, r, f): 


O(r,,°2,0 => [¥O, vy, te QO...) dé dO (9.9) 


P(r,,r2,0 =) | POE ye oP r,, 1) dé dQ (9.10) 


a 
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As in § 1, we shall assume the existence of the limits 
yO = lim PO, v, t) P= 1,2 (9.11) 


t~+0 


and determine the S-matrix from the relation 
POW =VS | S*(6,v, 0, 6", v", «P(E, w’)d6' dQ’ i,i’=1,2 (9.12) 


Let us express ® in the Heisenberg representation. Let o(r,,r2, &, v) 
be the eigenfunction of the operator H corresponding to the eigenvalue @, 
so that 


Hp =€o (9.13) 


where v is a unit vector which enters into the argument of @ (as in § 1). 
The eigenvalues & of the operator H form a continuous spectrum lying 
in the same region as the continuous spectrum of H, and H,, that is, 
—|E,| < & < o. Moreover, the region & < —|Ep5| includes a discrete level 
of the negative hydrogen ion, but this level is of no significance in the present 
context. 

In order to determine the function g, we must also specify its asymptotic 
behaviour. We shall suppose that an electron of energy 3X¢ is incident on 
an atom in the state a, the position vectors of the incident and atomic 
electrons being respectively r, and r,. Accordingly, the function @ will 
asymptotically tend to g{)(r,,r,) aS r2 > oc. Therefore the function ¢ 
depends on the parameters defining the initial state of the atom, and this 
will be implied in the ensuing analysis. 

Let us write ® in the form 


O(r,,r,t) = { eS WE vy) o(r, 12,8, v)de dO (9.14) 


The relation between ‘¥®?) and can be established by comparing equation 
(9.14) with (9.9) or (9.10), and is given by 


POE, y, t) =| eWE—EM WOE y BY w)Y(E', wv) db" dQ’ (9.15) 


where 
WE WE V) =[OP rr ENO aw )dri dry  i=1,2 (9.16) 
Let us transform equation (9.16) by substituting 
P(1 1, 12,8", ¥) = Pag +X (7) 

where, in accordance with (9.6), 

6’ = tk>t+E,, 
We first consider the case i = 1. Since 

(H—&')o =0 (9.18) 
and 


(HI, — 8)? = ( (9.19) 
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we obtain 
(H —&')y = —(H — H 19%,” (9.20) 
or using the explicit form of H and H, 
(H—&')y = —[Vo(r2)+ V4, rr) | op (9.21) 


By analogy with the method used in § 1, let us define the resolvents 
G=(H-6)"!, G,=(H,-6)7', i=1,2,... 
as the limits 
G(6 + i0) = lim (H—€@-—ie)™! 


r ai (9.22) 
G(& + i0) = lim (H;-—6—ie)7! 
e> +0 
With the aid of equation (9.22) we have from (9.21) 
4 = — lim G(6' + ie) [Vo(r2)+ V1, 72) Je) (9.23) 
e+>+0 
On substituting this into (9.17), 
gp = op) — lim G(6' + ie) [Vo(r2) + V(ry, 62) Jo) (9.24) 
e+>+0 


Again, following the method of Chapter 1, the latter relation may be written 
in the form 


@ = — lim ieG(6' + ie)p@? (9.25) 
cc +0 
The equation for the resolvent 
G = G,—G,[Vo(r2)+ Vr. 2) JG (9.26) 
can now be used to transform (9.24) so that it reads 
1 
P= Poy) — H,—é'—i0 LVo(r2)+ Vir. 12) ]@ (9.27) 


This result can be used to reduce the expression for w{ to the form 
TO(8,v, &', vo) 


WCB, ¥, 67, 90) = HE — 8) 509 ¥9) Saag = 


(9.28) 


where 
TUE, v, 6", Vo) 

=| OC aS, v)LVo(ra)+V(ri,r2)] O01, 42,6, %o) dr, dr, (9.29) 
Substituting equation (9.28) into (9.15) and repeating the steps analogous 
to those given in § 1 we obtain 

POE, v) = VE, v)ba2 


(9.30) 
WE, ¥) = POE, v)—2nf | TLE. ¥, 6, v0) PAE, Vo) dQ’ 


The S-matrix clement is therefore given by 
S''(8,v, a, 8, ¥9,%) = 5(6’ —&)[5(V6 — v) 6,,,— 271 T,X(E,V,8,%)] (9.31) 
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The presence of the factor 6(é’—&) in (9.31) ensures the conservation of 
energy (6 = 6’). In accordance with equation (9.6), the energy of the final 
state is © = 4k7+£,, and therefore 


Lk2 +E, =442+E,, (9.32) 


Consider now the case / = 2 and suppose that the function q 1s, as before, 
defined by equation (9.24). The vector W%? describes the so-called exchange 
scattering. It is now convenient to express the resolvent. G in (9.25) in terms 
of the resolvent G, as in (9.26) rather than in terms of G,: 


G= G, =a G2[Vo(r1)+ Vir, r,)|G (9.33) 
Substituting (9.33) into (9.25) we obtain 
gy =— lim ieG, o@ — lim G,[Vo(r7,)+ V(rq, r2) Je (9.34) 
e>+0 E> +0 
We shall show that 
lim 6G,(& + ie)e’ = 0 (9.35) 


To do this let us expand @{0(r,, r2) in terms of go: 


gM=V~ i a(8')o? dé’ dQ’ (9.36) 


If the atomic function w,(r,) entering into y refers to the discrete spectrum, 
then a,(€") is a continuous function of é&’. But if the atomic function P(r) 
refers to the continuous spectrum, then the coefficient a,(é’) has a 6-function 
type of singularity. However for this part of the spectrum the symbol 
)’, Tepresents integration over the continuous spectrum with respect to «, 
which removes the singularity and leads to a finite expression. Let us operate 
with G,(€+ie) = (H,—E€—ie)~' on both sides of equation (9.36). The 
result is 


— ~ 92 de" dQ’ 


G,Q) = = Efe (O)% 
The integral 
a&’ 
J = | ————_ 12’ al&')o? 
i Pe |: alé')p, 


has a finite value, so that lim eJ = 0 (we note that if we evaluate GC,0~2? 
e-0 


instead of GS, the result becomes G,y2 =—g/ie, so that 
lim ieG,y® = g&). In view of equation (9.35) we have trom (9.34) 
e> +O 
\ 
Q=- ie a aye es Viry.r2)|@ (9.37) 


It follows that g is the solution of the nonhomogeneous equation with the 
resolvent G, and the solution of the homogeneous equation with the 
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resolvent G,.¢ Substitution of equation (9.37) into (9.16) yields 


T2(8, v, 8’, Vo) 


218 v,8',v9) = —- — 9.38 
wi (8,Vv,E', Vo) ue (9.38) 


where 
T2218, ¥; GE", Vo) 


=| GO ra EMO AVC 2] 91 12.6%, V0) dri dry (9.39) 


Let us now substitute (9.38) into (9.15) and take t to +00. Since according 
to equation (9.30) 
Y(E.v) = PO (6. y) 


we have for discrete w and a 
P(8,v) =0 


9.40 
PONE, v) = —2ni T.2(6, v, 6, Wo) PO) (6, v9) dQ ee) 


aoz — 20 


For continuous «, corresponding to integration over the continuous spectrum 
of the atomic states in (9.10), the element 7{2? of the 7-matrix regarded 
as a function of « has a singularity of the form 6(E,—4k$). Consequently, 
for continuous a, ‘P2(&, v) will also have a singularity. We shall not 
consider the nature of the singularity here. 

Substituting equation (9.40) into (9.12) we obtain the following expression 
for the S-matrix element: 


S21(6,v,a, 6", Vso) = — 22i15(6’ — €) T2218, v, 6", Vo) (9.41) 


aoa 


In order to introduce the symmetry properties of the wave functions let us 
construct the symmetric and antisymmetric combinations 


I 


Y= — (WE) + pr) (9.42) 
no 
| 
Yo= ars i) (9.43) 
\ 


In view of equations (9.9) and (9.10) it is evident that the state vectors given 
by the latter two expressions correspond to the symmetric and antisymmetric 
functions ®. 

Transitions from the initial to the final state for ‘¥, and ‘V, are described 
by the symmetrized S-matrix which can be constructed from S'! and S?! 
by analogy with (9.42) and (9.43): 


St =(S'!+ $71) (9.44) 
S7 = (S'!— $71) (9.45) 


+ This point is discussed on the basis of the time-independent theory in (26). 
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From (9.41) and (9.31) we have 
S*(E,v, 0,6", Vos%q) = (6 —F’) [5(V9 — V)S gay — 271 T(E, V,E',Vo)] (9.46) 

where 

Toe tag, Be (9.47) 
The same result can be obtained in another way. Instead of the functions 
oS and g we can use as the basis in the interaction representation the 
symmetric and antisymmetric combinations (op) + a?) //2 normalized to 
6(8 —&')6(v—v’), and introduce the symmetrized function 

| 


QO, =——= lim ieG(6 + ie)(90 +. 0©) (9.48) 
\ 2 e> +0 


which may be rewritten in the form 
1 
a ae (9.49) 


where P,, is the exchange operator. Consider the integral 


l 
/2 | (Py, £055) * G4 dr, dr, 
Using equation (9.49) we have 


I j | 
5 | @@ Hoye: dr, dr, = 5 | Pi verc +P,,)9 dr, dr, (9.50) 


However, 


| (L4 Py») "(C1 +P,,)pdrydr, = | p+ P,2)°9* dr, dr, 


= 2] p(LtPy2)9* dry dry (9.51) 


zo 


and therefore 

I ) : ) 
15 | (P+ 9)" 0. dr, dr, = | (oh + 0% Pedr, dr, 
‘ 


=| ph *odr, dr, | gy * pdr, dr, (9.52) 


This corresponds to the representation of the matrix 77, in the form given 
by (9.47). It is possible to establish a relation between 75, and the effective 
cross-section by analogy with $1. But it must be remembered that the 
symmetrized wave functions correspond to the superposition of two physi- 
cally identical situations: scattering of the first electron by the atom con- 
taining the second electron, and scattcring of the second electron by the 
atom containing the first electron. The expression for the differential cross- 
section for elastic scattering is of the form 
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= (2n)* — (9.53) 


do 


ee alta ale 


while the expression for the a. cross-section for inelastic scattering 
in which the atom undergoes a oer from the state a, to the state « is 


dox, = 2dQ (9.54) 


= (2n)*5 


‘Gas 


When comparing theory with spain it must be remembered that in 
practice one 1s usually dealing with the scattering of an unpolarized beam 
of electrons forming a statistical mixture. Accordingly, equation (9.54) 
must be averaged over the triplet and singlet states. The result is 


(d6y5a> = tdoz,,+4do;,, (9.55) 
There 1s an important relation between the elastic and inelastic scattering 
cross-sections which is a consequence of the unitarity of the S-matrix and 


forms a generalization of the optical theorem given by (1.48). The matrices 
S* and S~ are unitary, i.e. they obey a relation analogous to equation (1.7): 


yy { S**(6,a,v, 6’, w’, v)S*(8, a, v, 6", 0", v") dé dQ 
= 5(6'—8") S(v"—v’) 5g" (9.56) 


With the aid of (9.46) we can use the unitary property to deduce a relation 
analogous to (1.48), namely 


—Im T,~,,(0) = 2 ‘ | \T.¢,|7 dQ * (9.57) 


In this formula the summation is carried out only over those levels which 
can be excited at a given energy in accordance with equation (9.32). When 
the initial state of the atoms is the ground state and the energy of the incident 
electron 1s insufficient for excitation of the atom, the sum in (9.57) consists 
of one tcrm only. 

The S-matrix is also invariant with respect to time reversal. This leads 
to a generalization of equation (1.8) and may be written in the form 


Sz(E,¥, 6", Vo) = SZ(6', — Vo, &, —¥) (9.58) 
Accordingly, we have the following relation for the 7-matrix elements 
TE, Vy é, Vo) = T,5(é, —Vo; &, _ v) (9.59) 


It is understood that in the time-reversed transition described by the right- 
hand sides of (9.58) and (9.59), the signs of the orbital and spin angular 
momenta have been reversed. The transition of an atom from one state 
to another which we have been considering is, in fact, invariant with respect 
to reversal of all the components of the angular and linear momenta. It 
follows that we can use equations (9.59) and (9.54) to write down the 
detailed-balance relation: 

ko dOyo = KZ doy, (9.60) 
(Sf, Vo, &, Po). 


+ 72 (0) represents Tt 
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10. EQUATIONS FOR THE WAVE FUNCTION OF THE 
SCATTERED ELECTRON 


If the eigenfunction g(r,,1r,) of the energy operator for the atom + 
electron system is expanded in terms of the atomic functions @,(r,), the set 
of expansion coefficients will represent the wave function of the scattered 
electron. 

In this section we shall establish the equations and asymptotic behaviour 
for the wave function of the scattered electron. 

To start with, let us change the normalization by multiplying g@ and @4 
by (2z)#ko*. The corresponding functions will be denoted by w and w,. 
For the & functions we have the integral equation 

Wry.) = — lim ieG(E& + ie) e'Xovo”? Wali) (10.1) 
e> +0 
As has been shown above, from equation (10.1) can be obtained two integral 
equations (one through the resolvent G, and the other through G,): 


W(r1,12) = eae Wall 1) — Gy(E+ iO) LVo(r2)+ V("q, ry) |p (10.2) 
Wry .r2) = —G(E+ 10) [Vo(ry) + V(ry, r.) |v (10.3) 
or, using Green's functions 


W(r 1,1) = ebovor? Walt'1)— 
—[GrGyra rt rd LMoOr + Vir Wri. rg) dry dry (10.4) 


e 


Wrst.) =— [Galt rari) VorDt Vr rawr dridr, (10.5) 


» 


The latter equations may be used to investigate the asymptotic behaviour 
of w(r,, r,), as follows. Consider the asymptotic behaviour of the Green’s 
functions G, and G,. These two functions are expressed in terms of the 
free-motion function Gy. By expanding G, and G, in terms of the eigen- 
functions of H, and H), it is easy to show that 


Gi(ry,r2.1},12) = ye Wir Wir )Golk,, ra.) (10.6) 
and 
Go(ry.12,1 1,02) = » Warr WE (12) Gol kas 11514) (10.7) 
where 
G(k,rr)== | — eftalr—e'| 
2n|r—r | (10.8) 


k, = /2(E,,-E,) +k3 
Whenr > r’ 
Go me = gikar eo” kane (10.9) 
ar 


For those values of x for which £,,+449 < £,. the function Go(r) falls off 
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as exp (—|k,|r). Substituting equation (10.9) into (10.6) and (10.7) and 
retaining only those terms in the sum in which A, 18 real, we have 


GD vale te) — 


T2772 Qn ka*>0 


Substituting equation (10.10) into (10.4) and on 
eikor2 eikar2 
W (ta) ~ (cas — Yad) +E Saws — aly) (101) 


r27- @ (k2> 0) r2 
where 


] —ikon.r>’ / r , ' / , ’ t , 
ive aa a see wiry )LVo(ra) + V(ri.ro) wry. 15) dry dry (10.12) 


er tenes (10.10) 


When r, > ©, 


W (11,12) ~ oy bao W (rz) (10.13) 


where 
I —ik, nr,’ / ’ rf , , / , 1 
dan = | ¢ see wir (Vor +V(r. ra) Wri. ra) | dry dr, (10.14) 


Comparison of equations (10.12), (10.14) with (9.29) and (9.39) yields the 
relation between f,,2. Za, and T{?, TS): 


T(E, 1,8, V9) = ~ On nye * Bk? f, (M1, Vo) (10.15) 
J 
Tyoa (Es, 8, Vo) = — TG KG KF Gaga, Vo) (10.16) 
° (27) 
Let us now expand w(r,, r,) in terms of the atomic functions: 
Wr1,r2) = » W(r1)F (12) (10.17) 


where the symbol ) includes integration over the continuous spectrum. 


We substitute equations (10.17) and (10.6) into (10.4), multiply by wz and 
integrate. For the coefficients F,(r,) which represent the wave function of 
the scattered electron we then have the set of equations 


F(a) =e"? Bae DY, | WEE V0rD) + V0, 75) 
Wael) Golkas Ys 05)Fa(rg) dry dry (10.18) 


where the integration is with respect to r;. Let 


Via(12) =| we(rs)[Vo(r2)+ V(ry, ro) |W.(ry) dry (10.19) 
Using the explicit form of V sa and V(r,,r,) given by (9.2), we obtain 
Via(r2) = — = bat |" Ya(riWa(ra) aU (10.20) 
Iro—ry| — r,| 
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Substituting (10.20) into (10.18), we obtain 
F,(1) = oF Saat D | Golkas st Var Fa(rdr’ (10.21) 


where the subscript 2 has been omitted. This system of integral equations is 

equivalent to the set of differential equations 
AF, +(k?-—2V,,)F, =2 > Vj Fy (10.22) 

a’ fa 
with the asymptotic condition 

Kar 

: e'* 
Fy = 0%" 6 et Seon — r— 0 (10.23) 

r 


The system given by (10.22) can also be obtained by substituting the expan- 
sion given by (10.17) into 


[¥iCDH-S Wry n)dr, =0 (10.24) 


It is shown in [17] that the latter relation 1s a consequence of the fact that 
the functional 


T=) | VEOH Obert) dr dry (10.25) 


is stationary with respect to arbitrary variations of F, which vanish when 
r, — oo. On this basis we shall suppose that equation (10.24) is valid even 
when w(r,, r.) 1s not an exact solution of (W7—é&)W = 0 but represents an 
approximate expression, for example, the sum of a finite number of terms 
in (10.17). This leads to the equations for F, which give the best approxi- 
mation in the variational sense. 

Consider the part of the sum in (10.17) which corresponds to the con- 
tinuous spectrum. It is given by the integral 


Jrssr2) = | Fe, k, ra W(x, ry) dx (10.26) 

There is also a relation between A and x which follows from the law of 
conservation of energy 

k* +? = k§+2E,, (10.27) 


It turns out that F(x, &, r) in (10.24) has a singularity at «2 = kg [27]. Let 
us consider the nature of this singularity in greater detail and expand F, 
in terms of w: 


aes >» Cay Wa! (10.28) 
Substitution of this expansion into (10.22) yields 
3 Cyn(A tke —2Vi,.) Wy ae 2 », CBA Vp YW, (10.29) 
a’ A,p Rx 


The function y,- satisfies the equation 


At2—2 (E, re -) Wy (10.30) 
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From the latter equation and from equation (10.21) we have 


1 
YY oe —2E Wy = 2 yy (5 Cpa Wi, (10.31) 
a’ B, A lr—r | pa 
Multiplying Wz, and integrating, 
1 
Calk; —2E,,) =2 yy CBA (—>) Wy (10.32) 
BA lr—r | Ba 


Substituting « = % and using (10.27), 
Cag(h) (ko— so =2 2 CRA | Wa (7, )Wo (r2)Walr )Wi(r2) Ir, — mont dr, (10.33) 
"2 


In order to eae caitlin connected with the specific asymptotic 
behaviour of continuous-spectrum Coulomb functions, we shall cut off 
the Coulomb interaction by introducing the factor exp (—mr) where @ is 
a finite positive number which can be as small as desired. The right-hand 
side of (10.33) will then be finite for kg = xk, so that 


const 
Cag) = Poe (10.34) 

Therefore, F(x, r) can be written in the form 
F(x, r) = PaaS WAT) + B(r) (10.35) 


Where A(r) is a function having no Fait at ky = K. 

When we substitute equation (10.35) into (10.26) we must determine the 
path of integration with respect to x. This can be found from the asymptotic 
form of W(r,,r2) given by equation (10.15) when r,; — oo. We substitute 
(10.35) into (10.26) and separate the angular variables. Assuming that r, 
is large, the radial part W,(k, r,) has the asymptotic form 


| 
Wir) ~ sin («0 -> ae 5 (10.36) 


so that 


lx 
ers ne sin (ns 2 + 5 
Wal?) J eo ede | 08) 


Ji(r1, 12) ~ 


can be reduced to the form 


oo 


1 ert 
J ~ const adr | a dk (10.38) 
The path of integration can now be chosen so that the pole kK = ko 1s excluded 
by taking a small circle around it in the upper half-plane, while the pole 
Ko = —kpo is excluded using a circle in the lower half-plane of complex k. 
E.A.C. S 
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The result is 


ikory 


(10.39) 


J ~ const W,,(r2) 
ry 


In the absence of a purely Coulomb field the divergence will be of the form 
(k§—K?)~1** instead of (k6—xK’)71, and instead of the asymptotic form 
(10.36) it will have the Coulomb form: 


1 l 1 
Wy =-sin er — 4 “In dur + are (I +] -*)] 
r 2 «K K 


The analysis is then very much more complicated and has not in fact been 
completed. However, the general result will be similar to that in the simpler 
case of a truncated Coulomb interaction. 

Let us now consider symmetrized function Y,. As with equation (9.49), 


; l 
this function can be obtained by applying the operator a + P,,) to 
the function w 


1 1 
Ws Se a (10.40) 


We now deduce the integral equation for w,, transforming the first compo- 
nent on the right-hand side of equation (10.40) using the inhomogeneous 
equation (10.2), whilst the second component will be transformed with the 
aid of the homogeneous equation (10.3). Since P,,G, = G, and 


Py2Vo(r2) = Volri), 
1 { t ! 
vi = ae Wall 1) -| G1(115%2511519) 
X[Vor)t+ Vir re) Wari ra)dri dr, (10.41) 
Using the expansion (10.6) for the Green’s function, 
LF es 
ba = sectoral) =F vale) | vO) 
N a 
x LVo(ra)+ Vr. 12) IGolka ta rawr, re)dry dr, (10.42) 


from which we obtain the asymptotic behaviour of w, for large r,: 


1 »tkor2 oikar. 
Y= alr) (eetn + dig) + D addled 5 (1043) 
V2 V2r, 


a2#26 Vi 2r, 


/2 
\ at n.r’> , , , , t , , wf af 
dana = — = . ey Mora) + V(r ro) Wa rg) dry dry (10.44) 


l 
Writing Ww, in the form wy, = NE (Wry, ro) + Wrz, ",)], and bearing in 


mind (10.13) and (10.15), we have 
Gaia = Sana Gaon (10.45) 
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We might expand wy in terms of the atomic functions, writing it in the form of 
(10.19). However, in this form the symmetry properties of yw, are difficult 
to see (they are ensured by that part of the expansion which refers to the 
continuous spectrum). Instead we writc yw, in a form in which the symmetry 
propcrties are immediately evident: 


i 
Wa = wr DCPs (radar Fa (ri Walr2)] (10.46) 
NV a 


This representation has an important property: it is invariant with rcspect 
to the transformation 


Ee) D Dine 


where Dj is symmetric or antisymmetric, depending on which of the two 
signs 1s taken: 

Dip = + Dja 
This result may be used to remove the singularity in F(r) in the integral 
over the continuous spectrum which is contained in (10.46) (this does not, 
however, definc FZ unambiguously). We now proceed to derive the equations 
for FZ. Substituting equation (10.46) into (10.42), multiplying by wi(r,) 
and integrating with respect to r,, we obtain 


F(a) =F Y Wala) | WEF H (rs) dry + et yay — 
LJ WED LMolrs) + Vr aPC) 
X Golkas tas ra) dry drs FY | Wry) Volrs) + Vr, 15)] 
«Vol FEC)GoCkes P25) dr} d's (10.47) 


This equation can be transformed to a more convenient form. Consider 
the following part of the last term on the right-hand side of equation (10.47): 


J= ED | Golka ta, ra Vora Welrs) deg | wr Fs(ri) dr, (10.48) 
Using equation 


a tA, Wa’ + Vo(ra)Wa = Ey Wa’ 
we can rewrite (10.48) in the form 


J = FY [GAL+EDWADGolke ras radrs [UREA )dr, (10.49) 


and hence, bearing in mind the explicit form of Gp given by (2.7), we have 
after some rearrangement 


J=£Y phalra) + [ Gk -Ex)Golkas tas ra)dra] | WIC )FECY) dr, (10.50) 


The term 
Wir )Fe(ry) dry 


+ », War) 4 
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is of exactly the same form as the first term on the right-hand side of equa- 
tion (10.47) but with the opposite sign; therefore they cancel out. It follows 
that we can reduce (10.47) to the form 


Fi (r) ie te Sao >. | Golka 12,1) Vig) Fa (13) dr5 
YL [Golkerarraddrs | Mri Fe@i dr, (10.51) 


where V,,- is given by (10.20) and IV,4. represents the kernel of the exchange 
operator 


Wray, 2) = tWi(r (Vir, r2) + Ey — tk] par) (10.52) 
On reversing the order of integration in (10.51) and substituting 


Ki(r2, ry) = Viag(h2)Go(k,, Vo, r5)— | Wari, r5)Golka; Vo,1 ry) dr; (10. 53) 
we have finally 


FE (12) = e072 5... + Y | Kiera, r2)FEC15) dr (10.54) 


The integral equations (10.51) and (10.54) are equivalent to the following 
system of integro-differential equations 


AFP 4k =2V Fo =2 | Wt(r’, r)F,(r’) dr’ 
=2y ( V_Ft+ | Ww, Ft dr’) (10.55) 


a’ #a 
with the asymptotic conditions 


thar 
Fz (r) ae elkovorr 5 ot Vinal Ms Yo) — 
In the ensuing analysis we shall, in view of equation (10.55), look upon 
the term V,,Fi + [ Wi. Fi dr’ for « #' as expressing the relation 
between the states « and a’. 
Utilizing equations (10.54) and (10.16) we have 


(10.56) 


Ti (E,n n,é, v o) = 


I 
~ (2n)? kékigi.(n, Vo) (10.57) 
Substitution of (10.57) into (9.54) yields the following expression for the 
inelastic scattering cross-section 


dot, = Z laé,|? dQ (10.58) 


The differential cross-section for elastic scattering can be obtained from 
this by substituting « = a: 


—— 7 dQ (10.59) 
These expressions can also be obtained directly from the asymptotic 
form (10.56) by dividing the flux in the solid angle dQ (i.e. k,!gi2|7 dQ), by 


the flux density in the incident wave ko. 


do pee 
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On substituting (10.57) into the unitarity condition (9.57), 
4 k 
oF im gt,(0) = ¥ —2 | laz,|2 dO (10.60) 
ko a kG 


The quantity 
7 dQ (10.61) 


k 
ot = Sok. Di | lade 
a a KO 


is the total cross-section for all the elastic and inelastic processes possible 
at the given energy. Therefore, 


4nIm qz,,(0) = ko* (10.62) 


11. EQUATIONS FOR THE PARTIAL WAVES 


Let us separate the angular variables in equations (10.55) and (10.56) 
and expand fF; in terms of the spherical harmonics as in § 2. In order to 
distinguish the quantum numbers associated with the angular momenta 
of the atomic and incident electrons, we write /,, 4, for the atomic and 
l,, ft) for the incident electrons, respectively. The symbol » represents the 
principal quantum number and the symbols + will be omitted. With all 
the indices written out in full, the expansion is 


Top 


Priya?) = > V4n(2l, +1) iP ¥,,,,0, 9) (11.1) 
I2 Pr 


The summation over jt, 1S unnecessary in view of the conservation of the 

z-component of the angular momentum of the atom + electron system: 
ly ti>y = const 

Since the z-component of the angular momentum of the incident clectron 

in the initial state is v3 = 0 it follows that n#,+. = nf. If we now sub- 

stitute equation (11.1) into (10.56), multiply by Y“ie9(0, ) and integrate, 

we obtain the following asymptotic expression for fy7%? : 


I Pee 1901 el( kar 3") 
2H2 27 gy a 2U52H2 = 
nly pty k Sin Kno? b) Onno t 1° Dei P Fn! pnts wy ¢| I .2) 
no r 
wheré 
SOlapt2 ‘ * 
Fn Ont yy = | ¥ 130 Anot rusts 21 %o> n) Nis dQ. (1 1.3) 


Next, let us substitute equation (11.1) into (10.55), multiply by Y/%.,,(0, @) 
and integrate with respect to the angles. It will also be convenient to 
substitute 


gear’ = | Yisu(O, OVarlt) Yity,(0, 9) dO (11.4) 


niypeym' ly py’ I’ p2 
and 
WM OP =r | Yaa Wart Yu y(0,9)dQ.d2’ (11.5) 


nly peg yay I2'p2 


In view of the presence of |r—r’|~* in W,,,, the form of (11.5) will be 
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different forr <r’ andrz2vr’. In order to take this into account, we shall 
use the symbol W“ for r < r’ andwW™ for r > r’. The system of equations 
(10.55) will then be of the form 
Cea + (ia - hte) jon 
dr? : na 
a 
2 LVR a SAB + | OER St de] 1) 
n’ly'ta ey’ 6 
where f,742(0) = 0 when r = 0. 

The function /,742 describes the state of both electrons on the /;, sy. 
l,, 4. Tepresentation. It is convenient to replace this by the /,./,. L, Af 
representation, where L is the total angular momentum and AZ its z-compo- 
nent (this is known as the total angular momentum representation [28)). 
Since in the nonrelativistic theory the angular momentum is a constant of 
motion, on the total angular momentum representation the system of 


equations (11.6) splits into a series of unconnected sub-systems. The transi- 
tion from one representation to the other ts carried out with the aid of 


the transformation 
LM CEM lope 
nil, ~~ » Cr ied aan (11.7) 
Hi. 2 


where C;-™,,, are the Clebsch-Gordon coefficients. Substituting equation 
(11.7) into the asymptotic expression (11.6), 


| oe on 
ra an i. sin (i, an ») =") Cie 50190 9 nno ate 


no 


5 ee 
oil ter 2 ) 


13 y Cbs 
+ » nol opOnt yyy © Ligylo2 (11.8) 
M1, H2 r 
We now divide both parts of (11.8) by C75 10190 and substitute 
1 

LM LM 
ilyly Eyly (1 1.9) 
Hlye2 C/°.,9%o ney 


On the total angular-momentum representation the function ‘amas has the 
standard asymptotic form 
4: _ tan 
at see ) 


I ox 
LM : Be as 2 LM 
Mii” gt sin Kno? = a Onn Stitg + Fol Ont yl, i (1 1.10) 
K no ~ } 
where 
LM LM —: LM isis sere 
J ngl1 nt! Cio n0100 = yy Cr ie Frit went sp (1 1.1 1) 
Bey, $2 
Using the orthogonality of the Clebsch-Gordon coefficients 
LM ClM 
2 Crayton © ly “py lo’ 2’ an Orsty fata’ Op yg? Onan’ (11.12) 


we can solve saute (11.11) to yield 


io 0 
Lop2ly pry a LM LM LM 
Gnolu8nly uy — > C76 ,0100 Folly Ciypstans (1 1.1 3) 
L 
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Consider now the transformation of the matrix elements ¥7342'242.,. 


and W7i¥2'242,.., In accordance with the general formulae for the transfor- 


mation of matrices 


Aritoty ‘le’ = 2 Cia die Cie aiseA (11.14) 
HMieawy p2 
Commencing with , we have, using the definition (10.20), 
1 1 
Via(r) are = Onn’ 4,1," + (=) (11.15) 
. Jr—r'| moyen ly ey’ 
We expand |r—r’|~' in terms of the Legendre polynomials 
lr—r'|o! = ) y,(r, r’)P, (cos a) (11.16) 
A=0 
where w 1s the angle between r and r’ and 
ps 
cre 
yArsr') = (11.17) 
r? 
nr nee ded 


We then have 


J tf / 
Via") = Onn’ Ont’ a \ Onur | YA ~')P,(cos WO) YC ” ~ ydQ 
A 
(11.18) 


where 
(Vantyn’ty’ = {Ru (r’ alr, r'yRan (rr? dr’ (11.19) 
0 


(R,1, and R,-,,, are the radial wave functions for the atomic electron.) 
Substitution of (11.18) into (11.4) yields 


J 


yolapzt2’ 2’ Soca 
t mlyyeyn'ly ey! ~ Onn’ Ont t 


5 My (Yantin'ty’ { Yasils (~) ) ea ~)P (cos a) Yin (O', p’) 
x Yiu (0', 0") dQdQ’ (11.20) 


We now use the transformation (11.14), for simplicity denoting by v the 
group of indices 

v=nl, 1, 
As a result we obtain 


I 
rLM _ + LM LM 
¥ vy 7 Oy + » (Yadntyn’ty’ > Chistous Cty’uy'ta' na" 
A 


Hits’ 22" 


. af ¥* (0, 0) Yury (8, @)P,(cosw) ¥*, (0', 9’ 
x Yin, 0", 9')dQdQ’ (11.21) 
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The expression 
h,= Loan Coe | Yiou9; P)Y1y'u,(8, )Pa(cos w)¥N (8, Q) 

ee x Ying (0%, @')dQdQ’ (11.22) 
which enters into (11.21) is independent of Af and can be expressed in 
terms of the Racah coefficient 1: 
hoi, b= (- 2 0" 444-1) wd, i, bb” 

X Cf G10 Chono(2ly + 1)F(205 + 1)7(21, + 1422, +1)? (11.23) 

Therefore, 


VAT) = — = Ow a yy (Yann Mal, l,, lh, I), L) (11.24) 
A 


For large r we have 


oe oes pee 


so that (11.24) is a multipole expansion of ¥,. The quantity /, has a 
non-zero value only under the following conditions: 
(1) The triads /,, /,, 2 and /;, /;, L satisfy the usual rules for the com- 
bination of angular momenta. 
(2) /,, 1, and /;, /; satisfy the condition expressing the conservation of 
parity, 1.e. 
(— 1 te oe (ye te (11.25) 
(3) The parameter A satisfies the inequalities 
L-i<4z<d,4+1 
Eesha Cert (11.26) 
n—| <A < (I, +15) 
which follow from the properties of Cr oie and Cais when 4 = 0 


lo = 611,'0 (11.27) 


Tol’ 


The quantity /: is tabulated in [28] for certain values of /,, /,. /;, /2, and A. 
We can transform y’ in a similar way. Starting with (10.52) and using 
(10.60) we obtain the following result 


Ww =H Y glistor th aL) 


x [5,0(Eo — SKY +7 IR CR’) (11.28) 
where 
9A, 12, L) = (— IP 2A 427 Ww A) 
x CP 9170 Chon o(2h + 1) (21, + DF(21, + 14(214, +1)? (11.29) 
and g, # 0 when 
Wy,-| <4<(,4+h): 


The quantity g, is also tabulated tn [28]. 


1-li] <4 < (1,41) 
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It follows from (11.17) that when r < r’ 


i? A 
Wry (Pi n=t) gy 50, (E, = =] 5 =ni| Ra (Rai Ar) (11.30) 
A 


Z 


whilst when r > r 
tA+]1 


Wee r) = a > Gi 5B — 3h) a é i | Rn Rn, (O) (I 1.31) 
2 P 


On the total angular momentum representation, the integro-differential 
equations (11.6) for the partial waves assume the form 


2fL oo 
ca ke - BFE st = 25 rece + [WhO DLO) ir 


dr? 
(11.32) 


Henceforward we shall omit the superscript L whenever this leads to 
ambiguities. 

The possible values of /, for different 7,, /, and Z are given in Table 1 
for (—1)''*"2 = (— 1)" and in Table 2 for (—1)"'"*"2 = (—1)*". 


TABLE | 
ny3 
L 1; 0 2; 0 22" 3; 0 3; | So2 
0 0 0 l 0 | 2 
l l | 0; 2 l 0;2 1; 3 
Z 2 2 1; 3 2 1;3 0;2:4 
3 3 3 2;4 3 2a b2 325 
TABLE 2 
mil 
L 2 3; 1 
0 
l l 1 
2 2 2 
3 1; 3 1; 3 


The system of integro-differential equations (11.32) can be transformed 
into a system of integral equations: 


fo = Cryetts) + Of Male for) de' +S | Nurs rr) dr! (11.33) 
vy’ 0 v’ 0 
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where ¢,,,, !@,, and N,,- are defined by 
M,y(r, r’) = 2[v,(r)u(r’) aie v(r’ u(r) J%,,-(r') 


+2 { Lo(us(r) 2,0) 
x PW’, AW Dr", rr’) dr” (11.34) 


N,.(r, r’) =2 J [o(nuy(r’)—o(r" u(r", r')dr” (11.35) 
Coy = Syoy $2). | [ater ete | [He ror’) ir folr')dr' (11.36) 
v0 6 


u,(r) = Jz J14.3(k,, 7) 


mk, r 
(7) =i JE HIP ky 0) 


The integral equation for the partial amplitude is of the form 


(11.37) 


div = 2D | OO (r+ [WD i fdr yde' (11.38) 
vO 0 


We can carry out a transformation on (11.33), analogous to (2.39). In 
particular, we can introduce new functions x,, which are related to f, by 


A) — » Cyoy Nyy (1) (1 1.39) 
We then have 
je Ode ts inven Voar ee ine Xyye dr’ (11.40) 


v’ 0 0 
From equation (11.36) are obtained the following equations for c¢,,,: 


Cyov = 0p. 2 os Cvov! { Doveas ] vr") Wr", r’) in| Xyy(r’) dr’ 
? (11.41) 


On substituting (11.39) into (11.38), we obtain an expression for q,,, in terms 
of the functions x,,-: 


Qvov =e: 2 y C voy" | eral +| UO Wy, r’) ir Xyryelr’) dr’ 
= 0 e (11.42) 
Let 
Lyry = Sy 2D | (rH wt") + i Wd, FOC dar'|s tye) dr’ 
2 (11.43) 


The matrix L,-, 1s a generalization of the Jost function given by (2.41). 
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We can use the properties of Bessel functions to write the matrix element 
L,-, in the form 


bo hak, Pik? Bake (11.44) 
where 
Ayn, = Oye, + 20 LPG kT! 4 | (oe | Wind Xyye dr 
: : (11.45) 
Byoy = 2k Pk! 2» { bi. +f Wo yy ty, ir’) Meade (11.46) 
“0 


The quantities A,-, and B,-, are functions of k2 and k2.. Use can be made of 
(9.32), which represents the conservation of energy, to express all the k? 
in terms of kg. The diagonal elements are of a form analogous to (2.45): 


Ly, = Ay +ik2'2*"B (11.47) 
The system given by equation (11.41) may be written in the form 
> Oyo" | = Oy (I 1.48) 
and the integral expression (11.42) in the form 
Ivov = a Cy oye 2k B (11.49) 


In (11.48) the summation is carried out over all the states v”, whereas 
in (11.49) it is carried out only over those states which are consistent with 
the conservation of energy. If we neglect virtual transitions and retain in 
(11.48) only those states over which the summation is carried out in (11.49), 
the number of constants in (11.48) and (11.49) will be the same. 

If we consider all the possible values of vo, the quantities c,,, will form a 
matrix, which we denote by c. Equation (11.48) can then be written in 
the matrix form 

cL= 1! (11.50) 


As with c,,,, the amplitudes q,,, depend on the subscript vy, and form the 
matrix q,,,. Let us introduce the matrix k whose elements are k,d,,. With 
its aid the expression (11.49) may be written in the form 


q = —ck'Bk! (11.51) 
Substituting c = L~', we have 
q =—L™'k’Bk' (11.52) 
This is a generalization of (2.47) which expresses the scattering amplitude 
in terms of the Jost function. 
From the unitary relationships (10.60) we can obtain the corresponding 
relationship for the partial amplitudes. When the atom is in the ground 
state prior to collision we have 


Im qoo- kolaool7 = 


(11.53) 
and /? = 0,/2 = L. 
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If the energy of the incident electron is too low for excitation to occur, 
the sum on the right-hand side will vanish and we return to (2.29). Let 


Es Bi- I 
= 11.54 
doo ik, ( ) 
The partial cross-section for elastic scattering is then 
7 
Foo = (2L4 1) 75 |B 1? (11.55) 
¢) 
and the partial cross-section for all the inelastic processes 1s 
Tt 
Y, 90, = 2QL4+1) 5 (1-|B,]’) (11.56) 
v#0 ko 


which follows from (11.12). Hence the partial cross-section for all the 


inelastic processes is at all times smaller than 2(2L+4+ 1)/k6. 


12. ELASTIC AND INELASTIC SCATTERING CROSS-SECTIONS. 
BEHAVIOUR NEAR THE THRESHOLD 


Let us now return to the general properties of the amplitudes. To begin 
with, we shall transform (11.52) to a more convenient form, using (11.44). 
In matrix notation, the latter expression is of the form 


L = k'Ak7!+ik'Bk'*! (12.1) 
Substituting 
A =—BM (12.2) 
into (12.1), we have 
L = —k’B(Mk~!—ik!*') (12.3) 
Using (11.52) together with the latter expression, we obtain 
q = [k’B(Mk !—ik'**)]~*k'Bk! (12.4) 


or 
q = (Mk! —ik'*!)7 BOK 7 'k'BE! 
and hence finally 
q = k'(M—ik7!*")7'k! (12.5) 
Bearing in mind (10.57), we obtain the following expression for the 
T-matrix 
(22)*T = k'*3(M —ik7?*!)7'k +4 (12.6) 
This representation of the T-matrix is used in the general theory of multi- 
channel processes [29]. It is evident from the preceding discussion that this 
representation 1s a natural consequence of the properties of the Jost matrix. 


The formula (12.2) establishes the relation between the M-matrix and the 
Jost matrix. The matrix M is real and symmetric. which ensures that 


T*=T =n tT * (12.7) 
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We recall one further representation [30] also used in the theory of multi- 
channel processes involving the so-called reaction matrix K, which is given 
in terms of M by the expression 


M = k'*2K 7 Ii *4 (12.8) 
The amplitude q is 
q=k~?#K(1—iK)7'k7? (12.9) 
and the T-matrix 
(2n)*T = K(1—iK)7! (12.10) 
It follows from (12.2) and (12.5) that the reaction matrix 
K = —k'’*32BA7~!k'*2 (12.11) 


The expression (12.5) 1s convenient for the investigation of the dependence 
of the cross-section on energy, particularly near the excitation threshold. 

As an example, consider a simple two-channel system [29] and suppose 
that the incident particle is in the s-state both before and after the collision. 
In spite of the fact that this is a much simplified example, it does exhibit 
certain general features characteristic of real multi-channel systems. 

By (12.5) the elements of the matrix q are of the form 


My .—1k, 


C= D (12.12) 
M,> 
= eee 12.13 
Gi2 = 421 dD ( 3) 
M ,,—ik, 
ee ee 12.14 
422 D ( ) 
where 
D =(M,,—ik,)(M..—ik,)— Mj, (12.15) 
k, = Vk2—-2AE (12.16) 
and AE is the excitation threshold. Below the threshold let 
ky = ix = if2AE—k? (12.17) 
The effective cross-sections are then given by 
Key 42 
on = 40 ail? ae ee eee (12.18) 


When the energy of the incident particle lies above the threshold, 1.e. 
when £, = 4ki > AE, the scattering and excitation cross-sections are 
respectively 


M2,+k3 
O1; = 4 ———_—______**§ = ____ __________ (12.19) 
(M4,M.—My.—ky ky)° + (Ky Ma. +hk,M,,) 
: M2 | 
O.=4 K2 2 (12.20) 


= 
ky (MiyMa.—My.—ky ky)? +(ky Maa tk, My,)° 
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When the energy of the incident particle lies below the threshold, 1.e. 
when E, < AE, we have 


(M,,+K) 
a,, =42n—-——_—___- _—_.._.-"T-—-—"————; (12.21) 
= [M.1(M..+K)— M7)? +ki(M,. +k)’ 
which can also be written in the form 
4 
= = sin? 5 (12.22) 
ky 
where 
k2(M c)? 
sin? 6 = MiMi) (12.23) 


[M.1(M..+k)—M{2]? +kj(M2.+K)? 
Consider the behaviour of the cross-sections near threshold. 


To within linear terms in k, in the amplitudes |g, ,| we have above the 
threshold 


sa 4n|q6? |? (1 — 2k kala’) (12.24) 
k|q49/’ (0) (0){2 
Cae ay ee k.|q4\9| ) (12.25) 
1 


where the superscript 0 indicates that the particular quantity is taken where 
k, = 0, 1.e. at threshold. 
Below threshold 


O11, = 4nlqs?|7(L-— 2k 1gh9 |? cot 6) (12.26) 
O11 
7 
Ew zs Evn - 
(a) (b) 


Fic. 7. Elastic scattering cross-sections near the excitation threshold: (a) cot 6 < 0, 
(b) cot d > O (Ej, Is the threshold energy) 


Comparison of equations (12.24) and (12.26) shows that do,,/dE, has 
a singularity at threshold.t At energies above threshold, o,, decreases with 
increasing £,, and the character of the variation depends on the sign of 
cot 6 below threshold. If cot 5 > 0, o,, decreases with decreasing E£,. 
But, if cot 6 < 0, o,, increases with decreasing E,. 

The variation of o,, with £, is illustrated in Fig. 7. Apart from the 
singularity at threshold, the cross-section ¢,, can also exhibit a resonance 


+ Elastic scattering at the threshold for inelastic scattering has been investigated by 
many authors (see, for example, {31, 32]). 
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maximum in the sub-threshold region. As can be seen from equations 
(12.22) and (12.23), the maximum will occur when sin? 6 = 1, i.e. when 


M,,(M,,+n)—M?, =0 (12.27) 


The width of the resonance peak depends on the ratio M,,/M,,. It can be 
shown [29] that the peak is sharp when |M, ,/M,,| > 1. 

Figure 8 shows a typical variation of o with £ for the two kinds of 
threshold singularity and in the presence of resonance. 


(a) (b) 


Fic. 8. Elastic scattering resonance near the excitation threshold: (a) cotdéd < 0, 
(b) cot d > 0 


Numerical calculations carried out for a potential well [29] show that 
for some values of the parameters one does in fact obtain resonance peaks 
near the threshold. 

Equations (12.24}(12.26) are valid in the immediate neighbourhood 
of the threshold. It is possible to obtain approximate formulae which hold 
in a broader region by retaining in the amplitudes g,, terms which are 
quadratic in k,. With this intention, the quantities A, and B, in the Jost 
matrix (12.1) are expanded in powers of k2. An approximate formula 
for o,, has been obtained [33] in this way. A basically similar approach 
is used in [29], where the matrix M,, is expanded in powers of k3. Calcu- 
lating o,, by means of equation (12.21), in which the matrix elements M;, 
are represented by the two-term expansion 


Milk) = Mi(0)+ Rig” (12.28) 


it is possible to carry out an approximate study of the behaviour of o near 
resonance. This method is used in [34] to investigate the resonance in the 
elastic cross-section of the hydrogen atom for electron scattering, assuming 
a single excitation channel (excitation of the 2s-state). 

If the angular momentum of the particle after collision has a non-zero 
value and is equal to, say, /, then according to (12.5) the excitation amplitude 
91> near the threshold varies as kj. The excitation cross-section near the 
threshold is 


O12 ~ k3'*! (12.29) 
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The elastic scattering cross-section at the threshold will then no longer 
have a singularity of the form shown in Fig. 7a or 75 since do,,/dE 1s 
continuous at the threshold. 

The resonance in ¢,, below the threshold can also exist when / # 0. 
In the presence of three or more connected channels, the resonances can 
appear not only in the elastic scattering cross-section but also in the inelastic 
cross-section between two neighbouring thresholds. 

The case when there is degeneracy, 1.e. two or more channels at the 
same energy, occupies a special position. It is encountered in collisions 
with hydrogen atoms. Since the matrix elements V’,,- include elements 
which fall off as r~* as r > 00, the behaviour of o near the threshold will 
be different from that indicated by equation (12.29) [35]. 


Chapter 4 


THE ELASTIC SCATTERING OF AN ELECTRON BY A 
HYDROGEN ATOM 


13. THE SCATTERING OF AN ELECTRON BY A HYDROGEN ATOM 
REGARDED AS A CENTRE OF FORCE 


As has already been pointed out, the elastic scattering of an electron by an 
atom can be described in various approximate ways as the superposition 
of three effects: the central-force interaction between the electron and the 
nucleus, the exchange interaction, and polarization of the atom by the 
electric field of the incident electron. All these effects are allowed for formally 
in the exact equations for the function Fy. If we take for F, the system of 
equations (10.55), the terms Vo9 Fy and { Woo9Fo dr represent the central- 
force and exchange interactions between the electron and the atom in the 
ground state, while the sum 


3 (Vou F,+ i W,, F, dr) 


represents the effect of polarization and the connection between the elastic 
and inelastic processes. By neglecting one of the factors it is possible to 
obtain simple equations which admit of a practical solution. Furthermore, 
the relative importance of the various factors may be estimated by comparing 
the results of different approximate calculations. 

The crudest approximation involves the rejection in the equation for Fo 
of all the terms containing F, and the term describing the exchange inter- 
action. In that case the system of equations reduces to a single equation. 
For simplicity, we shall omit the subscripts so that the equation for Fo 
becomes 


AF +(k*—2V)F =0 (13.1) 
According to equation (10.20), 
2 
V(r) =— + fe 4 


Substituting the expression for i wave ae of the hydrogen atom in 
the ground state 


I 
Wo — 2 e' 
Vn 
and completing the integration, we have 


V(r) =- me (1 + ;) (13.2) 
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The asymptotic expression for F, 1s 
eikr 
Pagers q— (13.3) 
r 


This approximation is equivalent to assuming that the atom behaves as a 
centre of force. 


On separating the angular variables 


is. G = “ = 2V | f,=0 (13.4) 


for which the systematic solution 1s 


id} ] 
f, ~ Fain (kr =F re i) 


2 


Equation (13.4) has been solved numerically for /=0 and /=1 for 
0O<k* < 1-75 [36,37]. According to [37], the scattering length is 
a =-—9-44, The large value of the scattering length and the negative sign 
indicate the presence of a virtual level near zero energy. 

A reasonable approximation to the numerical calculations [36] can be 
obtained by the method of trial potential (§ 6). We shall take the trial 
expression for Vo in the formt 


1 1 
aa r<R 
Vo = R r (13.5) 
0 r>R 
where R is determined from the condition given by equation (6.7). 
For the s-wave at zero energy, we have 
1 J 
3-2(5- 7) fo=0 r<R 
: (13.6) 


= r2R 
When r < R, the solution given by (13.6), which satisfies the condition 


fo(0) = 0, is 
fo= cre F(t ae 2, 24/27) (13.7) 


where F is the degenerate hypergeometric function and C 1s a normalizing 
constant. When r > R the function fo is of the form 


fo=rtq (13.8) 


t Expressions of this form were used in the 1930’s to investigate certain general 
properties of the phase 6 [38]. In this paper the parameter R was determined with the 
aid of the semi-empirical rules of Slater for the atomic radii (developed in the theory 
of atomic speetra). However, the phase o is very sensitive to changes in R, particularly 
at low energies. Therefore intuitive considerations leading merely to an order of magnitude 
for R are inadequate in the present context. 
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or on introducing the scattering length ag = —g, 

fo =1r—-ao (13.9) 
The continuity condition for the function and its derivative at r = R yields 
relations giving C and a, in terms of f/,(R) and fo(R). In particular, 
fo(R) 
F(R) 
The parameter R can be determined from equation (6.7). In the present 
case, this relation assumes the form 


Aag=R- 


R oO 
{(v- Vo) f2 dr +{ Vf2dr=0 (13.10) 


This equation has been solved numerically, and the result obtained for R 1s 
R= 1-43 (13.11) 


For this value of R, the magnitude of 2 [ |5V|r dr turns out to be equal to 
0:33. Therefore the condition given by (6.13) is satisfied and one can expect 
good agreement between wave functions calculated with V and Vo. 

The results of the calculations are as follows. The scattering length is 
a =— 10-3, quite close to the exact value of a = —9-44 [37]. 


14 f 
I2 fo 
lO 
8 
6 
05 1 2 3 4 5 


P 
Fic. 9. The form of the functions fo and ffor k = 0 


Figure 9 shows a plot of fg and a plot of the values of f found in [37] by 
numerical integration of the equation 
f"-2Vf=0 
In order to estimate the degree of approximation reached with the expres- 
sion for Vg with k 4 0, consider the s-wave for k # 0. The corresponding 
equation 1s of the form 
fo t(k? —2Vo) fo = 0 (13.12) 


6—2 
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The solution forr < R is 


1 
a= cre F(I —-, 2, 240 (13.13) 
pi 


2 
w= f2—e (13.14) 


For r > R we can write the solution in the form 


9 = Sin(kr+69) (13.15) 
Both the constant c and 6, can be determined by matching equations (13.13) 
and (13.15) atr = R. We note that if is an integer, the hypergeometric 
function becomes identical with the corresponding Laguerre polynomial 
and the calculation becomes particularly simple. For example, when 


2 
[a O40 (13.16) 


where 


and ju = 1, we have 


Io=cre” (13.17) 
A function f suitable for the comparison (A* = 0°45) is given in [36]. 


Figure 10 shows the values of fg and f. The phases 0, and o are, respectively, 
equal to 1.11 and 1.01. 


Fic. 10. The form of the functions fo and f for A? = 0:40 and k? = 0°45 respectively 


A considerable number of calculations have been carried out by the 
variational method. 


As in the analysis considered in § 7, one constructs the functional 
L.= | Wri re) HOW ra) dry dry (13.18) 


where // 1s the encrgy operator for the entire electron + atom system. The 
trial function yw, should contain the atomic wave function which is not varied. 
The expression 
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is stationary with respect to variations in the trial function, and when the 
exact wave function w is used the parameter Q is equal to the scattering 
amplitude. If the trial function contains nm parameters c, then by Kohn’s 
method these parameters are defined from the stationarity condition 


OL, OL 


= 0, 2 13.20 

de ag anne 

while the scattering amplitude is given by an expression analogous to (7.8): 
1 

te Ge (13.21) 
Tt 


If the atom is taken as equivalent to a centre of force, the trial function 
is of the form 


Willi, 72) = Vols) F (72) (13.22) 
where Wo, 1s the atomic function and F, 1s a function describing the scattered 
electron and containing the varied parameters. 

In order to obtain the partial scattering amplitude, qg,, the function F, 
in (13.22) must be replaced by 
f(r) 
Fa = Yio(9.9) (13.23) 
For example, consider the variational calculation of the scattering 


length in [37]. Various forms of trial function were studied including the 
two-parameter function 


f,=r—a(l—e *")+ bre" (13.24) 
The Kohn scattering length was found to be 
a = —9-38 (13.25) 


14. EXCHANGE EFFECTS 
A better approximation is obtained by retaining the exchange operator 
on the left-hand side of equation (10.55). This leads to the following 
integro-differential equation for Fg: 


AFE+(k2—-2V) Fz = 2 { Wer’ r)FE(r’) dr’ (14.1) 
where V is given by (13.2) and 
: J ae 
Wn) = (a5 — 5) vol ol? (14.2) 
Jr’—rl 2 2 


If we separate the angular variables using the notation of (11.5) and omit 
both the superscript + and subscript 0, then 


d*f, 


: 
dr- 


(+1 r 
+ G — te — 2M f= 2 | Wn nA’) dr’ (14.3) 
0 
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where 
eth k* 
Wir’ r)=+ 4rr’ e~"~? | (7, r’) — oo (14.4) 
subject to the conditions 
J,(0) = 0 
- - (14.5) 
f~ eee sin (19 ay +- 3) 
The scattering length in the triplet and singlet states is [35] 
a, = 8-058, © a_ =2:347 (14.6) 


The wave functions are substantially modified by the introduction of the 
exchange interaction. For example, Fig. 11 shows the functions f*(r) and 


Fic. 11. The form of the functions f + and ffor k = 0 


f(r) for k = 0 in the absence of exchange. The phases 6g, df and 5} were 
calculated in [39] forO0 <A <1. The results of the ealeulation are illus- 
trated in Fig. 12. Plots of 69, 6, and 6, which do not include exchange 
effects are given for comparison. It is evident that in many cases the 
phase 6 lies between 6* and 67 in the absence of exchange. 

It is shown in [17] (§ 15) that in sufficiently accurate variational calcu- 
lations the phase 6 in the absence of exchange should tend to the half-sum 
4(6* + 6°). 

The numerical solution of (14.3) must be carried out by the method of 
successive approximations owing to the presenee of the infinite integral 
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3 
(a) 
2 
50 
' 
bo 
56 
00501 O25 05 Re 
42 
0-4 5; 
(b) 
OO? 
6) : 
6, 
-O0:-2 


O25 O05 ne 


Fic. 12. The phase shifts 6+ and dp (a), 6; + and 6, (6), and 62+ and 6, (c) 


on the right-hand side. It is shown in [39] that nine iterations are necessary. 
In this connection it is convenient to transform equation (14.3) to a form 
permitting direct solution without iteration. This transformation was 
carried out in [40] for the integral equation (11.35) which is equivalent to 
(14.3), and later by basically the same method for equation (14.3) [41]. 

We now consider a version of this involving the integral equation for 
the partial wave / = 0. Equation (11.33) assumes the form 


sin kr 
k 


f=C ee | M(r,s)f(s) ds + | N(r, s)f(s) ds (14.7) 
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where 


2 2 
M(r,s) = ; sink(r—s)¥(s) + ral sink(r—N[WO(s, DN—-WO(s, D] dt (14.8) 


AY 
r 


2 
N(r,s) = al sink(r—1)wes, hat (14.9) 
O 
a bea) | oe V(r) + } W (r,s) es is faryar (14.10) 
0 0 
To begin with, we consider the singlet case. Let f = cx, so that 
° k : 
x= = : + | Mer s)x(s) ds +| N(r, s)x(s) ds (14.11) 
0) O 
! 
c= OS (14.12) 


1—2 } Ly er + } W (r,s) el is x(r) dr 
0 0 


Inspection of the expression for M(r,s) given by equation (14.9) shows 
that it is the product of the function of r: 


r 


8 
R(r) = = | sin k(r—tye 'tdt (14.13) 
0 
and the function of s: 
agp eo. we 
O(s)=se *(-—k*-—1 (14.14) 
S 
Le: 
N = R(r)Q(s) (14.15) 
Bearing this in mind we shall rewrite equation (14.11) in the form 
sin kr 
x= P + | Mads + R(re (14.16) 
O 
Co = } xQ ds (14.17) 
0 
Let 
XN =NX,+CoQXo (14.18) 


If we substitute (14.18) into (14.16) and require that x, satisfy the equation 


r 


1, 
x= ea Kr + | Mr s)x,(s) ds (14.19) 


0 
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we have the following expression for xo: 
2 
Xo = R(r) + -| M(r, s)x0(s) ds (14.20) 
0 


Equations (14.19) and (14.20) can be solved numerically without iteration. 
Therefore the problem formulated above is solved with the aid of trans- 
formation (14.18), which is possible because of the property indicated by 
(14.15). The unknown constant cg can be found in terms of x9 and x,. 
Substituting equation (14.18) into (14.17), we obtain 


[ Ox, ds 
Cy = — > (14.21) 
1—| Oxods 

The triplet case has a specific property associated with the fact that the 
antisymmetric function 


Wr 1s T2) =f Wolr2) -—f(r2) Wo" 1) (14.22) 
is invariant with respect to the transformation 
f(r,) > f(r,) + const Wo(r) (14.23) 


It may be easily verified that equation (14.16) is also invariant with respect 
to (14.23) in the triplet case. This allows one of the two constants c or Co 
to be chosen arbitrarily. If, for example, we set cg = 0 then 


sinkr 


k 


+] M(r, s)x(s) ds (14.24) 


We shall now show how to transform equation (14.3) to a form not requiring 
iteration. 


We have 
f"t+ G — “ — a f= 2[ wr, s)fds (14.25) 
The right-hand side of this equation may be die 
{wha = [ww Fas + [ werpas (14.26) 
Since 247°? sites saitea the form 
2) = A(r)B(r’) (14.27) 


we have 


2 WOFds = A(r)c (14.28) 
O 
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where 
c =| Bfas (14.29) 
Therefore, : 
f! + [ee a Ds ay f =cA(r) (14.30) 
Substituting 
S=x-+cy (14.31) 


we have two equations for x and y which can be solved without iteration: 
ldi+1 
ae G au ( 5 ) 
r 


— - x . ds (14.32) 


_! l ze 


The unknown constant c can be found as bates by substituting (14.31) 
into (14.29). The result is 


C =| Bs is[( ~ [By ts (14.34) 


In variational calculations, exchange effects are taken into account 
by using a symmetrized trial function y, in the form 


W, = Fry )Wolr2) EF ilra)Wolrs) (14.35) 


For example, we mention a variational calculation of the scattering length 
carried out in [37], where the trial function f(r) = rF(r) was chosen in the 
form given by (13.24). The scattering lengths were found to be 


a, = 8-11, a_ = 2°38 (14.36) 


15. THE EFFECT OF POLARIZATION OF THE ATOM ON THE 
SCATTERING OF AN ELECTRON 


The effect of the polarization of an atom on the scattering of an electron 
may be investigated by including the term 


3 (VorF. +] WoaF dr) 


in the equation for Fo [see (10.55)]. In spite of the fact that this term is the 
sum over all the excited states a, the polarization effect is not additive since 
in addition to the terms 


VioFo + | WaoFodr 


in the equation for F, there are also the terms 


) (4 Pik | WF, dr) 


a’ #2 
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which represent the coupling between the different states. The effect is 
additive only when this coupling is neglected and the coupling 0 > @ is 
taken into account. 

We shall show that in the limit of large r the effect of excited states is 
equivalent to an additional interaction between the electron and the atom 
which can be characterized by a potential energy of the form U = —f/r* 
[27], where f is the polarizability of the atom. We shall confine our attention 
to the approximation in which only the 0 — «& coupling is taken into account. 
Under these conditions the system (10.55) becomes 


AF) +(k2—2V0)Fo—2 { Woo Fo dr’ =25. (VoaF a ie { WF, dr’) (15.1) 


a#QO 
AF, +(k2—2V,,)F,—2 { W..F,dr' = 2(VaoFo fe { W.oFo dr’) (15.2) 


If the energy of the incident electron 1s such that excitation can take place, 
the right-hand side of (15.1) contains both the effect of polarization and the 
effect of coupling between excitation and scattering due to the unitary 
relation. But if the energy of the incident electron is insufficient for excita- 
tion, then the right-hand side of (15.1) contains only the polarization term. 
We shall begin with the latter case, i.e. we shall suppose that kz < 0 for 
all « £ 0 and only k§ > 0. Consider equation (15.2) for large r. If k, < 0, 
the solution of the homogeneous equation AF,+k2F = 0 which is finite 
for r > oo will have the asymptotic form exp (—|k,|r). Under these con- 
ditions, the asymptotic behaviour of (15.2) is determined by the right-hand 
term. Since 


k2 = k2+.2(E,—E,) (15.3) 
we have 
AF ,+ (k3—2V,.)F,—2 | Walt", r)F(r’) dr’ + 2(Eo —E,)F, 
29) 5F 29 { W.o(r', KF o(r’)dr’ (15.4) 


For large r the expansion of V,, in terms of the reciprocal powers of r 
begins with r~'~!, where / is the angular momentum characterizing the state « 
(when / = 0, V,, falls off exponentially with increasing r). The kernel of 
the exchange operator W,, falls off exponentially as r + oo. Therefore for 
large r the right-hand side of equation (15.4) has the asymptotic form 


oikor 
yin d io ne q 
r 


The leading term on the left-hand side of equation (15.4) for large r is then 
(E,—E,)F,, so that we may assume that the asymptotic expression for F, is 


Viol 
F,~- eae 
bg ig 


(15.5) 
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It is easy to show that for this asymptotic behaviour of F, the expression 
(A+k2—2V,,)F,—2 | WF dr’ 


in equation (15.4) may be neglected when r —> od. In point of fact, according 
to equation (15.4), the term in F, which exhibits the slowest variation is 


r~!~! exp (ikgv.r), whilst in 


(A+k2—2V,,)F,-2 | W,,Fdr’ 
the corresponding term is r7'~ 
(15.1) and let 


2 exp (ikpv.r). We substitute (15.5) into 


Vi o(r)|? 
U(r)=—-) = (15.6) 
Wor nV io(r 
Wee) = — el (15.7) 


We then have 
AF 9 +(k3~2Vog—2U)Fp = 2 | (Woot W,)Fo dr" (15.8) 


In this equation U enters as an addition to the potential energy due to 
excited states, while W, is an addition to the exchange operator. V9. JV 
and W, fall off exponentially as r > oo. As regards U, we need only con- 
sider the slowly varying part, i.e. the dipole terms (/ = 1). For large r, 
we then have 

2 S 


2~ ale 
V0 


(the factor 2/3 appears as a result of the integration with respect to angles 
within the matrix element), so that 


ya | 


(15.9) 


a) poe 
The quantity B=s y = is the polarizability of the hydrogen 
atO0*aq = "0 


atom. Therefore, for large r, the function U(r) tends to —f/r*, which is 
the interaction energy between the electron and the induced dipole. Cal- 
culations [27] show that about 66% of f is due to the 2p-state. The discrete 
states contribute 81-4% in all, the remainder being due to the continuous 
spectrum. The numerical value of f is 9/2. 

Let us now return to the case where excitation is possible, 1.e. one or 
more of the A? are positive. We can substitute 


: ad V0 Fo 
| oe (15.10) 
E. ~E, 
in equation (15.2) for the corresponding «. The leading term in the 


ik, 


asymptotic form of F, will be the diverging wave g, ae On substituting 
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(15.10) into (15.1), we have the following equation for Fj(r) 
AP Ake OVOP 20 | (Woo + W,)Fo ar’ 
=25 (VoF, ‘ { WF, dr’ (15.11) 


a#t¥O0 
When r —> o, the right-hand side of this equation expresses the coupling 
between elastic and inelastic scattering channels. Similarly to equation 
(15.8), the left-hand side includes the polarization interaction U. 

The presence of the slowly decreasing potential (~ r~*) in (15.8) or 
(15.11) has an important effect on the asymptotic behaviour of the wave 
function. This can most easily be established in the limiting case of zero 
energy. For large r, the exchange terms fall off exponentially and are of 
little significance as far as the asymptotic expressions are concerned. We 
shall therefore consider the equation 


AF-—V(r)F =0 


The function F has the asymptotic form 


const 
Fr~ti+ 
For large r we shall set 
AF = V(r) 
and hence 
haf i V(r’) os 
An lr—r | 


Completing the integration with respect to the angles we obtain 


1 
F=1- “| Vr’? dr’ -| Vr' dr’ 
r 
0 r 
or 
1 2 / I 2 
F=1-—-—|Vr'*dr' —| Vr’ dr’ +—- 4 Vr’ dr’ 
r r 
0 r r 
which amounts to the same thing. The function 
' ' I 12 , 
J=|Vr'dr' —-| Vr’*dr (15.12) 
r 


determines the nature of the approach to the asymptotic form. 
If V falls off as exp (—/r) when r > o0, J will also fall off exponentially: 


-; l 
J=e “(1 +=] 
Ar 
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If, on the other hand, the polarization potential V ~ r~* is present, then 
1 
= 3p 
In the absence of the polarization potential we have for large r 


Fw 1 +-——-+conste7* (15.13) 


const 
r 


i.e. the approach to the asymptotic form is exponential. In the presence of 
the polarization potential, 


const const 
Pettey (15.14) 


r 


and the asymptotic function is approached slowly in accordance with 
the r~* law. 

It is possible to introduce polarization formally at finite distances, 
not only for r— co. However, it will then no longer be a function of 7 
as given by equation (15.6) but an integral operator [42, 43]. We shall 
show this for the case of the simplified system of equations in which the 
exchange terms have been neglected. Thus we have 


AF, +(k2—2V,,)F, = 2V 9 Fo (15.15) 
AF o+(kg—2Vo0)F 0 = 2 me Vout (15.16) 


We shall confine our attention to the case kz < 0. Suppose that G, is the 
resolvent defined by 


G, = (A+k?+i0—2V,,)7! (15.17) 
The solution of (15.15) is then 
F, = 2G6,V oF (15.18) 
Substituting this into (15.16) and introducing the operator 
U = 2 Vor. Vio (15.19) 
we have 
AF) +(k2—2Vo9—-20) Fy = 0 (15.20) 


Let us now consider the numerical calculations which include polariza- 
tion effects. Calculations of the phase shifts for L = 0, 1, 2 with allowance 
for the 1s—2s coupling in the interval 0 < k* < 0-75 are given in [44] 
(k? = 0-75 corresponds to the threshold for the excitation of the 2s-level). 
Calculations in which the ls—2s—2p coupling was taken into account 
are reported in [45]. When the 1s—2s coupling is allowed for, it is found 
that 

ay = 7:5, QaL= 2°33 
while the result for 1s—2s—2p coupling is 


a, = 6-74, a_ = 1-89 
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The phase shifts 6g and 69 are plotted in Fig. 13. Calculations without 
polarization and with polarization allowed for by the variational method 
(see below) are given for comparison. 


O02 O04 O06 O08 Oo O02 04 O06 O8 


k2 k* 


(a) (b) 


Fic. 13. The phase shifts 69+ (a) and d)~ (6): 1—without allowance for polarization, 
2—with ls — 2s coupling, 3—with ls — 2s — 2p coupling, 4—variational method with 
polarization 


The relative importance of the various couplings may be deduced from 
Table 3 which gives 6g for two values of k?. As can be seen, Is—2s and 
ls—2p coupling give rise to roughly the same contributions. The simul- 
taneous inclusion of 2s and 2p states leads to a much greater contribution 
than the sum of the ls—2s and ls—2p contributions. Therefore, the 2s—2p 
coupling is significant and the effect of polarization is not additive. 


TABLE 3 


* ra 
States included 85 for k? equal to 


0:55 0-60 
Is 0-7004 0-6704 
Is — 25 0:7352 0:7115 
ls — 2p 0-7338 0-7040 
Is — 25 — 2p 0-7846 0-7707 
ls — 2s — 2p — 3s 0:7894 0:7738 
ls — 2s — 2p — 3p 0:7920 0:7738 
ls — 2s — 2p — 3d 0:7841 0:7694 
ls — 2s — 2p —3s — 3p 0-7975 0:7814 
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Another feature evident from Table 3 1s that the coupling with the 3s, 
3p and 3d states adds very little to the contribution due to the ls—2s—2p 
coupling. On the other hand, it cannot be considered that the inclusion of 
only the 2s and 2p states is sufficient since the contribution to static polariza- 
bility of the hydrogen atom due to the 2p state is 667%. It should therefore 
be assumed that the convergence is very slow and a large number of excited 
states must be taken into account to improve the approximation which is 
achieved with the Is—2s—2p coupling only. 


O O02 04 06 08 !10 12 14 
ke 


Fic. 14. Peak in the elastic scattering cross-section [44] 


As can be seen from Fig. 13, the phase shift increases very rapidly near 
the threshold for the excitation of the 2s-level. The variation in the phase 
shift in the region of rapid increase can be described by the Breit-Wigner 
formula with resonance energy E = 9-61 eV (A? = 0-711) and a full width 
at half height [C = 0:109 eV [45]. As can be seen, this is a very narrow 
resonance whose presence is predicted by the general theory (§ 12). In 
addition, there is also a threshold peak. Figure 14 shows the scattering 
Cross-section as a function of energy with a clearly defined pcak (44). 
Analogous results, i.e. the narrow resonancc, were reportcd in [34]. 
It should be noted that detailed analysis of the variation in the 
phase near the 2s excitation threshold was madc difficult in [44, 45] 
because the convergence of the iteration method used to solve the 
integro-differential equations deteriorated as threshold was approached. 
It is evident from Fig. 13 that the phase shift 69 is not modified to any great 
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extent when the ls—2s—2p coupling is taken into account. Figure 15 
shows k cot 6¢ as a function of k*. As can be seen, k cot dg is appreciably 
non-linear for small £7. This is due to the terms & and k? In k in equation 
(3.24). In the case of k cot 5g the departure from linearity is very slight. 


+ 
16) 


k cotd 


0'05 O10 0:15 
‘2 
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Fic. 15. Dependence of k cot d9*+ (a) and k cot dg~ (6) on k?: 1—with ls — 2s — 2p 
coupling, 2—without polarization 


This is due to the fact that the coefficient of A 1s inversely proportional to 
the square of the scattering length and is therefore much smaller in the 
singlet state than in the triplet state. 

Figure 16 shows the phase shifts df and 5;.+ Here again one can see 
a rapid increase in the phase near the threshold for the excitation of the 


(b) 
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Fic. 16. The phase shifts 6,* (a), d:~ (6), d2* (c) and 627 (d): 1—without polarization, 
2—with ls — 2s coupling, 3—with 1s — 2s — 2p coupling, 4—by the method of polarized 
orbitals 


+ For comparison, the results of calculations including the polarization effects and 
using a different method are also given (see below). 
E.A.C. 7 
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2s-level. As has been pointed out in [45], the phase shifts 67 and 67 approach 
each other at low energies and are adequately represented by equation (3.25) 
if B is replaced by 66% of the true polarizability of the hydrogen atom, 
i.e. by the contribution due to the 2p state. Thus, when k* = 0-1 the values 
of 5 calculated in [45] are 67 = 0-0084 and 5; = 0-0096, whereas equa- 
tion (3.25) yields 6, = 0:0089. 

There is another way of including the polarization effects. In this method 
the basic functions for the expansion of the wave function of the atom + 
electron system are taken to be the wave functions of the atomic electron in 
the field of two centres of force, i.e. the nucleus and the incident electron, 
rather than the atomic wave function as in (10.17) or (10.46) [46]. This 
method is a development of the so-called method of perturbed stationary 
states described in the book of Mott and Massey. In [46] it is referred to as 
the method of polarized orbitals. 

We first consider the elastic scattering in the absence of exchange. We 
shall write the wave function for the atom + electron system in the form 


W(r1, 2) = F(ry)P(ry, 12) (15.21) 


where F(r,) describes the incident electron and ®(r,, r,) describes the atomic 
electron in the field of the nucleus and the incident electron. The function 
@(r,, r,) can in its turn be written as the sum 


(11,12) = Wolr2) +9015 12) (15.22) 
where y,(r,) 1s the wave function of the atomic electron in the isolated atom 
and g(r,, 1.) 1s the change in the atomic function at r, duc to the incident 
electron at r,. The equation for F can be obtained from cquation (10.24): 


{ vo(r2)(H-EW(ry, ra)dry = 0 (15.23) 
where & = 4k?7+£, and the cnergy of the ground statc of the atom is 
Eo = —}i, 

We substitute the expression for Hf given by (9.1) and the expression for 


W(r,, 72) given by (15.22) into (15.23). We shall, moreover, require that 
g(r,,"2) and w,(r,) should be orthogonal for allr,: 


Ge r,)Wo(r2) dr, = 0 
We then have 
AF +(k?-2V—2V,)F =0 (15.24) 


where V is the potential cnergy of the clectron in the field of the undisturbed 
atom given by (13.2) and V, is the polarization correction to the potential 
energy and is given by 


Wolr 2g » 12) 


Ir; 1] 


V,(r;) = dr, (15.25) 


The problem is therefore reduced to the determination of the function 
ry, 2). 
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The equation for g(r,, r.) can be obtained by substituting (15.21) and 
(15.22) into (H—&)W(r,, r.) = 0 and using (15.24) and the equation satisfied 
by wWo(r,). In the present method, the determination of g(r,, r.) is carried 
out on the adiabatic approximation in which the position of the incident 
electron r, 1s regarded as fixed and all the derivatives of g(r,, r.) with respect 
to r,; are equated to zero. Moreover, if itis assumed that g is a small correction 
to Wo, the determination of g(r,,r,) can be simplified still further and g 
can be neglected in all the terms into which the sum W)+g enters. We then 
have the equation for g 


1 l jt 
(442 += — 4) o¢rs) = [= = ver) | oles) (15.26) 

r2 Ir;—r.| ry 
An accurate solution of this equation can only be found by numerical 
integration. However, the asymptotic behaviour of g for r, > r, can be 
found in analytical form. Whenr, — oo the potential V(r,) falls off expo- 
nentially and we can neglect V(r,) on the right-hand side of (15.26). Next, 
if we expand |r, —r,|~* in terms of the spherical harmonics and retain the 

dipole term only, we have 
1 1 

— x 20s 0 (15.27) 


lr, — 1 Py Oy 


where @ is the angle between r, and r,. When these simplifications are 
introduced, equation (15.26) is replaced by the asymptotic equation 


1 r 
(44: += = 4) golrista) = 273 volts) cosd (15.28) 
2 1 
where go is the asymptotic function. The solution of this equation 1s 
cos0d _. 
Jolr1, 12) = > 7=e *(ry +47) (15.29) 


This solution can be augmented by the solution of the homogeneous equation 
obtained from (15.28) by rejecting the right-hand side. However, since g 
and Wo are orthogonal, the additional function is equal to zero. 

As regards the approximate solution of (15.26), we must first mention 
[47] which gives a variational method of solving this equation. The method 
is based on the fact that the function 


t= { o(rsrs)|(d2 += = 1) ore) ~ ole) | dra (15.30) 
r3 nr 
is a minimum with respect to an arbitrary variation of g (for any r,), and 
its stationary value is equal to the polarization potential V, given by equation 
(15.25). Another approximation is described in [46], where the polarization 
potential is calculated by replacing g by its asymptotic form gg for r,; 2 rz 
7) 
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and assuming g = Oforr, < r,: 
g = &(r1,12)9o 
1 riers (15.31) 
E= 
0 ri <r, 
Having determined g(r,,r,) by any of these approximate methods we can 
calculate the polarization correction to the electron—atom interaction energy 
from equation (15.25). 


In order to take the exchange effect into account we must construct the 
symmetrized function 


Wry, 12) = (ry) (Wolr2)+9(1, ryt Fe aLWolri)+9(r2,11)] (15.32) 
where g Satisfies (15.26) as before. Substituting equation (15.32) into (15.23), 
we obtain an equation of the form 


AF +(k?—2V—2V,)F = 2 | (W+W,)F dr’ (15.33) 


where W, is the polarization correction to the exchange operator. As before, 
the angular variables of this equation can now be separated and equations 
with radial functions obtained. 
The equations derived from the approximate expression (15.31) for g 
are [46]: 
fe. 6) 


d (d+ 
Tat |e P -ar| i fore. n gar 
r 


= 22 ht NKFPA) (1534) 


where ¥ and /’,, are given by equations (13.2) and (14.4). 

The right-hand side contains terms describing the polarization effect. 
The first term contains 9/2r*, i.e. it is the asymptotic form of the polarization 
potential with which we are already familiar. The remaining terms are: 


N =— fe *(r? +3774 9r? +42r? +-22r+22) (15.35) 


P(f) = 4 +k)6,, foes 


l ee e "fir)dr 

GN ee ea ee 

" armas’ | pie? 

(1+1)r'*! e~"fir)dr 
(21 +-1)(21 4-3) pees 


r 


| e~'(r+4r?) 


oH) r (15.36) 


+ 4r°e7' di, { 


r 
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The equations (15.34) are transformed to a form not requiring iteration 
and solved numerically for / = 0, 1,2 and k between 0 and 1. The wave 
functions are calculated up to r = 20 after which they are assumed to be 
identical with the asymptotic form 


It 
sin (11 irs + y 


Figure 16 shows plots of 6¢ and 5. The scattering lengths are found to be 
a” = 5:8,a~ = 1-9, whereas the corresponding values without the polariza- 
tion corrections are respectively 8-1 and 2:35. However, subsequent calcu- 
lations [48] show that even for r > 20 there 1s an appreciable contribution 
due to the polarization potential. The change in the phase shift 6) between 
r= 20 and r = © can most easily be determined from equation (3.35) 
which for / = 01s of the form 
i, =— ae eT eee (15.37) 
dr k 

For r > 20 all the exchange terms in (15.34) are practically zero and 
this equation is valid. At such large distances we need only retain the 
asymptotic part £/r* in the potential V. If, in sin? (Ar+6,)), we replace 
Oo(r) by its limiting value 6 (co) which is valid for small changes in the 
phase shift, we have 


2 pein? kr+0 
0o(00) — do(R) == | in (15.38) 
R 
Hence, to within terms of the order of k, 
1 Ba? a? 1 Tp 
= — =(-— — f-, — k— 2Ink 
-[86(00)-34(R)) = (F 55 — Bos + Bg) +R + of? Ink) (15.39) 


where a is the exact scattering length. This formula was derived in another 
way in [48]. 

We note that equation (15.39) contains the linear term Axf/3 which is 
identical with the linear term in equation (3.24). When & = 0 we obtain the 
relation between the accurate scattering length a and the scattering length 
a(R) which is obtained from the asymptotic form for R = 20: 


lat al 
a=a(R)-P (555-73 +R) (15.40) 
Calculations yield @~ = 1:74. The correction to the singlet scattering 


length a* is less significant because a~ is small. 

In variational calculations, polarization is taken into account by including 
in the trial function terms containing r,, = |r,—r,|, so that W(r,, r,) is not 
reduced to the symmetrized combination of the products of functions of r, 
and functions of r, of the form given by equation (14.35). 

Determination of the scattering length is particularly important in 
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variational calculations. In the case of scattering of an electron by an atom, 
it turns out, just as in the simpler case of scattering by a central force (§ 7), 
that the scattering length calculated by Kohn’s method is the upper bound 
on the true scattering length [49]. This provides us with a criterion for 
estimating the accuracy of the variational calculations. In the triplet state 
we have a@~ < ag where ag is Kohn’s scattering length: 


eae ae) i Ww (H—E,W, dt (15.41) 


H is the energy operator for the electron + atom system, Eo 1s the ground 
state energy of the hydrogen atom, y, is the trial function whose asymptotic 
form 1s 


Wr = (1—P12)Wolrs) (“ = . (15.42) 


and Wo is the atomic wave function. 


In the singlet state, we have the following inequality which is analogous 
to (7.29): 


<a, 42 [vi (H—E,)W,; dt — woe ver E,W; a (15.43) 


where W’ and E’ are the eigenfunction and energy eigenvalue of the bound 
state of the electron + atom system (it is assumed that there 1s only one 
such state). The function yw’ must be a sufficiently good approximation to 
the eigenfunction for the negative hydrogen ion in order to ensure that 


E’= | WH’ dt < Eo (15.44) 
The inequality given by (15.43) can also be written in the form 
a* <a; (15.45) 


provided that one uses the function w varied by Kohn’s method and including 
the bound-state function w’, so that 


Y,=W,+ Ap’ (15.46) 
where A 1s the varied parameter. 

It is shown in [49] that even if the function w does not have the structure 
indicated by equation (15.46), the inequality given by (15.45) 1s still obeyed 
under certain conditions. 

We now briefly review the numerical calculations. The most systematic 
are given in [50]. The s-wave phase shift was calculated from the wave 
function W for the electron + atom system of the form 


w= pty 
where 


rie 


=> r.) (15.47) 
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and 
1 k _k 
jj ee ine a (= 5n) (15.48) 
k, k 
ag etic 2 
= Yo Cyomne 2? MoCTritrird (15.49) 
I,m,n>0 


. i , 
The function > tan 6 was defined as the stationary value of the functional 


k 
- tg6=4—[ W(H—B)p dr, dr, (15.50) 


The set of numbers /, m,n in equation (15.49) was chosen so that the 
sum /+ 2+ varied between 2 and 7 in the triplet state and between 1 and 6 
in the singlet state. Correspondingly, the coefficients C,,, , formed matrices 
of rank 3, 7, 13, 22, 34 and SO. 


l 
Figure 17 shows the stationary state values of — ; tan Oo for the triplet 


s-wave with k = 0-4 (about 2 eV) as a function of the scale parameter kx. 
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Fic. 17. Dependence of — (1/k) tan d9- on the scale parameter « for k = 0-4; the 
numbers indicated on the curves correspond to matrices of the following rank: 34 (1), 
15 (2), 7 (3) and 3 (4) 


When xk = 1:55, the curve corresponding to the matrix of rank 7 exhibits 
a discontinuity. Similar properties are exhibited by all matrices of higher 
rank but they are so narrow that they have not been resolved. 

When the number of parameters is large, the curves become smoother 
and approach each other. If it is assumed that as the number of parameters 
is increased all the curves tend to a single common curve, then it is possible 
to determine the corresponding phase shifts by extrapolation. Figure 12 
shows the extrapolated phase shifts for the triplet and singlet states. 
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For the calculation of the scattering length the function g was chosen 
in the form 


7 a 2:25 re teps 
p= (lt p)2e" [1-9 f,-=Ff- BE feos] (15.51 
ry ry “ry 
where 
1 
n-1 
— y,i/l P dp — pry, — 

fr=ry (7—))! : (n = 1,2, 3) 

0 


The dependence of g on r, was chosen so as to reflect correctly the effect 
of the polarization of the atom on the asymptotic form of the wave function. 
Figure 18 shows the scattering length as a function of the parameter x. 
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Fic. 18. Scattering lengths a, and a_ as functions of the parameter x. 


Table 4 gives a summary of the scattering lengths calculated by the various 
methods. 


TABLE 4 


ls — 2s —2p Method of polarized Variational 


coupling orbitals method 
a, 6:74 5°8 5:96] 
Q_ 1:89 1-74 1-768 


16. DISPERSION RELATIONS FOR THE SCATTERING AMPLITUDE IN 
THE CASE OF THE HYDROGEN ATOM 


The dispersion relations for the scattering amplitude in the case of the 
scattering of an electron by the hydrogen atom can be deduced by means of 
an analysis analogous to § 4 [5], 52]. As in § 4 the dispersion relations 
are based on the properties of the Green’s function G(E+£,+i0) for the 
atom + electron system, where £ is the energy of the incident electron, 
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and Eo is the energy of the hydrogen atom in the ground state. We shall 
assume, by analogy with the scattering of a particle by a centre of force, 
that G(E+£,+i0) is an analytic function in the plane of the complex 
variable E with a cut along the positive real axis, and has simple poles on 
the negative real axis at the points E = —|E,| which correspond to the bound 
states of the negative hydrogen ion H~. The residues at these points are 
equal to ui(ry, rz (ri, 12) where w,(r,,r2) is the eigenfunction of H™. 
Moreover, there are branch points on the real positive axis which correspond 
to thresholds for the excitation of the atomic levels. 

Using the expressions given by (10.14) and (10.16) for the direct and 
exchange amplitudes, and substituting the expression for w through the 
resolvent G into these expressions, we have from (10.1) 


Wri.) = ene Wolr1)— G(E+ Eo + 10) [Vo(r2) + V(r,, r)] em Wolr,) 
(16.1) 


Since we are interested in the forward scattering amplitude, we shall set 
v =n. This yields 


Fafor a | Mr User) lVolr) + V0.1] 


x G(E+ Ey +i0)[Vo(r5) +V(rj, 75)] 2 Wo(r) dr, dri (16.2) 
and 


1 —ikv.r’ / t t , 
9 = 90+ 5 | eH USMC + VC 3) 
x G(E+ Ey +10) [Vo(rg) + Vir, 75) ] et"? Wo(r)) dri dr (16.3) 


where 
fo=— 5 | [Vales + V4. rN Mole DP art dr (16.4) 
T 
oe =| (Vor) + Very, ro) ] er TO Wor )Wo(ra) dr dr, (16.5) 


The functions fo and gy are real and can readily be given in an explicit form. 
For f—fo and g—gy we have expressions analogous to (4.16): 


Im f(E’ 
Ref =fo +- to acta i a ee (16.6) 
and 
l Im g(E’) C; 
= —Pf | dE’——— - d 16.7 
ee 20ene fe E-E + E+|E) er 
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where 
d,;= | ert (rs) [Vor +V (ry, ro) ]uj(r1, r2) dry dre 

x [mmole [VOCE VL MGCL rDdr{ dry (16.8) 
Cc; = | dr, dre“ Who(ry) [Vo(r2) + V(ry, re) ]uj(ry, r2) 


xf ari are "Wor )LVolr + MM U(r.) 06.9) 


and # represents the prineipal value of the integral, while x; = J 2|E;|. 
The singlet ground state of the negative hydrogen atom (binding energy 

~ (0-7 eV) is well known at the present time. There are grounds for supposing 

that H~ does not have exeited states. We shall continue our analysis on the 

assumption that there is only one singlet state. In this state the funetion 

u,(r,, f2) is symmetric and therefore we have from (16.9) and (16.8) 

2 


et Wo(rs) [Volri) + V(r," 2) Jur, 12) dr, dr, (16.10) 


Substituting the symmetrie and antisymmetric combinations of f and g 
and using (10.45) we have 


k'a ©) ag 2c 
Req* = qo +73 =39 |e ~ E+E (16.11) 
1 


1 
SS 
: 2n 


kia (E’) 
E'—E 


Req” =qo + efor’ (16.12) 


As was to be expected, the singlet ae state of the negative hydrogen ion 
has an effeet on the singlet scattering amplitude g*~ but has no effeet on the 
triplet amplitude g°. For comparison with experiment one must use a 
formula for the eross-seetion <a> which is equal to 347 +407. It is readily 
shown that 


k 
7 <9) = Im(f—49) (16.13) 
TU 
The dispersion relation for f—4g is 


Uk yy c 


“PE ~XE+ ED (16.14) 


Re(J—39) = (S-4ao + a3 fae 


At zero energy (16.14) yields the ésitowine expression for the seattering 
lengths: 


fa” t+4at* = (f—4g)o —- sa oof ay (16.15) 


2|E| 
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An interesting example of the application of (16.15) is the analysis of the 
two contradictory experiments on the scattering of electrons by hydrogen 
atoms [53, 54] which is given in [52]. The analysis is based on the comparison 


of the values of i <a» dk from (16.15) and the same quantity determined 
1) 


from experimental data. 
First we note that because a* and a™ are positive we have immediately 
the inequality 
ce @) 


c 
[ <o> dk < 2n? (aE a $40] (16.16) 


0 


If we now calculate fo and gy from (16.4) and (16.5) we have 
fo-390 =-2 (16.17) 


The quantity c/2|E,| cannot be calculated exactly since the wave function 
for the ground state of the H™ ion is not known. It is pointed out in [51] 
that the approximatc values available for this quantity vary by about 20%, 
depending upon which approximation is used for the wave function of H™. 
Using calculatcd values for the quantities on the right-hand side of (16.16) 
wc have 


{ <a) dk < 7-52 (16.18) 


On substituting the values of a* and a™ calculated in (46, 48] into (16.15) 
we obtain 


| <a) dk = 4-77 (16.19) 


oO 
Let us now return to the determination of | <a> dk with the aid of 


O 
experimental data. This was done in [51] by the following procedure. At 
low energies (<S 10 eV) use was made of the experimental data on elastic 
scattering of clectrons by hydrogen atoms. At high energies (2 500 eV) 
the determination was based on the values calculated on the Born approxi- 
mation. The cross-section <a> in the intermediate region was dctermined 


ie @) 
by interpolation. The uncertainty in | <a> dk due to interpolation errors 


0 
can be estimated quite easily, and the final result is 


foe) 


{ 
= | <a) dk = 10:3+0°5 (16.20) 
0 
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when we use the results of [53], and 


1 
2n 


3 | <> dk = 5-1+0°5 (16.21) 
0 
when we use the results of [54]. 


As can be seen, (16.20) lies well above (16.18), whereas (16.21) is in 
agreement with (16.18) and (16.19). 


17. ANGULAR DISTRIBUTION OF ELECTRONS AND THE SCATTERING 
CROSS-SECTION. COMPARISON WITH EXPERIMENT 
In this section we shall considcr the differential scattering cross-sections 
averaged over the spins and calculated from the phase shifts 6¢, df and 
6; with allowance for the 1s—2s—2p coupling for four values of the 
energy [45]. The corresponding curves are shown in Fig. 19. The points 


in this figure represent the experimental data on the angular distribution 
as reported in [54, 55]. 


2O 
oN 
L 10 
0°29 0 000 O% 3: 8eV 
1-O 5 
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4 
Fic. 19. The differential scattering cross- Fic. 20. The elastic scattering cross- 
section section 
i Ecole 3 da~ <a> = jo_ + ho, 
dQ 4d2  4adQ in units of 2a2 


in units of zap? 


One notes the characteristic form of the angular distribution of the 
scattered clectrons. At low energics (~ 3eV) forward and backward 
scattering predominate. This effect is duc to the exchange interaction. 
As a result of exchange the phase shifts 5s and 5; are much greater than 
the values 6, and 6, which are obtained without allowing for the exchange. 
Figure 20 shows the total clastic scattering cross-scction (or morc preciscly 
the sum of the three first partial cross-sections). The various experimental 
points are given for comparison [45]. 


Chapter 5 


EXCITATION OF THE HYDROGEN ATOM BY 
ELECTRON IMPACT 


18. APPROXIMATE METHOD FOR THE CALCULATION OF THE 
EFFECTIVE EXCITATION CROSS-SECTION 


THE various approximate methods can conveniently be discussed on the 
basis of the system of equations given by (10.55). If we confine our attention 
to direct transitions between the initial and final states and neglect virtual 
transitions to intermediate states and also the coupling between different 
real transitions in the case of multichannel processes, then the system given 
by (10.55) reduces to the following two equations: 


AF ,+(k2 —2Vq)Fg—2 | Waa Fade" = 2Vq9 Fo +2 | Wao Fo dr (18.1) 


AFy+(k2—2Vo,)Fo —2 i WyoFodr’ = 2(Voa Fat WF, dr’ (18.2) 


The quantities V,,, Wi. Voo, Woo and W,,, fall off exponentially as r > o. 
The quantity V,, falls off as r~'~' when r > oo (except in the case of tran- 
sitions between s-states when V,, is also found to be exponential). If we 
retain only the slowly decreasing interactions represented by V,, and neglect 
all others we obtain 


AF +k2F, =2V,oF 9 (18.3) 
AFotkeFo = 2v,F, (18.4) 


The simplest and crudest approximation is the Born approximation in 
which V, is regarded as a small perturbation. The term V5, F, on the right- 
hand side of (18.4) is then neglected and a plane wave is substituted for Fo 
on the right-hand side of (18.3). As a result we havet 


AF, +k2F, = 2V io elhov" (18.5) 


The amplitude g, on the Born approximation is 


dor =— 2 e~ (kan — kov)r’ V0 dr’ (18.6) 
2 

+ There exists also the so-called Born-Oppenheimer approximation in which (18.5) 
is augmented by the exchange term f Wao e'*e-r’ dr’ where the function Fo is replaced by 
a plane wave. However, numerical calculations have shown that this leads to a con- 
siderable overestimate of the cross-section. In many cases the cross-section calculated on 
this approximation exceeds the theoretical limit by a large factor. A modification which 

removes much of the discrepancy is described in [105]. 
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Let us introduce the vector 
K=k,n—kopv (18.7) 
which represents the momentum transfer. In this notation we have 
| -iK.r 
i.== a e Vio dr (18.8) 


We can now use the expression for V,,, given by (10.20) to transform (18.8) 
to the form (cf (5.25)) 


2 -iK.r 2 —i 
doa = K2 ee Wo dr = ZAG we) (18.9) 
The total excitation cross-section is given by 


Cor = an | |q_|” sin 6 dO 
0 
0 


The integration with respect to 0 may be replaced by integration with 
respect to K. From (18.2) we have 


K? = k?+ki—2k,k,cos0 (18.10) 
Since 
k? = k6—2AE (18.11) 
where 
AE = £5, £5 (18.12) 
it follows that 
K* = 2(AE—k,k, cos 0) (18.13) 
and hence 
KdK = kok, sin @ dé (18.14) 
Therefore, 
kotk, : 
cua | Me™ al (18.15) 
ko—ke 


When the momentum transfer is small, the matrix element (e'*-”),, can be 
expanded into a series in powers of AK. The first term will vanish because 
the wave functions of the excited and ground states of the atom are 
orthogonal. For optically allowed transitions (e'*-")., ~ K(r).9, which 
may be referred to as the dipole approximation. The quantity (r),, in its 
turn is related to the oscillator strength for the transition 0 > a [56]. In 
general, the dipole approximation for (e'*-"),. cannot be used in the cntire 
range of integration in (18.15) since at the uppcr limit the quantity A may 
not, in fact, be small [57]. However, the cxpression given by (18.15) may be 
modified by introducing instead of the uppcr limit A4)+k,, an cffective 
number k’ which is chosen so that the dipole approximation can be used in 
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the entire range of integration. The number k’ can, for example, be deter- 
mined from the condition that for sufficiently high energies the magnitude 
of the cross-section o calculated from (18.15) (without expansion in the 
powers of K) should be practically identical with the value of o calculated 
on the dipole approximation [58]. On substituting the dipole approximation 
for (e'*-") , into (18.15) we have 


(18.16) 


Since, on the dipole approximation, the cross-section o is a logarithmic 
function of k’, the latter need not be determined to a high degree of accuracy. 
On the dipole approximation, therefore, the excitation cross-section do, 
can be expressed in terms of the oscillator strength. 

On the Born approximation the expression for the partial amplitude 
is obtained with the aid of the partial-wave expansion (18.8). On the total 
angular momentum representation we then have 


dw =—2 | ailkor)Y wathy (kyr) dr (18.17) 
0 


We note that even at low energies (near the threshold), the amplitude 
(18.17) for a given transition is not a monotonically decreasing function of 
the angular momentum /, of the incident electron. In contrast to elastic 
scattering, when the s-wave has a dominating effect at low energies while 
the contribution due to partial waves with / # 0 is small and decreases 
rapidly with /, it turns out that in the present case the dominant partial 
wave is that corresponding to a value of /, which is not the smallest value. 
This phenomenon was established, for example, in [56] as a result of calcu- 
lations of the excitation of p-levels in He. The reason why a partial wave 
corresponding to an angular momentum /, may give a bigger contribution 
than waves with /,—1 and /, +1 can be deduced by analysing the properties 
of the integrand in (18.17). For example, when the atom undergoes a transi- 
tion from the s-state to the p-state, the matrix element ¥,,-, vanishes at 
r = 0 and r— o, while its modulus |”\,,-;| has a maximum for r of the 
order of a few atomic units of length. The functions w,, and w,; vanish at 
r = 0 and then monotonically increase with increasing r up to the first 
extremum (the position of the extremum depends on 7, and k). As r Is 
increased still further, the two functions exhibit an oscillatory behaviour. 
At r = 0, and in the oscillatory region, u,, and w,; contribute very little to 
the integral. In order that the value of the integral as a whole should be 
sufficiently large, itis necessary that the first maxima of w,, and u,; should 
lic in the neighbourhood of the maximum of |¥%,,-,]._ This is achieved only 
at a particular optimum value of /,. 

A modification of the Born approximation has been given in [59]. In 
this work the idea was to calculate the reaction matrix K given by (12.11) 
on the Born approximation, i.e. assume that 


K =—k'ttB Kitt (18.18) 
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where By is obtained from B by replacing the functions x,,- in (10.72) by 
the terms u,d,,- outside the integral. The exchange interaction was 
neglected. We shall denote the approximate reaction matrix given by (18.18) 
by K©® and the approximate amplitude by q‘®. According to (12.9): 


q® = k~FK%%(] — 1K)" 1k7F (18.19) 


In point of fact this approximation is a generalization of the formula 
given by (5.18) which was derived for elastic scattering by a centre of force. 
If the matrix elements of K‘® are small, then q’? ~ k7~?K'™k7? and we 
have the Born approximation. However, if the interaction is sufficiently 
strong, the Born amplitude given by (18.17) may exceed the limit set by the 
unitarity relation for the accurate amplitude, while the expression given 
by (18.19) will never exceed this limit. Therefore, the latter expression may 
be looked upon as a reasonable extrapolation of the Born approximation 
to the region of strong interaction. It may be used to allow, at least in part, 
for the effect of coupling of the type a > a’, a > a’ ~ a’ and also the com- 
petition between the a> a’ and a— a” transitions. In this method it is 
no longer possible to sum over all the partial waves by reducing the entire 
sum to a simple expression of the form of (18.8) as on the Born method. 
However, in some cases it is sufficient to take into account a few partial 
waves, and then the summation problem does not arise. When a large 
number of partial waves enter into the calculations it is occasionally possible 
to use the fact that, beginning with a certain value /,, the amplitude calcu- 
lated from (18.19) is practically identical with that found on the Born 
approximation. We then have the approximate relation 


lo 
g=0"+) (6-9/7) 
1=0 


When /, is very large, this method is not particularly productive and 
other methods of estimating this part of the sum must be found. In [60] 
the specific case of the excitation of the hydrogen atom is used to discuss 
in detail the various possible methods of summing over all the partial waves. 

The effect of the atomic field and the exchange interaction on the electron 
can be calculated on the distorted wave approximation. In this method, 
the quantities V,,, W,,, Vo9 and Wo, are retaincd in the system given by 
(18.1) and (18.2). Moreover, V,, and W,, are retained in (18.1), but in 


(18.2) they are neglected. In this way one obtains the uncoupled equations 
AF, +(k2—2V,,)F,—2 i WF, dr’ = Vig Fo+2 { W.oF dr’ (18.20) 


AF) +(k3—2Voo)Fo—2 | WooFodr’ =0 (18.21) 
The asymptotic forms of F, and Fo are 
Fie gitar yg Oo 
i ’ (18.22) 


Fa = doa 
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Let G,(r, r') be the Green’s function satisfying the equation 
AG, +(k2—2V,,)G,—2 | W.(r,r")G,(r", r’)dr" = 5(r’—r) (18.23) 


and having the following asymptotic form for r — oo 


1 ikar 
G,~ ——F {r') (18.24) 
4n r 
where #, is a function to be determined. We can then write 
Po | G(r, r’) (Vo Fo+ i W.oF o dr") dr’ (18.25) 
In view of (18.24) we have for large r 
1 ikar 
Fy~- 5" | F.( Voor +] Wao Fo dr") de (18.26) 
nm or 
and hence 
1 
doe =~ 52 | Fal YoFe + | WaoFod") dr (18.27) 
Tt 


In order to elucidate the significance of the function F, let us express G 
in terms of Gp, which is given by 
1 eikalr—r’| 


On jr—r'| 
The connection between G and Go is contained tn the integral equation 


G=G, -| Vig’ G(r, vr Go(r", vr’) dr” 
-| Wir’, G(r, r")Go(r", r') dr’ dr” (18.29) 


Gyo= (18.28) 


For large r 
Gog~-—-—e hr (18.30) 


Assuming that, for large r, G is given by (18.24) we have 
F (vr) = e7 ikem’ _ 9 | V_.GoF dr’ —2 | WG ¥ dr'dr" (18.31) 


We thus obtain an integral equation which is equivalent to the integro- 
differential equation 


AF, +(k2—2V,,)F,—-2 | W,,F, dr =0 (18.32) 


with the asymptotic behaviour 
ikgr 
e x 


F . ~ eW tkanr Qax— (18.33) 


As can be seen from these formulae the function ¥ describes the elastic 

scattering of an electron by the excited atom where prior to the collision 

the electron has the momentum k,n. The distorted wave method involves 
E.A.C. 8 
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the resolution of F and ¥ into partial waves. The partial amplitude g,,. 
in the distorted wave approximation is proportional to the integral 


[A. (7) >. vir fy, Lr) ak | Wy. vill’, r)fy, ir’) dr’ dr 


0 0 


where f,., describes the elastic scattcring with angular momentum /, by 
an atom in the ground state, where /,,, represents the elastic scattering 
with angular momentum /, by an cxcited atom in the state m, /,.. As in the 
simplest casc of the Born approximation, g,y 1s not a monotonic function 
of the angular momentum /7, of the incident electron. There are spccial 
values of /, which play a dominant role. These values arc somewhat different 
from those which enter into the Born approximation, since the zeros and the 
extrema of the functions f,,, are displaced relative to the corresponding 
points for the functions w;,(r). 

One occasionally uses a simplified version of the distortcd wave approxi- 
mation in which the exchange interaction 1s neglected, so that (18.20) and 
(18.21) are replaced by 


AF, +(kz —2V,,)F, = 2Vi0F 9 
AF o +(k* —2Vo0)F 6 = O 


We note that unless factors other than the distortion arc also taken into 
account (for example, exchange and coupling through intermediate states), 
the final result may not in fact constitute an improvement on the Born 
approximation, but may on the contrary lead to an cven larger discrepancy 
between the approximate and true amplitudes. 

A review of the various versions of thcse methods from the point of 
view of variational principles is given in [17]. 

When all the terms in (18.1)-(18.2) are included, the corresponding 
approximation 1s known as the strong coupling approximation. However, 
this approximation takcs into account only the strong coupling O—a and 
ignores the effcct of O—a’—«” ctc. coupling. Nevertheless, the cffect of the 
latter may be appreciable, so that it is better to consider a system of a number 
of equations which take into account the coupling within a group of states. 
Both the distorted wave and the strong coupling approximation involve 
extensive numerical work which is usually carricd out by computers. A 
detailed discussion of the various problems associated with the sclection 
of the best computational scheme is given in [61]. 

In addition to numcrical calculations, usc can also be made of variational 
methods. In their simplest form the variational methods are employed 
within the framework of the distorted wave approximation. As was shown 
above, in the distorted wave approximation the transition amplitude is 
expressed in terms of functions describing clastic scattcring by the atom 
in the ground and cxcited statcs. These functions may be represented by 
approximate analytic functions containing paramcters to be varied. Therc- 
forc, the application of the variational method rcduces to the determination 


(18.34) 
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of the optimum parameters of trial functions describing elastic scattering— 
a problem with which we are already familiar [17]. If the excited state is 
an s-State, the trial function will in general have the same structure as in 
the case of elastic scattering in the ground state. But if the excited state has 
a non-zero angular momentum, then the variational calculation becomes 
so involved that it is much easier to carry out a numerical integration of 
the equation. 

The application of variational methods becomes more difficult in the 
case of the strong coupling approximation. Comparisons of different 
methods, and also the comparison of variational calculations with the exact 
solution (using a model which admits of an analytical solution on the strong 
coupling approximation) are given in [62]. 

Below we give the results of some calculations concerned with the excita- 
tion of the hydrogen atom to different states; these will be compared with 
experimental results. 


19. EXCITATION OF THE 25-STATE 


The excitation of the 2s-state of the hydrogen atom has been extensively 
investigated. It is the testing ground for the various methods which are 
available for such calculations. Thus, calculations have been carried out 
using the modified Born approximation [60], the distorted wave approxi- 
mation [63], the strong coupling approximation (1s—2s) [64], and finally 
calculations with allowance for the 1s—2s—2p coupling [61, 65]. 

It is evident from Table 1 that the latter calculation involves the solution 
of three equations for L = 0 and four equations for each of the angular 
momenta L > 0. The calculations yield the excitation cross-section for the 
2s-state and for the 2p-state. Partial cross-sections for 0 < L < 4 are cal- 
culated in [65] and for 0 < L < 61n [61]. The total cross-section can be 
obtained by adding the partial cross-sections for large L, calculated 
on the modified Born approximation. Figure 21, which is based on the 
data reported in [65], shows the excitation cross-section for the 2s-state 
averaged over the spins and computed on different approximations. The 
partial cross-sections for angular momenta L between O and 2 are shown 
separately for the triplet and singlet states in Fig. 22. Inspection of these 
figures shows that the results obtained when the 1s—2s coupling is taken 
into account are not very different from those obtained when ls—2s—2p 
coupling is allowed for. The addition of the 2p coupling has little effect 
on the excitation cross-section in the triplet state, except for the region 
near the threshold; but it has a more pronounced effect in the singlet state. 
Moreover, the cross-section near threshold is increased when 2p is taken 
into account, whereas away from the threshold it is reduced. The partial 
cross-sections for large L exhibit a greater relative change when the coupling 
with 2p is added than do the cross-sections corresponding to small L. 

The excitation cross-section for the 2s-state has been determined experi- 
mentally [66]. Direct measurements of the 1s—2s transition cross-section 

g—2 
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Fic. 21. Cross-section for the excitation of the 2s-state (in units of mao”): solid curves— 

distorted wave approximation (sum of partial cross-sections for L = 0, J, 2); crosses— 

with ls — 2s — 2p coupling (all partial waves); open circles—strong 1s — 2s coupling 

approximation (sum over L = 0,1, 2); triangles—with ls — 2s — 2p coupling (sum 
over L = 0, J, 2) 
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Fic. 22. Partial cross-sections for the excitation of the 2s-state in units of za o?: (a)—sym- 
metric case, (6)—antisymmetric case (solid curves—distorted wave approximation; open 
circles—with strong 1s — 2s coupling; crosses—with ls — 2s — 2p coupling) 
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are not readily obtainable, and therefore in the experiment metastable 
hydrogen atoms in the 2s-state were subjected to the action of an electric 
field which induced their transitions to the 2p-state. The radiation due to 
the 2p— 1s transition was measured, and from the intensity of this radiation 
an estimate was made of the magnitude of o(is—2p). In order to obtain 
the absolute magnitude of o it was necessary to introduce a calibration. 
For example, it was assumed in [66] that at 45eV the cross-section o is 
practically equal to the cross-section obtained on the Born approximation. 
Another method of calibration was used in [67]. 

The results reported in [66] and [67] differ by a factor of about 2. The 
reason for the discrepancy is obscure, and therefore the comparison between 
theory and experiment is somewhat preliminary in nature. Nevertheless, 
there is a good agreement between the calculations of [65] and the experiment 
described in [66]. Figure 23 shows the experimental curve and the theoretical 
results from [65]. 


E (ev) 


Fic. 23. Calculated cross-sections for the excitation of the 2s-state and experimental 
results (solid curve—experimental; crosses—calculation including 1s — 2s — 2p coupling. 
The cross-section is given in units of 2a? 


20. EXCITATION OF THE 2p-STATE 


The cross-section for the excitation of the 2p-state has been calculated 
on the Born approximation [60] by the distorted wave method [68] and with 
allowance for the 1s—2s—2p coupling [65]. As in the calculation of the 
excitation of the 2s-state, the partial cross-sections with allowance for the 
ls—2s—2p coupling were determined for 0 < LZ < 6 and were then aug- 
mented by the cross-sections calculated on the Born approximation for 
L> 6. 

Figure 24 shows the partial cross-sections for L = 0,1,2 calculated 
by different methods. The formula given by (11.13) was used to determine 
the amplitudes for the excitation of the 2p, and 2p,, states. The excitation 
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Fic. 24. Partial cross-sections for the excitation of the 2p-state (in units of zap7): 

(a) symmetric case, (b) antisymmetric case (solid curves—distorted wave approximation; 

open circles—with strong coupling Is — 2s; crosses—calculation with coupling. For 
L = 1, 2, the cross-sections arc reduced, taking /, = L — 1) 
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Fic. 25. Polarization of the radiation Fic. 26. Comparison of measured and 
cmitted as a result of the 2p — Is transition calculated values for ao): 1—experiment; 
(curve — calculations; points — experi- 2—calculated with 1s — 2s — 2p coupling; 

mental) 3—Born approxunation (all in units of a9?) 


cross-sections for these states, o(2p,) and o(2pp), were then caleulated from 
these amplitudes. In addition the tollowing quantities were calculated: 
1. The total cross-section for the excitation of the 2p-state: 


o(2p) = o(2po)+ 20(2p4) (20.1) 

2. The cross-section a, determined experimentally from the number 

of photons emitted by the excited atoms in the direetion perpendicular to 
the direction of the electron beam. Thus cross-section is given [69] by : 


a, =0-9180(2p,)+0°246a(2po) (20.2) 
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3. The polarization of the radiation emitted at right angles to the 
direction of the electron beam [69]: 


_ _ o(2p;) I 
os E ror) 5(2p,) ee) 
7+11 5(2p.) 


The calculated cross-sections are given in Table 5 in units of za?. The 
calculated values of P and o are compared with experimental values in Figs. 
25 and 26. 


TABLE 5 

k o(2p,) o(2po) o(2p) Oy P 

0-9 0-:0683 0:2234 0:3599 0:3853 0:2009 
1 0:0997 0:6005 0:7909 0:8820 0:2835 
1:1 0:1395 0:8149 1:0939 1:2047 0-2799 
1:2 0:1631 0:8539 1°1801 1 -2934 02667 
1°5 0:-2244 0:6423 1:0911 1-1596 0-1800 
2 0:2395 0:3927 0:8716 0:8968 0-0855 


21. EXCITATION OF STATES WITH 7 > 2 


The excitation of states with n = 2, 3 and 4 has been calculated on the 
Born approximation in [70]. The effective cross-sections summed over 
J and yp for given 7 are given in Table 6. 


TABLE 6 
E Cross-sections (in units of zap”) for the transitions 
(eV) [=D i—3 i—4 i—S5 i—6 
13:5 1:2818 0-1787 0:0509 0-:0202 0-0092 
19-5 1°5354 0-2782 0-:0999 0:0476 0-:0265 
26:5 1°4993 0-2798 0-1021 0:0490 0:0274 
34°8 13886 0:2600 0-0951 0:0457 0:0256 


The interpolation method was used to determine the excitation cross- 
sections for states with mn = 5 and n = 6. This method ts based on the fact 
that for large quantum numbers » an energy change of dE is related to the 
change dn in the quantum number by 


dn~ nrdE (21.1) 
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The excitation cross-section of an atom in the interval dn near n is 

do =o, dn (21.2 
or, on substituting (21.1) 

do =o,n°dE (21.3) 


In the case of ionization of the atom, the probability that the ejected 
electron will be found in an energy interval de may be written in the form 


do =0o,deé (21.4) 
It is therefore expected that n°o, will be a continuous function of 1 which 


goes smoothly into ¢, when n > o0.° Figure 27 shows n°a, as a function 
of the energy of the atomic electron. Negative values correspond to discrete 


-O2 -Oll O O'| O-2 
E (eV) 


Fic. 27. Dependence of n°a, on the energy of the atomic electron (¢ in units of 2do7). 
The figures against the curves indicate the energy of the incident electron in eV) 


levels, positive values to ionization. In the region of ionization, o, was 
calculated on the Born approximation. The broken curves indicate the 
interpolation region. The excitation cross-sections calculated in this way 
for states with n = 5 and n = 6 are also given in Table 6, but we only 
cover the low energy region. In [70] the cross-sections are calculated up to 
energies of about 1000 eV. 


22. TRANSITIONS OF AN EXCITED ATOM INTO HIGHER LEVELS AS A 
RESULT OF A COLLISION WITH AN ELECTRON 


In the preceding sections we considered transitions of an atom from 
the ground to the excited state produced by electron impact. In some cases 
it is necessary to know the cross-sections for the excitation of atoms from 
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excited states to higher excited states. A number of such transitions have 
been treated on the Born approximation. 

The cross-section for the transition from the 2s-state to all the levels 
with m between 3 and 10 summed over / and yp are calculated in [71]. The 
calculations were based on the general formula for values of n up to n = 7, 
whereas for n = 8, 9 and 10 use was made of the interpolation method 
described in § 21 in connection with the ls—n transition. The calculations 
in [71] are concerned with the energy range between 4°85 and 120 eV. 
Some of the results are given in Table 7. As is expected, the calculated 
cross-sections are much greater than the cross-sections for the corresponding 
transitions from the ground state. 


TABLE 7 

E Cross-sections (in units of zao?) for the transitions 

(OV) “96-34 O94 Is 5 DG eT Weak Oe 9 Ds — 10 
4°85 65-019 12-330 4-658 2°312 1°331 0-843 0:569 0-404 
8-65 49-441 9-454 3:612 1-808 1-047 0:664 0-450 0-320 
13-5 37:°667 7+150 2°726 1-364 0-790 0-499 0-339 0:241 
26°4 23:703 4-446 1-687 0-842 0-487 0-309 0-209 0-149 
34°6 19-488 3°640 1-379 0-687 0-397 0-252 0-170 0-121 


Transitions from the 2p) and 2p,, states to levels with nm between 3 and 
10 were calculated using an analogous method in [72] in the energy range 
between 4-85 and about 120 eV. Some of the results are given in Tables 8a 
and 8b. It is possible to use these data to determine the excitation cross- 
section from the 2p-state averaged over up: 


o(2p) = $0(2po) + $0(2p +1) (22.1) 
TABLE 8A 
E Cross-sections (in units of zao”) for the transitions 
(CV) 2p, -—3 2po—4 2po —5 2po — 6 2po —7 2po — 8 2po — 9 2po — 10 
4:85 74-837 14:-426 5:259 2-499 1-409 0-892 0-601 0-424 
8:65 64:529 12:323 4:749 2-347 °1:348 0-861 0-586 0-417 
13:5 52:482 9-826 3:790 1-859 1:072 0-684 0:461 0-328 
26:4 35-403 6:352 2-408 1:176 0-677 40:432 0-288 0-207 
34-6 29°727 5:246 1:974 0:963 0°553 0:348 0:234 0-165 


a i a Sr I a IT a 
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TABLE 8B 
E Cross-sections (in units of a9) for the transitions 
(VY). “Opa Dp ed. On 5 Ds 6 Op 7 iy = 9p — 10 
4-85 76°S15 13°346 4-995 2512 1-414 0-890 0:600 0-424 
8:65 57°44 10:692 4-052 2:040 1-166 0-740 0-505 0-358 
13:5 45:05 8-181 3:077 1:537 0-883 0:560 0-378 0-270 
26:4 29:278 5102. 1°891 0-935 0-537 0-342 0-229 0-160 
34-6 24°313 4-173 1-538 0-758 0-435 0-268 0-180 0-128 

0-8 | 2 345 810 20 3040 
E (eV) 

Fic. 28. Cross-sections for the transitions 3s — 4p (curve 1), 3p — 4d (curve 2) and 


3d — 4f (curve 3). The cross-sections are given in units of zap? 


The cross-sections for the following transitions were calculated in a number 
of papers [73, 74, 75] (see Figs. 28-34): 

3s — (4s, 4p, 4d, 4f, Sp) 

3p — (4s, 4p, 4d, 4f, 5d) 

3d — (4s, 4p, 4d, 4f, 5f) 

4s — (Sp, 6p) 

4p 5d 

4d 3f 

4f > (5g, 69) 

S(s, pdf, g) > O(p, df, 9, h) 


+ The broken curves in these figures represent the dipole approximation. 
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O06: { 2 345 810 
E (eV) 


Fic. 29. Cross-sections for the transitions 3s — 4s (curve 1), 4s — 4d (curve 2), 3p — 4p 
(curve 3), 3p — 4f (curve 4), and 3d — 4d (curve 5). The cross-sections are given in units 
of TMAo” 


0-6! 1 2 345 B10 2030 


E (eV) 


Fic. 30. The cross-sections (in units of 2a *) for the transitions 3s — 4f (curve 1), 
3p — 4s (curve 2), 3d — 4s (curve 3), and 3d — 4d (curve 4). The scale for curve 1 is 
reduced by a factor of 30 


1 2 3°45 8 10 [5 
E (eV) 
Fic. 31. The cross-sections (in units of zao”) for the transitions 3s — Sp (curve 1), 
3p — Sd (curve 2), and 3d — Sf (curve 3) 


05 :10 (82 3 5 810 1520 40 


E (eV) 


Fic. 32. Cross-sections for the transitions 4s — Sf (curve 1), 4p — Sd (curve 5), 4d — Sf 
(curve 3), and 4f — Sg (curve 4). The cross-sections are given in units of 2a? 


500 
400 


300 


200 
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E (eV) 


Fic. 33. Cross-Sections for the transitions 4s — 6p (curve 1), and 4f— 6g (curve 2). 
The cross-sections are given in units of zap? 
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E (eV) 


Fic. 34. Cross-sections for the transitions 5s — 6p (curve 1), 5p — 6d (curve 2), 5d — 6f 
(curve 3), 5f— 6g (curve 4), and 5g — 6h (curve 5). The cross-sections are given in 
units Of 7a? 


As can be seen from Figs. 28-34, the cross-sections for these transitions 
are much greater than the cross-sections for the corresponding transitions 
from the ground state. The cross-sections for optically allowed transitions 
involving an increase in / are greater than the cross-sections for optically 
allowed transitions involving a decrease in /. 


Chapter 6 


THE COLLISION OF AN ELECTRON WITH A HELIUM ION 


23. THE ELASTIC SCATTERING OF AN ELECTRON BY A HELIUM ION 
AND THE EXCITATION OF THE 2s STATE OF THE HELIUM ION 


THE elastic scattering of an electron by a helium ion differs from the elastic 
scattering of an electron by a hydrogen atom owing to the presence of the 
unscreened Coulomb field oc I/r acting on the scattered clectron. This 
Coulomb field modifies very significantly the scattering process and the 
form of the wave function, and its effect 1s not simply additive to that due 
to the hydrogen atom (the two effects are additive only on the Born approxi- 
mation). By analogy with the seattering of an electron by a hydrogen atom 
we shall begin with the approximation on which the helium ion is regarded 
as a centre of foree. We can then use the theory given at the end of § 2 and 
the problem is reduced to the determination of the phases ,1,. 

In [76] the phase jtg 1s calculated by a variational method, and also 
(at two energies) by a numerical integration of the equation 


d? 2 
can Gi age aw) f= 0 (23.1) 
dr r 
The formula for w 1s 
1 2 
ee Md — dr 
r lr—r | 


where w = 8!x~t exp (—2r) is the wave function for the ground state of 
the He* ion, so that 


1 
Woo — (2 + = e74r (23.2) 
r 
In [76] the funetion fis normalized so that it has the asymptotic form 
f= F+tanwyc, r+ © (23.3) 


where F is a regular Coulomb function having the asymptotic form 
, 1 
F~sin (ir ~5 + , in 2kr+1) (23.4) 
and G is a non-regular function with the asymptotie form 


1 
G ~ cos (19 — - + L In 2kr+1) (23.5) 
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According to (2.68) and (2.71) the scattering amplitude g may be written 
as the sum 
erin 


q = ac + (*— I) (23.6) 


where qc is the Coulomb scattering amplitude. As a measure of the departure 
of the actual scattering from the purely Coulomb scattering, the authors 
of [76] used the quantity 


een eee (23.7) 
dc 
From (23.6) 
ae a, : 
= {1 + — ee? —(e?#_]) (23.8) 
Qc "Dy ik 
When exchange is taken into account, (23.1) is replaced by 
d’f 
7 oat + (k2 +22) f= 2] we r') f(r’) dr’ (23.9) 
In the present case 
t =2 eS ol k? I : 
Wri’) =+32re “'r'e a ee) (23.10) 


Equation (23.9) was also solved numerically for two values of the energy. 
Moreover, the phase p.* was found with the aid of a variational method. 
The trial function was taken in the form 


f=F+(at+be*")(1—e-*")G (23.11) 


The trial function for the electron plus ion system containing the inter- 
electron distance r,, was also employed. The quantities R, and R_, 
calculated with allowance for exchange, were averaged using the formula 


R=4R,+4R_ (23.12) 


The results of the variational calculation were found to be strongly dependent 
on the particular method employed, especially at low energies in the sym- 
metric case. This suggests that the trial function was not very satisfactory. 
The numerical values of jt are given in Table 9, and a comparison of the 
results of the numerical integration and the variational calculation is given 
in Table 10. The function R(@) at an energy of 48-84 eV is shown in Fig. 35. 
A calculation of the excitation of 2s and 2p states of the helium ion with 
allowance for the 1s—2s—2p coupling is reported in [77] and [77a]. The 
calculation is quite analogous to that for the hydrogen atom. The partial 
wave equations are similar to the equations derived in Section 11, except 
that the functions u, and v, are the corresponding continuous-spectrum 
wave functions for the Coulomb field. The presence of the unscreened 
Coulomb field has an important effect on the behaviour of the amplitud 
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Fic. 35. <A plot of the function R(@) for E = 48-84eV: 1—without exchange and 
polarization, 2—with exchange, 3—with exchange and polarization 


TABLE 9 
Lu ue 
E (eV) with aan Kohn’s Hulthén’s 
exchange method method 
3°28 0-571 0-959 —0:28] 0-419 
8-25 0-554 0-824 0:038 0-404 
15-74 0-531 0-807 0-073 0-365 
24-90 0-509 0-762 0-110 0-333 
48-4 0-460 0-645 0-240 0-315 
65°49 0-443 0-604 0-226 0-307 
TABLE 10 
rv p Z 
without exchange e . 
E oes ; 
Numeri- Numeri- 
(eV) elas Naor cal Kohn's Hulthén's cal ae 
eaie eae integra- method method integra- male 
Pee tion tion re 
48-4 0-468 0-460 0-288 0-240 0-315 0-648 0:645 


65°49 0-443 0:443 0-283 0:226 0:307 0-604 0-604 
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Fic. 35a. (Reproduced with permission from Burke e¢ al. [77] ) 


near the threshold. In particular, it is found that as the energy tends to the 
threshold value, the excitation amplitude increases without limit as k7?. 
This is a consequence of the behaviour of the normalization factor Co, 
given by (3.30), as k > o. It is the reason why the excitation cross- 
section o = 4n(k,/k,)|q|* remains finite when k = 0. 

Figure 35a shows the excitation cross-section for the 2s and 2p levels 
of Het based on the data reported in [77]. 


Chapter 7 
THE COLLISION OF AN ELECTRON WITH A HELIUM ATOM 


24. THE WAVE FUNCTION FOR A THREE-ELECTRON SYSTEM 


Tue form of the wave function with a given spin and the required symmetry 
with respect to the spin and coordinate variables for a system of three or 
more electrons is not so apparent as in the case of a two-electron system. 

There are a number of ways of constructing such wave functions in the 
form of a suitably symmetrized sum of products of coordinate functions 
and spin functions. The various methods of doing this are reviewed in [78] 
from the point of view of group theory. 

Consider the wave function ®(r,,0,,..., ry, Gy) Of an N-electron system. 
It should be antisymmetric with respect to the interchange of any pair 
(r;,0;) and (r,, o,). We shall require that it should be the eigenfunction of 
the square of the spin angular momentum, so that 


$?@ = s(s+1)® (24.1) 
V. A. Fock [79] showed that an antisymmetric function satisfying (24.1) can 
be expressed in terms of the coordinate function ‘P(r,, ..., rglfeais ++) PN) 


and the spin function y(¢,,..., 0,|6,44,; -.-,6y). Henceforth we shall write 
WU, ...,A|A+1,..., M) and y(J,..., Al//41,..., N) instead of V(r, |r.) and 
y(o,|oy), or simply V, x, whenever this does not introduce an ambiguity. 
The vertical line separating 1,...,A from A+1,..., N represents the division 
of the electrons into two swarms with opposite spin components. The 
number of electrons A which make up the swarm 1,..., & is related to the 
spin quantum number by 


N 
k=>=s (24.2) 


The functions ¥ and y have the following symmetry properties with respect 
to the exchange of the arguments: 

1. The function is antisymmetric with respect to the interchange of 
the arguments before or after the vertical line; this may be expressed in 
terms of the exchange operator P; ; in the form 


(+P; )¥ =0, ij<k or ij2k (24.3) 
2. WI, ..., AJA +1, ..., NM) exhibits cyclic symmetry, i.e. 


( — y Pa.) Y=0 (24.4) 
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3. The function y is symmetric with respect to the interchange of the 
arguments before or after the vertical line, i.e. 


(1-P,)y=0, ij<k or ij>k (24.5) 


Consider the expression 


—_ [ki(N—k)! 
® am oe » (—1)*P(A,, veey Aint vuey Ay) 
. Atscrer Ake 

KYA psssey An Ap eissens Ay) (24.6) 
where /,,..., Ay is a set of N unequal numbers from among the sequence 

L 2 we ay 
1,..., NV, and A is the parity of the permutation ( 4 ; ) In view 
“1 2 oeoe “N 


of (24.3)-(24.5), the function (24.6) is antisymmetric with respect to the 
interchange of r;, 0; and r,, o,. Next, using the well-known formula 


$2 = N—4N'+ 2 it (24.7) 


which relates the operator §? to the spin exchange operator P{, and bearing 
in mind (24.2) and (24.4), it can be shown [79] that ‘¥ is an eigenfunction 
of $* with the eigenvalue s(s+1). If x is chosen to be of the form 


NPs cate Py Ops ara when (24.8) 
where 


a=6,, and B=0,-, (24.9) 


then ® will also be an eigenfunction of S, with the eigenvalue s. The 
coefficient [((N—k)!k1/N!]* in (24.6) ensures that when and y are nor- 
malized to unity, the function ® is also normalized to unity. 

In the case of the three-electron system there are two possible spin states, 
corresponding to s = 4 and s = 3. According to (24.2) the coordinate 
function corresponding to the first case is ‘¥(1|23), while the wave function 
corresponding to the second case is ‘¥(|123). The properties expressed by 
(24.3) and (24.4) can then be written in the following form: 


(1+P,,,)¥((123)=0, o,f =1,2,3 (24.10) 
(1+P,3)'¥(1|23) = 0 (24.11) 
(1—P,,—P,3)¥(1|23) = 0 (24.12) 


The problem now arises as to how one can construct the two functions 
Y(|123) and ‘Y(1|23) so that they will have the properties given by 
(24.10)-(24.12). It turns out that when an arbitrary function of three variables 
(123) is symmetrized in accordance with Young’s scheme, we obtain a 
function having the properties of ‘¥(1]23) or ‘¥(|123) depending on the 
form of the scheme.t This can be seen most easily in the case illustrated 


+ This was shown for an arbitrary number of electrons in (80). 
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in Fig. 36(a) [57]. This scheme involves the action of the operator 
(1-—P,,)(1—P,3;)(1-—P,3) on the function (123). It is easy to see that the 
function 


@, = (1—P 42) —P,3)(1— P23) (123) (24.13) 
satisfies (24.10) since, for example, 
(1+P,.)U-P,.)=( — Pi) = 0 


It is also possible to carry out the symmetrization in accordance with the 
scheme illustrated in Fig. 36(b). This corresponds to the application of the 
operator (l1—P,3;)(1+ P,,). Consider the function 


Q2 = (1 —P3)(1 + P12) 90123) (24.14) 
It satisfies the condition given by (24.11) since 
(1+ P;)U—P23;) =0 
It also satisfies (24.12) since (as can be shown by direct calculation) 
(1—P,,—P,3)A—P23)14+ P;2) e(123) = 0 (24.15) 


, (a) (b) 


Fic. 36. Young's scheme: (a) for the function 9, = w(l1 23); (b) for the function 
G2 = y(I2 3) 


It follows that the functions ‘¥(1]23) and ‘¥(|123) can be obtained from an 
arbitrary function of three variables by means of the Young symmetrization 
scheme of Figs. 36(a) and 36(b) respectively, namely, Y(|123) = @, and 
¥(1123)= @s. 

The above analysis is concerned with the general case of a three-electron 
system. Since we shall be dealing with the collision of an electron with a 
helium atom, we shall consider three-electron functions (1 | 23) and ‘¥(| 123) 
of a special form, constructed by combining two-electron and one-electron 
functions. Since in the ground state the spin of the helium atom Is Zero, 
the only function which is of interest in practice is ‘¥(1|23). The function 
‘¥(|123) can only be useful in the case of transitions between triplet states 
of the helium atom as a result of a collision with an electron, and we shall 
therefore confine our attention to the function ‘¥(1 | 23). 

The function (1/23) can include two forms of two-electron functions, 
namely, symmetric w(1|2) and antisymmetric &(|12). This is a consequence 
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of the law of addition of angular momenta. The state of a three-electron 
system corresponding to spin 4 may be formed by combining a one-electron 
state with both a singlet and a triplet state of a two-electron system. This 
can also be expressed in Young’s terminology: the scheme shown in 
Fig. 36(b) may be constructed in two ways, as shown in Fig. 37. Accordingly, 
there are two expressions for ‘¥(1|23), namely, 


Y (1[23) = c,[W(1|2)F(3) — (13) F(2)] (24.16) 
Y (1/23) = c[W({12)F(3) — W({13)F(2) + 2 (]32)F(1)] (24.17) 


2}-qi} [21] 


(a) (b) 


Fic. 37. Two ways of constructing Young’s scheme for the function g(1|2 3): (a) addition 
of a cell to the scheme for yw (|2 3); (b) addition of a cell to the scheme for y(1|2) 


It can be easily verified that these expressions exhibit all the properties which 
the function ‘¥(1|23) should show. A general formula for the construction 
of (n+ 1)-electron functions from n-electron and 1-electron functions will 
be given in Chapter 8. The expressions given by (24.16) and (24.17) can be 
obtained from this formula as a special case. The coefficients c, and c, 
are the normalization coefficients. Suppose that F refers to the continuous 
spectrum and 1s normalized as a plane wave, i.c. 


i F*(k, r)F(k’, rn) dr = (22)95(k—k’) (24.18) 


Suppose further that w(r,|r.) and w(|r,, r,) refer to the discrete spectrum, 
and are normalized to unity: 


{ Werle)? ar, dr, =1, [Wri ra)? ary dr, =1 (24.19) 


If it is required that ‘V(r, |r2, r3) should also be normalized as a plane wave, 
1c: 
i ww. dr, dr, dr; = (2n)2d(k—k’) (24.20) 


then we must set 


6, (24.21) 


1 
6 
The connection between these coefficients and the normalization of the spin 
functions will be established in § 25. 
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25. EQUATIONS FOR THE WAVE FUNCTION OF THE 
SCATTERED ELECTRON 
Starting from equations (24.16) and (24.17) we can write the coordinate 
wave function for the electron plus helium atom system in the form 


1 
V(r |ro, r= BY Wars |e) Fairs) — War rs) F4i(r2)] 


1 
oS y LWadlry, r)F 42(r3)— Wari, r3)F 42(r2) 
/6 A 


+2W,(\rs, ro)F 4(r3)] (25.1) 


where wWy(r, |r.) and wWy(|r,,r2) are the coordinate functions respectively 
of the singlet and triplet statcs of the helium atom. In the simplest approxi- 
mation, they can be constructed from the one-electron functions: 


YW 4(l|2) = W(1)e(2)+¥(2)eU) 

Wa({12) = WD) e(2)— e(2)¥U) 
The functions F,,(r) deseribe the state of the seattered elcetron. The sub- 
script A represents all the quantum numbers charaeterizing the state of the 
atom, and a = |, 2 represents the two possible combinations of the spins 
of the atom and the incident eleetron. 


The equations for F,, and F,, may be obtained from relationships 
analogous to (10.24): 


[ valra|ra)(H—8)M(ry|ra, 43) dry dry = 0 (25.2) 
[ Vara r2)(H-S)¥(r,|ra, 73) dr, dry = 0 (25.3) 


where H is the energy operator for the electron plus helium atom system 
and is given by 


l I l 2 2 2 

Henk Hd ak, oe ee Se , Nx = |rimny| 
Fig Fpp- Fog Fy: Ty V3 

(25.4) 


The energy & of the system is made up of the energy of the atom prior to 
the collision and the energy of the incident elcetron. 

Let us substitute cquations (25.1) and (25.4) into (25.2) and (25.3). 
We shall take into account the symmetry propertics and normalization of 
war, |r.) and w,(lry,r2), and will assume that the latter two functions 
Satisfy the cquations 


2 2 l 
| edn (- 5A eae 8 © clear he a) Wa(ry|ro) dr, dr, =0 
hy Fa Ir, —r| 
(25.5) 
* 2 2 l 
wiry. ry) i. bA, << +A, i eee a Wa(|ry. r,) dr, dr, = 0, 
Py OF) Ir; —r,| 


A’ £A (25.6) 
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These equations are satisfied if w, and W,- are the exact atomic functions, 
but they may also be satisfied by approximate atomic functions. 
Let 


6 = tkgat Ena (25.7) 


where E,, is the energy of the atom and 3k4, is the energy of the incident 
electron.t After some simple rearrangement, we have the following system 
of equations for F,,, 


AF 4,(¥3)+ k4AaF salts) = 2» [Vaaralts)a' F grg(¥3) 
+ | Wa Ataa'(" 1; 13)F grq(¥'y) dr, | (25.8) 


where the potential energy matrix V,,,, is given by 


2 | 
Vawdlts) =~ 54a +2 Paarl) i (25.9) 
3 


Ir3—r1| 
with 


PAA‘ =| War tra) a(r|P2) dr, 
Paaa =| WH ri tear te) dre 


The exchange operators W, 444° are of the form 
Wasratists) =— | Wars |r2) (Hoya (ras) dr (25.11) 
Waara(tasts) =— V3 | WiCrsra)\(H-Oba(|rst2) dr, (25.12) 
Waaraa('asts) =— V3 | Wri, r2)(H—- OW ars|ra)dr, (25.13) 


Waar22"is'3) = — [uidn. r2)(H — Oh ars, r,)dr, (25.14) 


The expressions for the exchange operators can be transformed further. 
To begin with, consider the operator —4A, in H. Since this operator 
is self-adjoint, we can replace wW,(r, |r.)A; F4,(r,) in equations (25.11)-(25.14) 
by Fy,(r,A,Wy(r,|r.). A similar replacement may be carried out in the 
terms containing the function W,(r,|r.). Next, if W,(r, |r.) and WyC{ry, ro) 
are the exact wave functions for the helium atom, we can employ the con- 
ventional equations for them. As a result, the exchange operators assume 
the form 


I 2 
Wa araa'(" r3) aaa (44: + ice ~~ =| P aataa’(" 1; r3) 


31 1 
+ | Peer Qtaetan Ta) (25.15) 


Py2 


(25.10) 


+ When w, are the approximate functions, the quantity E,, is the mathematical 
expectation of the energy of the atom calculated from this function. 
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where 
Paarsa(tasts) = | WRQry|radbalra|rs) drs (25.16) 
paaaatists) = | WAC. ra Wal |r2, 43) dra (25.17) 
Paaialtasts) = V3 | Wa(ra|ra barat) drs (25.18) 


Paa2(145 13) = 3 { WaClrs, ro )Wa(ra|rs) dr, (25.19) 


Panay 425 %3) = Warr alr a|rs) (25.20) 
Paai2(his¥2sts) = V3 Wary |radba|ro. rs) (25.21) 
Paarai(his tats) = V3 W4( [ry ra Wara|rs) (25.22) 
Paara2h 1) ¥25%3) = We([ris ra ara, 3) (25.23) 


In practice the equations must be simplified. In addition to the errors 
which appear as a result of the simplification and which are also encountered 
in problems concerned with the hydrogen atom we have here one further 
source of difficulty, namely, the fact that the atomic wave functions are only 
approximate. This leads to errors in V and W. 

The equations (25.8) must be augmented by the asymptotic conditions. 
We shall assume that 

eikaar 
F 4,(r) Or CANO ig 4 aaa ee? r—> © (25.24) 


Let us now consider the asymptotic behaviour of the complete wave 
function for the system (including the spin variables). We return to (25.1) 
and consider any one of the terms representing two discrete states of the 
atom in the sum over A: 


I 
VA, ro, r3)= V3 [Wary |r2) Fars) — Wari |rs)Fai(r2)] 


1 
+ Raae r2)F 42(r3)— Wary. r3)F 42(r2) 
\/6 
+2 4(\r3, r2)F4o(r,)) (25.25) 


In addition to this term, the complete wave function also includes 
YP 4(r3|ry, rz) and ‘VY ,(r,[r3, ",). According to (24.6) we have 


1 
D(r, 01,12 02,1303) = pg LV Arar 2, 13)B 1 224 
N 


+P 4(r3|r1, r)B3 aya. +P 4(ro|rs, ry )B2%,%3) (25.26) 


Asr, — 0, only those terms remain in which r, cntcrs into F,,; the remaining 
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terms all vanish, 1.e. 
2 
Y(r [ro r3)~ gv Alles r2)F 4(r;) 


1 1 
Y 4(r3|ry, 2) SO Fev alrsira)Faslry) — Valles r2)F 42(r;) 


1 1 
P(ra[r3, ry)~ Aral) Fasr) + TeV Ara r3)F 42(r1) 
If we collect these terms and substitute them into (25.26), we obtain 
(when r,; > o) 
«38, —a2 Bs 


3 P(r, 01,1202,1303) ~ Fash alr2|rs)or aa 


2 tas == 
i F gal Walle ha) nna 3 Pats 2 Bas (25.27) 


V6 


The functions 


ei 5 21(Ba%s— Bats (25.28) 
and 
t 
kz = 7g LPP aes — (es B +a B3)] (25.29) 


are the spin functions for the electron + helium atom system which are 
normalized to unity (normalization is ensured by choosing c, and c, in 
accordance with (24.21)). The function y, describes the state of the system 
in which the atom has zero spin, and the z-component of the spin of the 
incident electron along a special direction is +4. The function y, describes 
the superposition of the following two spin states: 

(a) The triplet state of the atom with the spin z-component equal to 1 
and the state of the incident electron with the spin z-component equal 
to —4. 

(b) The triplet state of the atom with the z-component equal to 0 and 
the state of the incident electron with z-component equal to +4. 

Therefore the wave function for the singlet state of the atom constructed 
in accordance with (24.6) is 


J 
P22 02,1363) = Waltalts) 75s B,—%, B3) (25.30) 


The expression 
Pa2(¥2 02,1303) = Wa(\ro, 13) X3 (25.31) 


is the wave function for the triplet state of the atom, constructed in accordance 
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with (24.6) with the spin z-component of the atom equal to I. The expression 


3B y+ 
O'42(¥2 02,1303) = Wal|r3, 2) ——— = a : (25.32) 
/ 2 


is also the function for the triplet state, but with zero spin z-component. 
In view of (25.30)-(25.32) and (25.24), we can rewrite (25.27) in the form 


olka 


J3 50 (7104,1%2 02,1303) ~ Par dar % 


F147 00 ry 
a ik aan | eikazrt 
+QO429a2P1 RE ——— + 429 42% RE (25.33) 
3 Oy 3 Oy 


The current density along the r, axis in configuration space yields the 
cross-section for the excitation of the singlet state 


do, = ——|q43|7 dQ (25.34) 


The excitation cross-section for the triplet state with spin z-component 
equal to | is 


dos) = Tl qa2|" dQ 


Whereas the excitation cross-section ie the triplet state with zero spin 
z-component is 


da§?) = : “a2 


ye, daal” dQ. 
On summing over all the possible spin —_, in the final statc, we find 
that 


Kao 
do, =—— 9.42 
ko 


2 dQ (25.35) 


is the excitation cross-section for the triplet state with arbitrary spin z-com- 
ponent after collision. 


26. ELASTIC SCATTERING OF AN ELECTRON BY A HELIUM ATOM 


When polarization is neglected the elastic scattering of an electron by a 
helium atom is described by 


AF +(k?~2V)F =2 | WF dr’ (26.1) 


where V is defined by (25.9), (25.10) and W by (25.11) or (25.15) through the 
ground state function W(r, |r.) (the subscripts of V, F and HW’ are omitted 
for simplicity). 

The function W(r,]r.) may be represented approximatcly as a product 
of the onc-electron functions 


W(ry|r2) = w(ry)u(r2) (26.2) 
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We then have 

ju(r’)|? dr’ 
jr—r"| 

The kernel of the exchange operator in (25.11) is 


22 2 1 
Wr’, r) = u(r) EN + 4A, + . + oe tae + 4ko+ Eo 


Vea (26.3) 


”" yt J n” Nt / 
+ | ue )( 44" 4 + — vor] — rr) u(r”) dr u(r’) (26.4) 


Using (25.15), we obtain 


2/,,! " 
Wr',r) = u(r) 44, + 4k he — — — | u(r’) (26.5) 
ry |r—r'| [r”—r’| 


The expressions given by (26.4) and (26.5) are, in general, different. However, 
they become identical in one important special case, namely, when u(r) 
satisfies the equation for the self-consistent field 


2 ur 
Au+2|-— | ——\ dr’ Jut2du =0 (26.6) 
r J |r—r'| 


and if instead of the exact energy E of the helium atom, equation (26.4) 
contains the energy calculated by the self-consistent field method, 1.e. 


7 wr’), 
= 24 — | ae dr dr (26.7) 


Equation (26.6) may be used to reduce (26.5) to the even simpler form 


W = u(r)u(r’) (34: + A+ ra) (26.8) 


which 1s analogous to the kernel of the exchange operator in the case of 
scattering of an electron by a hydrogen atom in the antisymmetric case [see 
equation (13.2). 

If the function u(r) does not satisfy the equation for the self-consistent 
field, so that equation (26.4) cannot be reduced to (26.5), the problem 
arises as to which of them should be used in calculations. The variant 
represented by (26.4) was used in the calculations considered below [81, 82]. 
However, it may be argued that (26.5) is preferable. The point is that (26.4) 
and (26.5) differ considerably at small r: in (26.4) the function W(r’, r) 
includes the term 2/r, whereas (26.5) does not include this term. When 
the replacement (r,, 13) > (@’,r) 1s introduced into the exact expression 
(25.15), the function W(r,, r3) does not contain 2/r. Thus equation (26.5) 
is closer to the exact exchange operator than (26.4). The term 2/r in (26.4) 
appears only because the wave function 1s not exact. 

Equations for the partial waves can be obtained by separating the 
angular variables in (26.1). Numerical integration of the equations for the s 
and p waves using the atomic functions for the self-consistent field was 
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carried out as far back as 1933 [83]. The dependence of the phases 6, and 0, 
on the wave number which was obtained in that paper is illustrated in 
Fig. 38, together with the values of these phascs without allowance for 
exchange. Morse and Allis used 69 and 6, to determine the angular dis- 
tribution and the total scattering cross-section; thcy compared these with the 
experimental data existing at that time (see Figs. 39, 40). It is evident from 
these figures that, in general, the theoretical calculations reproduce quite 
satisfactorily the observed dependence of the cross-section on the energy 


I80 
90 
3 
7 2 
do 
3 5, 
| 4 
O 1 2 3 
k 


Fic. 38. Phases for an electron scattered by a helium atom: 1, 3—with exchange; 2, 4— 
without exchange 


and the angular distribution, including the characteristic backward scattering 
maximum at 180° at low energies, which is similar to that observed in the 
case of the hydrogen atom (§ 17). Figure 39 can be used to dctcrmine the 
scattering length. It is found to be a2 1:5 which corresponds to 
o = 4na* ~ 28a. 

The scattering length was recently recalculated by numerical integration 
of the equation for s-waves of zcro energy [82]. The one-electron atomic 
function u(r) was chosen in the form 


u(r) = N(e~*+Ce 27") (26.9) 
where JN is the normalization factor and 
C = 0-6, Z = 1°4558 (26.10) 


This function is a good analytical approximation to the function for the 
sclf-consistent field: the calculation yields a = 1:442, corresponding to 
o = 26:lag. The cross-section at zero cnergy obtained by cxtrapolation of 
the experimental data was found to be o = 22. 

Thus the numerical calculations in [82] lead to a result which cxcccds 
the cxpcrimental valuc by only 18%. Part of this discrepancy is due to the 
fact that the polarization of the hclium atom by the scattcred clectron was 
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6) 10) 20 30 40 
E (eV) 


Fic. 39. Elastic scattering cross-section for an electron scattered by a helium atom (in 
units of 77a 7): curve—theoretical, points—experimental 


not included. Some of the error is due to the approximate nature of the wave 
function for the helium atom, including the fact that the approximate 
functions W(r,|r,) given by (26.2) do not contain the distance r,, between 
the electrons and therefore the correlation between the atomic electrons is 
ignored. 

The variational method was used in connection with the scattering of an 
electron by a helium atom in [81]. The one-electron atomic function was 
taken in the simple form 


Uy = Ne” (26.11) 


O 30 60 90 120 I50 I80 
@ (deg) 


Fic. 40. Angular distribution of electrons scattered by a helium atom at electron 

energies of 6eV (curve 1), 20 eV (curve 2) and 30 eV (curve 3). The broken curves 

represent theoretical calculations without exchange, the solid curves the theoretical 
values with exchange allowed for, and the points the experimental values 
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27 Z3\t 
Z=— N= (=) (26.12) 


The exchange operator was determined with the aid of (26.4), using the 
approximate expression for the energy 


where 


E, = ~2Z? (26.13) 
The trial function was chosen in the form 
f=sinkr+(at+be-“)(1—e7“")coskr (26.14) 


The phase 6,(k) obtained from these calculations is in good agreement 
with the results of numerical integration in spite of the fact that different 
one-electron atomic functions u(r) were used in the two cases. 

The scattering length determined by the variational method was 
a = 1-46, quite close to the result obtained by numerical integration. 

Finally, the phase 6) was determined in [84] by an approximate solution 
of the integral equation for the s-wave. The resulting scattering length was 
found to be somewhat lower than that obtained by the variational method. 


27. EXCITATION OF THE 2!S- AND 2°,S-STATES OF THE HELIUM ATOM 


The excitation of the 2'S- and 2°S-states of the helium atom was cal- 
culated in [85] by the distorted-wave approximation using the variational 
method; in [84] approximate expressions were obtained for the distorted 
waves with the aid of integral equations. 

Here we shall confine our attention to the calculations described in [85]. 
On this approximation, the equations for the wave function of the scattered 
electron are 


AF, +(k2—2Voo)Fo—2 { WooF dr’ =0 (27.1) 


AF, +(k2—2V,,)F,—2 | W,,F,dr’ =2 (V20 Fo ze | Woo Fo dr’) (27.2) 


AF, +(k}—2V,,)F:-2{ W,,F,dr' = 2] Wo Fodr' (27.3) 
with the asymptotic behaviour 
Se i tkor piki.2 
Fo =e" + ae Fij2 = 41,2 (27.4) 


The function Fy describes the scattering of an electron when the target 
atom is in the ground state; F, and F, describe the excitation of the 2°S- 
and 2'S-states respectively. Following the distorted-wave method, we 
neglect ,,F, and W,,F,, which connect the excitation of the triplet and 
singlet states, and the effect of the functions F, and F, on Fo. As can be 
seen from (27.3), excitation of the triplet state is only possible as a result 
of the exchange interaction IV,9. The quantities V,, and J¥’,, are determined 
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from (25.9) and (25.11). The wave function for the ground state is taken 
in the form 


Wo(ry|r2) = Ud(Ty)Uo(T2) (27.5) 
where wy is given by (26.11). The wave function for the triplet 2°.S-state 
is taken in the form 


Wilh, ro) = Ny[uy(ry)ua(r2)— uo(r,)us(r2)] (27.6) 
with u, and w, given by [86] 
u(r) =e 7 
u,(r)=re "’—ce-™ Ci) 
where 
C= 2. fp =0-61, y = 1:57 
and 
_ 3(a+y)° 
— (B+ayt 


The coefficient c is chosen so that the function uw,(r) is orthogonal to 
u,(vr), and the parameters f, y determined from the condition for minimum 
energy. The function for the singlet 2’\S-state is taken as 


Wo(r:|r2) = N[uy(ry)u2(r2) +u2(ry)uy(r2)] (27.8) 
where u,(r) and w,(r) are in the same form as for the 2°S-state [equation 
(27.7)] in order to simplify the calculations. The parameters of the one- 
electron functions u, and uw, are different for the triplet and singlet states 
when determined from the minimum energy condition. 

The function w,(r,|r,) given by (27.8) is not orthogonal to Wo(r, |r2). 
As a result V,, includes terms which fall off as r~* when r > oo. In reality, 
V4 should fall off exponentially when r > oo. For this reason, terms which 
behave as r~! when r > oo are rejected in V,,, and this removes the most 
undesirable consequence of the fact that Wo(r, |r.) and w,(r,|r,) are not 
orthogonal. There remains the error in the numerical value of the terms 
having the correct asymptotic behaviour, which is due to non-orthogonality. 
This error was not estimated in [85]. Only the s-waves were taken into 
account in the expansion of F, into the partial waves. Therefore the problem 
reduces to the determination of the three functions 


o(r) = rF,(r) 
f1,24") = rF (1) 


where F ,(r) and #.,(r) are the solutions of (27.2) and (27.3) without the 
right-hand side terms; /, and /, have the asymptotic behaviour 


(27.9) 


fa ™ : e= sin (k,r+Nq) (27.10) 


The amplitudes qg, and q, are expressed in terms of fo, “,; ,and #, by formulae 
analogous to (18.27). 
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The latter three functions are determined by a variational technique. 
The trial funetion fp 1s chosen to be analogous to (26.14). When the 
asymptotic behaviour given by (27.10) is taken into account, the trial function 
becomes 

ello 


fo= 


where Z = 27/16. 
The trial functions for /, and /, are taken to be 


(L+a?)~?[sinkgr+(at+be *)(l—e-“)coskgr] (27.11) 


ini 
Pie = (1423 “(1+ Ce)sinkyr+4,(l—e)eosk,r] (27.12) 
1 


and similarly for “,. Calculations show that (27.12) is a very poor approxi- 
mation for “, at low energies near k, = 0, and therefore n, and q, are 
calculated for A, 2 0-3. The corresponding effective cross-section is shown 
in Fig. 41. For 4, and g, the calculation is carried out with &, > 0-1, and 


4 


o (107* 1108) 
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Fic. 41. Partial cross-sections for the excitation of the 24.8 (curve 1) and 2°S (curve 2) 
states of helium at / = 0. The energy is measured from the threshold 


a determination is made of the scattering length for the 2° S-state of the atom 
(although the numerical result is not given in [85]). 

In order to illustrate the importance of exchange, the parameters n, 
and q, have been calculated without introducing this effect. Figure 42 shows 
the dependence of 42, y, and 4,2, on & without exchange. The broken 
eurve shows the part determined on the basis of the estimated scattering 
length, i.e. the slope of the (A) curve for k = 0. As can be scen from the 
figure, even with relatively small changes corresponding to k 2 0:5, the 
ordcr of magnitude of the phase and its dcpendence on k is very nearly the 
same with and without exchange. However, near the threshold (k ~ 0-1-0:2), 
the inclusion of exchange has a radical effeet on the variation of the phases: 
the scattering length in the 2°S-state increases very rapidly as a result of 
exchange. The phase 7, approaches 32/2 when k = 0-05, 1.c. very near 
the threshold. In other words, seattering by a hclium atom in the 2?S-state 
is very close to resonancc scattering. This is associated with the existence 
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of a negative metastable He™ state which lies about 0-1 eV above the 
2°S-state of He. 

The calculated excitation cross-section in the 2°S-state near the threshold 
is shown in Fig. 43, which also gives the experimental values [87]. As can 
be seen, the calculations show a much sharper peak near the threshold than Is 
observed in experiment. This is not surprising when it is recalled that 


Fic. 42. The effect of exchange on the phases of trial functions: 1—phase 7,2 with 
exchange, 2—y72 and 3—”, 
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Fic. 43. The variation of the cross-section for the excitation of the 2°S-state of helium 
near the threshold: 1—calculated [75], 2—experimental [77] 


scattering by a helium atom in the 2°S-state is very near to resonance 
scattering. Under such conditions, the phase 7, and therefore the calculated 
cross-section for the 2°S-state, are very sensitive to small changes in the form 
of the wave function inside the atom, so that the variational method using 
simple trial functions is unsatisfactory. No exact numerical calculations 
have as yet been carried out. 


E.A.C. 10 
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A possible interpretation of the behaviour of the 2°S-state excitation 
cross-section near the threshold is given in [88], where it is suggested that 
inelastic electron scattering giving rise to excitation of the 2°S-state is 
analogous to nuclear reactions which proceed through the formation of a 
compound system (in the present case He” in the 2S, state). If this is so, 
the effective cross-section for s-scattering near the threshold should be 
described by the Breit-Wigner formula 


7 
OSS SS Se ee SS eS 
ko (E-E,)? +4(Po +P)’ 
where E, is the energy of the compound system, I, is the width for the decay 
of the ion into an atom in the ‘S-state and an electron, and I, is the width 


for the decay of the atom into an atom in the 2°S-state and an electron. 
The widths [y and I, may be written in the form 


we = 2ko Vos I, 2k, V7 (27.14) 
IT, and I, are the so-called reduced widths and are independent of energy. 


(27.13) 
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Fic. 44. Cross-section for the excitation of the 2°S-state of helium on the compound 
system model: curve—Breit-Wigner formula with the parameters given by (27.15), 
points—experimental values [77]. Excitation threshold E,, = 19°79 eV 


Figure 44 shows a plot of o(£) as given by (27.13), with the following 
empirical parameters 


E, = 20:19 eV 
yg = 2:2.10°? eV.cm (27.15) 
yt = 2:5.107'2eV.cm 
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These values of the parameters are not the only ones to give adequate 
agreement with experiment. Another selection of the parameters is also 
discussed in [88]. 

For comparison, Fig. 44 also shows the experimental points, and as can 
be seen the agreement is quite good. The magnitude of yp , corresponds 
to a width T ~ | eV or a lifetime ~ 107~!° sec, which is of the same order 
as the half-life for the auto-ionization of He in the 2s*'S) state (about 
2:5x 107! sec). The experimental o(£) curve exhibits one further feature. 
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Fic. 45. Cross-sectitons for the excitation of the 2°S-state of helium: 1—experimental, 
2—calculated for /=Oand 1. E is the energy measured from the threshold 


It is found that near the threshold the values of o lie very nearly on a straight 
line and therefore in this region o may be regarded as a linear function of 
the energy measured from the threshold: 


og ~ (E-E,) 
However, it follows from the general expression o = (k,/ky) J |q| 7dQ that 


if g remains constant for k, = 0, and this is certainly the case for neutral 
atoms, then o should be proportional to k,, 1.e. 


oO™~ (Ea 
It appears that the region near threshold where g is approximately indepen- 


dent of energy and o ~ ./E—E, is very narrow and is not detected owing 
to the low resolving power of the instruments. 
Direct comparison of the calculations of the partial cross-section for 
= (0 with the experimental data at energies of a few eV above the threshold 
has no meaning because the partial wave with / = 1 introduces a considerable 
contribution at such energies. The contribution due to the p-wave was 
estimated in [85] with V,, and }V,, neglected. The results are given in Fig. 45 
together with experimental data. 


28. THE 2!S—2°,S TRANSITION DURING A COLLISION WITH AN ELECTRON 


The transition from the metastable 2'S-state to the lower 2°S-state in 
which the excess energy is given to the incident electron was discussed in 
10—2 
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[89], in which the system of equations for the s-waves was numerically 
integrated on the strong coupling approximation. The equations were 


ae r r 
2+ (2 adhi-2 | Warhrdr’ =2 | Wash dr’ (28.1) 
0 0 
aay 2 va , r r 
Git wh-2 | # 11i14r = 2 o/s i2f2 dr (28.2) 
0 0 
with the asymptotic behaviour 
ae ; 
fi ame pe kar+4@e ear ti cm qy eikir (28.3) 


2 


where f, describes the incident electron with energy 343 and /f, describes 


the electron after collision, with energy 34] = 4k3+AE. The partial cross- 
section for the process 1s 


= 4nla,)?— (28.4) 
Z 


Y 4, and #,, were calculated from the general formulae (25.9) and (25.11) 
using the analytical expression (27.7) for the atomic functions. Table 11 
gives the results of the calculations, with the cross-sections o in units of ao. 


TABLE 11 
ky E (eV) Co 
0-044 0-026 180 
0-050 0-034 14] 
0-060 0-048 113 
0-100 0-134 64-7 
0-300 1-210 14-6 


The calculated cross-section at 0-026 eV is smaller than that observed 
by a factor of about six. This discrepancy may be due in part to the 
fact that the functions for the 2'S- and 1'S-states are not orthogonal. 
However, this would hardly lead to an error of this magnitude. Another 
possible explanation of the discrepancy is the formation of the metastable 
He™ ion. However, no quantitative estimates of these effects are available. 


29. EXCITATION OF SINGLET 1'P-STATES (11 = 2, 3, 4) 


Excitation of the singlet P-states was analysed in [90] on the strong- 
coupling approximation, but without allowance for exchange effects. 

In order to elucidate the relative importance of the various partial waves, 
a calculation was made of the cross-sections corresponding to total angular 
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momenta L = 0, 1, 2and 3. In the group of indices v’ = (n’, /;, /3) referring 
to the initial state, /; = 0 for all the transitions under consideration while 
> = L. In the group of indices v = (,/,,/,) referring to the final state, 
/, = 1, while /, can assume the values |L—1|</, < L+1. When in 
accordance with the laws of composition of angular momenta there were 
more than two equations, they were split into independent pairs. Tables 
12-14 give the calculated values of the effective excitation cross-sections for 
the levels 11S—3'P (E,, = 23-308 eV), 1'S+2!P (Ey, = 21-211 eV), and 
11S = 4'P (Ey, = 23-734 eV) in units of 2a/100. As can be seen from these 


TABLE 12 
ky E (eV) C10 Oo. O21 O12 O32 023 o 
0-05 23:342 0:°446 0-0°20+ 0-172 00°15 0:0446 — 0-649 


0:10 23-444 0-245 0-0743 Jez 0-:0°35  0:0°56 0-0°37 1-808 
0:20 23-852 0-166 0-210 S22 00712 0-606 0-0412 4-503 
0-40 25:484 0-132 0-430 Dede 0-0726 0-745 0:0712 4-031 
0-60 27:204 0-102 0-366 2°12 0-0732  =1-135 0-0747 3°731 
0:80 32:012 0-074 0-298 P52 0-:0728 1-286 0:07978 3-191 
1-0 36:908  0:0503 0-233 1:02 0-:0717 = 1-220 0:0139 2°532 


TABLE 13 
ky E (eV) C10 Oo. O21 Oj2 
0-05 21:245 1-13 0:0°43 0-036 — 
0:10 21°347 0:976 0:0713 0:29] 0:0&24 
0-20 21:755 0-601 0:024 3°72 0:0°68 
0-60 25°107 0-274 1-28 9-16 0:0743 
1-00 34°811 0-128 0-810 4-44 0:07420 

TABLE 14 
ky E (eV) C10 Co1 O21 O12 
0:05 23:°768 0-139 0:0°25 0-250 0:0°61 
0-10 23°870 0-092 0-011 1-075 0:0710 
0:20 24-278 0-075 0-541 1-210 0:0723 
0-60 27-630 0-050 0-172 0-820 0:0727 
1-00 37°334 0-026 — 0-406 0:0717 


t Superscripts indicate number of zeros after the decimal point, for example 
0:0°43 = 0-00043. 
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tables, at the energy of only | eV above the threshold, the main contri- 
bution to the excitation cross-section is due to the partial cross-section g;. 
In some regions the latter is much greater than all the remaining cross- 
sections. The origin of this effect was discussed earlier in § 18. 


TABLE 15 
O01 T10 
V 
Ms ew) S.-C. d.-w. S.-C. d.-w. 
0-1 23-444 00°43 «0-044 0-245 0-233 
0:6 27°204 0-366 0-346 0-102 0-096 
1:0 36:908 0-233 0-186 0-050 0-040 


The relative importance of strong coupling may be estimated from 
Table 15 which gives the cross-sections for the excitation of the 3'P-level 
(in units of 2a3/100) calculated by the distorted-wave (d.-w.) and the strong- 
coupling (s.-c.) approximations. 


Chapter 8 


GENERAL EQUATIONS DESCRIBING THE COLLISION OF AN 
ELECTRON WITH AN ARBITRARY ATOM OR ION 


30. THE WAVE FUNCTION FOR THE ATOM + ELECTRON SYSTEM 


CONSIDER an (7+1)-electron system consisting of an n-electron atom (ion) 
and an incident electron. According to § 24, the coordinate part of the 


wave function P(r,, ..., rr. [ Peet. +++> fn 1) IS antisymmetric, i.e. 
(+P, )¥ =0 (30.1) 
fi,j < k or i,j > k; it is also cyclically symmetric, i.e. 
n+l 
(1 —- P..] =O (30.2) 
i=k+1 


The number of terms on the left of the vertical line is related to the spin s 
(which is assumed to be fixed): 


k= a = (30.3) 


By analogy with the treatment of the helium atom, we shall construct 

the function 

Pisces Ppineou hia) 
from the products of n-electron atomic functions w and one-electron functions 
F, so that (30.1) and (30.2) are satisfied. 

For a given square of the spin angular momentum of the (7+ 1)-electron 
system, there are two possible spin states of the sub-system of n-electrons, 
namely 

Ss, =s—4, Sy =st4 
The corresponding functions are 
wid, vedi [Fi $1, nh) and w,(, Leeysda| dat Io e2eq dt) 
According to equation (30.3), 
Jp ein-G-HsHk  j,sin-Gt+h)ak-i (30.4) 

We shall therefore consider the products of the n-electron function wW 
and one-electron functions F of the form 

Widl,...,klk+1,..., mF (n+l) 
and 
Wil, ...,kK—-Ik, om F{(n+1) 
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Consider the sums 


n+1 
y= y A; Ping W1,...,k[K+1,..., mF y(n +1) (30.5) 
isk+1 
and 
k nt+1 
Y,=) 2b Pie Pine Holl, Kk, n)F(nt+1) (30.6) 


t=1 i= 
with the coefficients a; and 6, ; chosen so that (30. 1) and (30.2) are satisfied. 
It is understood that application of the exchange operator results in the 
following transformations 


Pp tl 2 ace Gaveis Oy ccd) Sf ls Zara Dewre Opctee)) 
and 


Pe Ils cide he wa Dy eee IG Ms easing Ds atay))) 
Nines D yratos yng GO laxctns Di cats Oyuetag Ml) 
i.e. the a-th and b-th arguments can be exchanged independently of their 
position within the brackets. 
We note that all the steps necessary to obtain a; and 5, ; may be performed 
in compact form by using the commutation relations 
PoP ee = ~~ Pel = Piet ac (30.7) 
These relations can easily be verified. For example, 
Pap Py-flabe) = Pa, flach) = f(bea) 
Pa Pap flabe) = Pa-f(bac) = f(bea) 
The coefficients a; and 5; , are defined unambiguously to within a common 
factor which is unimportant and can be included in the normalization. 
We shall not go into details but will merely write down the final result and 
then verify (30.1) and (30.2) [91]. The expressions for ‘¥, and ‘¥, are of the 
form 


=<; (1 ae), Pras) Vall vay k{k+ l,..., mF ,(n4+]) (30.8) 


k-1 n 
Wa =ea(1- Pus }[t- by Pane #124 2)P hat 
t= 1 i=k+1 
xWo(1, 4 K— 1k, mF2(n+1) (30.9) 


The coefficients c, and c, can be defined from the normalization conditions. 
Suppose that yw, and w, belong to the discrcte spectrum while F, and F, 
belong to the continuous spectrum, and that YW, , arc normalized to unity, 
but that F, and F, are normalized as a planc wave to (2n)°?d(kK+4’). If 
we require that Y, and , should also be normalized to (27)*d(A +4’), 
then 


cy =(n—k+1)73 


30.10 
c, =k 4[(n—2k+2)*+(n-—k41)]73 ( ) 


8. COLLISION OF AN ELECTRON WITH AN ARBITRARY ATOM OR ION 141 


We shall now show that Y, satisfies (30.1) and (30.2) (the normalization 
factors will be omitted from the proof). 

1. The anti-symmetry condition is clearly satisfied up to the vertical 
line, since the arguments of , and w, up to the line are identical. 

2. In order to establish the anti-symmetry condition after the vertical 
linc let us apply the operator (1+P;, ;) and use equation (30.7). The result 


(1 - Y Pao U4 PW Pa lj#Antli 


— 2 Pratl + Pid Fi, j=nt+l 
#i 
which 1s zero because 

(L+Pi dW, =0 


3. The cyclic symmetry may be established as follows. From (30.7) it 
can be shown that 


n+ n n n 
(1 $F! n,)(1— $ Pras) =(1- SP ies)(1- fra) 
jJrukt+il i=k+1 jrzk jJ=kt+1 


it follows that 


The proof for ‘¥, is quite similar. 
1. The anti-symmetry condition up to the vertical line can be proved by 
operating with (1+P;, ;) on ‘V,: 


k-1 - 
( = Pas) [ 2 Pine +(0=2K+2)P rae 
t= juktl 
XCL+Pi,)2F2 i j#ntl 
(1+P, )'¥,= a 
(1 7 s P| |(1 = 2k +2)Pe ne (E+ Pi, 
| t=1 | 
| —YPuneiltPid]vaFa, fant 


which is zero because 
(+P, )W2=A+P, JW. = 9 
2. When i, 7 # A, the anti-symmetry condition after the vertical line is 
obviously satisfied. When j = k we have 


(1 + P; ,)\P2 — = 2 Pax! a y Eipaaa UP BIE AE xa (+P, )W2F2 
t# i 


Which is equal to zero because 
(+P, )W2 = 0 
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3. The cyclic symmetry condition may be verified by using the properties 
of w,. After some transformations we have 


n+1 n n 
(Si ru)e (Hof roll 8) 
j=Hkt+i1 j=k j=kti 


-¥ P(t 3 P,.j) = (= 2k+2) WF, 


Next, it may be shown es 
(1 Sms) EPa(t- SPs) ve = (1—2k+ 2, 
j=kti anit 
and therefore 
n+l 
(1 — \ Pas] Y, — 0 
J=k+1 
From the general formulae (30.8) and (30.9) we have, as a special case 

of the expression used in § 24, 


l 
¥ (1[23) = a [w(1|2)F(3)— w1|3)F(2)] 


| 
‘Y2(1|23) = er [y(|12)F (3) —W(]L3)F(2) + 2W(]32)F(1)] 
‘N 


We can use the above results to write down an expression for the coordi- 
nate function of the atom + electron system in terms of the atomic and 
one-electron functions. Let A be the set of numbers indicating the energy 
of the atom, its orbital angular momentum and the z-component of the 
angular momentum. We shall write Y as the sum 


=F (Pat ha) (30.11) 


Each of the terms '¥4, 2 willconsist of expressions of the form of (30.8) or (30.9). 
If we omit the common factor of Wy, 5, which will be unimportant in 
the ensuing anaes we have 


Y, = (1- y Putt) Yak . 


i=k+1 

(n—2k+2)2 \- ( ka 
re ae k l—) P 
ore ( n—k+l pa os 


x f - ¥ Pane tO 2+ DP ins W>Fy> (30.12) 


=k+1 


31. EQUATIONS FOR THE FUNCTIONS F,, 


In order to establish the equations for Fy,, we shall use a relation analo- 
gous to (10.24), Le. 


| Wi(H-S)¥ dry, ...,dry = 0, nn ee, (31.1) 


in which #/ is the energy operator for the electron + atom system. 
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We shall suppose that W4, is the exact atomic function so that 


(Ho — Eg Aa = 0 (31.2) 
The exact functions will be replaced by the approximate atomic functions 
only in the final expressions which will be of the density-matrix type. 
There is another way of deriving the equations for F',, [92, 93], in which 
it is supposed from the outset that w,, are the functions for the self-consistent 
field. 
Let us now write the energy operator for the atom + electron system 
in the form 
H = Ho—43A,4,+V(14,---5 tne) (31.3) 
where V is the energy of interaction of the (7+1)th electron with the 
nucleus and the remaining n electrons; it is given by 
c | 
+ hee aba a es 
Pn+i 2 rneimti 
Substituting the expression for ‘Y and H given by (30.11), (30.12) and (31.3), 
(31.1) and using the orthogonality and symmetry of w, we have, in view 
of (31.2) 


V=- 


(31.4) 


(A+k4a)F 40 = 2 ys VaA'a Fata Zi \ Wigeaz F arg’ (31.5) 
A’ A'‘a’ 
where 
Vaardltne i) =| Weare dis ves dr (31.6) 
is the interaction energy matrix, 
ka = 2(6— E4,) (31.7) 


and W is the exchange operator of the form 


Waal ier ie4) = Poa | Wall kas r, (Ay +kiq—2V) 


MW gai Te cys cheese PP nel har onwdl, 
+ aa | Waals ---|-+-5 tn) (Ant kag —2V) 
X Wye (11, alone rie DF acg(tn) aN 15 0005 UY, (31.8) 
Pao’ 2Nd gq are given by 
Pit = Pr2 = Par = 9 ; (31.9) 
Po2 = (K—1)(n—2k +3)(n—2k42)7! 
Qqy=n—k 
di. = kA(n—k)(n—2k43)(n—k +1)? f(n—2k 4-2)? +n-—k +1473 
n—k+1 73 
Jar = (nk) tea | 


Go. = —(n—k+1)(n—2k4+2)7' 
Using the explicit expression for ]”, the symmetry properties of w,, and the 
fact that A is self-adjoint, the potential energy matrix and the exchange 
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Operator may be reduced to the forms 


Z  f Paar) d | 
Vastgu|lo-n a 
AA‘a : +| ry (31.10) 
Wasea Fyre = [ Waareal" OE va") de (31.11) 


? / z J , 
Wearaalt’,) = (34 ei a ) Paseo r 
r lr—r | 
-| Paataa("'’s r’, r) 


=r 


dr” (31.12) 


where the p represent expressions of the density-matrix type. 

In order to write these expressions in the most compact form, we shall 
introduce a special symbol for the integral of the product W4,W,4’, over all 
the variables with the exception of a specified few. For example, the integral 
with respect to all the variables r,,...,7r, except for r; will be denoted by 
the symbol 


(rad da‘a’ =| Wal! ara’ AV 1, 060, AN AN 4 yy 0s UN, (31.13) 


If, after integration, r; in w,, is replaced by, say, r, while in W,-,: it is replaced 
by r’, this will be denoted by 
r r’ 
Pil} Aa \Ti/] A‘a’ 
In this notation we have 
; r’ r’ r’ 
meteree(").(0),,40-(7),(9) 
nJAt\F1/ At Ma7AL\Tn/ Al 
(31.14) 
r r’ 
pass) = (k= 1) (7 ,) (a +(n-k41)(7 ) 4 (31.15) 
n/ A’2 rn A2 n/ A’2 


paar) = air( ) | (31.16) 


r 
rn 


ry (: \ (31.17) 
) ( )., (31.18) 
r’ r 
r 2 x iat? eee we an rate 
? r 
; ( an Py ae 
+n k=) (F yO i (31.20) 


Paan2's") = 442 


r 


; r 
Paar s") = G22 


Paani ors) = 414 


Paars(hs!) = dry . 
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y! yr” r 
Paani yr = 412 |- D(7 dale ie 


iL f 4 7 f "7 
r pr r sr r” r r r 
ry VyJ ar \lR Va) a2 Pa-1 UnJ At \on-1 = 'n/ A’2 


Paar er = 421 [ke v(f. ules As 


ee aa ae r r’ r” r 
VE alas ote soko (" 7 a 
Pe Vat a2\lR 1 Pa-1 UTnJa2z2\Tn-1 "n/a 


" yr’ r” r 
ror’ rn) = Kl 
Paar2t"; )= 422 ( 1 Dy rae 
+= (7 4 (" ‘ | 
Pa-1 Unf A2\Tn-1 n/ a'2 
rr" rr 

+ k-—2 

P22 [ (7 ani an 


yr’ r” r r” 
+n ke)(7 ae 3 on (31.23) 


Using the expression for the exchange operator we can reduce (31.5) to 
(Q+KR Fae = 25 VaweFeat2D, | Wasae(t DFaule de’ 31.24) 
A’ A’a’ 


The asymptotic expression for F,, in the case of a collision with a 
neutral atom Is 


eikaar 
each 
Py 4 e' < "04g Aone + GAotgAa r (31.25) 
while in the case of a collision with an ion it 1s 
26 75 
ik aat + 2 In(kaat — kar) lL i€kaar+—e In 2k Aart 
Fag SOE | ‘ nae . 7 q ApagAa ~ i . aa ‘ ) (31.26) 
r 


where k,, = K,,V. 

The system of integro-differential equations (31.24) subject to the 
asymptotic conditions (31.25) or (31.26) may be transformed to a system of 
integral equations. Let G,, be the Green’s function which in the case of 
a collision with a neutral atom 1s of the form 

1 eikaalr—r’| 


G, —_—_— 


a 


while in the case of a collision with an ion of charge Z it satisfies the equation 


Z 
en Care 4 G,,=—-i(r—-r’) (31.28) 
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By analogy with (10.53), consider the kernel 


Z 
K 4 4'aa’ eos Gadr, r ‘) (Vv AA’a” bax) Oaa’ -| Galt, r’) WW, AA’‘aa’ (r’, Pr ‘)dr" 
r 
(31.29) 
(in the case of collisions with neutral atoms Z = 0). We then have 
Fee luta S | Karadts?)Eyelde (31.30) 
A’a’ 


where g,, represents either exp (ik ,,r) in the case of collision with a neutral 
atom, or else the Coulomb function with the asymptotic behaviour 


Z 
5 re ee i(kaar+ ——In 2kaar 
D4, sshd [A A0 a (kK Aa Kaa | + ee (31.31) 
It is known [57] that this function can also be written in the form 
ue DY eye Me 2 
o=e*T ale e'*" F “Liar ke) (31.32) 


where Fis the degenerate hypergeometric function. 

From equation (31.30) it is possible to obtain the expression for the 
amplitude g,4...4,. One must consider the asymptotic behaviour at large r. 
The Green’s function G,,(r, r’) for large r and finite r° has the following 
asymptotic properties: for collisions with neutral atoms 


ee ee se (31.33) 
whereas for collisions with ions 
1 i( aor + In 2k Aer ) . 
G eh e Kaa Y Aa (31.34) 
4n 
where 
nZ . 
= Z\ Z 
Dp = e2k r(1 +it) elk aar F | - a le —= (kr + kr) (31.35) 


Substituting (31.33) or (31.34) into (31.30) and bearing in mind (31.25) 
and — we have for collisions with neutral atoms: 


q AgaoAa = —? yf ee E aa’ Vaar a F 4: a Ar") 
rs dt" Wapaal?’, rFaal) dr’ (31.36) 


and for collisions with ions: 


l 7 ; 
Y AgagAa = es » | Pal ) Maes Oaa’ |r ) 
mit A’ Q’ 


+ Waratah 0’) F gah”) ir" dv + qe, 4a,544o Saag (31-37) 


e 
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Consider the asymptotic form of the complete wave function, including 
the spin variables. By analogy with § 25, we shall scparate out the terms 
which are finite when, say, r,4,— 00 (we shall thercfore be conccrned with 
that part of (30.11) which contains the atomic functions for the diseretc 
spectrum). After some transformations 


F ( ) Ree ae | ( ) 

@ ~ ) PD r, oa Fare he Aa a eee ) D> ve a, 

= Al* Al +1 +1 Jin ak 2)? ‘i k 1G F 42 + | +1 
n—2k+2 


a ) . 
ioe Ey ie » Dao gtas DBat1 (31.38) 
where a numerical factor has been omitted. In these expressions, M,, 1S 
the atomic wave function with spin s — $ made up in accordance with 
(24.6) from W,, and x = fy, ..-, Bus pags «+> Xs Py. 1S the atomic wave 
function with spin s+ 4 made up in accordance with (24.6) from W,, and 
X = Bis... Be-vs Les --+5 3 and @',, is made up in accordance with (24.6), 
also from W,., but with y given by 


X —— nk 1y4(1 + 3 Pay) Bu mengihl ps Oras sony O, 
jakt+ 


(Di. corresponds to spin s+} but the spin z-componcent for this function 
is s—4 and not s+}). 

As an cxample, consider the collision of an electron with an atom whosc 
spin prior to the collision is zero. In view of (30.3) we have from (31.36) 


ikaw 
e' 1” € 
D~ O,,1e mage +2, Pat Qaorar + [PY Ostia ns 


ae - es 


elk ar 
an yy D 42 VAgta2—2 F B (31.39) 


k+5 4 


where we have omitted the subscript +1. Henee, when there is no change 
in the spin of the atom, the differential cross-section 1s 


Kay 
K aol 
while if the spin change is allowed, we have on summing over the final 
spin states 


lao ar|? dQ (31.40) 


UO Ast AI = 


k 

I0 Agi a2 = 7 |Aaora2|” dQ (31.41) 

Agl 
If the atom has a non-zero spin, say, So, prior to collision, there are two 
possible valucs for the spin of the system, namely, s,; = so+4 and 
Sy = 5 9—3. This leads to two different problems. In one of them the 
initial state of the atom is described by a function My, , while in the other 
it is described by ®,,. In order to obtain the observed cross-section, the 
result must be averaged over the spins. Considcrable simplification 1s achieved 
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by considering thc valency electrons of the atom and the incident electron 
as moving in the field of the nucleus, and fixing all the internal electrons. 
In this approximation we shall usc the symmetrized function for the system 
consisting of the valency electrons of the atom plus the incident electron. 
The simplest example of this kind is the monovalent atom. According to 
V. A. Fock, the wave function for the valency electron satisfies the equation 


[-4A4+ U(r) ]W.+ Wy, = Eq (31.42) 


where U is the potential energy of interaction of the valency electrons with the 
nucleus and the inner electrons and . given by 


Z 2 
Use |e (31.43) 
r jr—r’ ie 
while Wis the exchange operator defined by 
Ww, = ee ors a rd) a, (31.44) 
Equation (31.32) may be written in the shorter form 
Hr, = EW (31.45) 


The coordinate part of the wave function for the system consisting of 
the valency electron plus the incident electron will satisfy the equation 


Hr) +H,(r’) + —| Y = EW (31.46) 


pr" 


By analogy with (10.46) we shall write in the form 
P=) LOWER War) (31.47) 


From 
[ vo |te HD+ Hale) + a) dr’ =0 


we have the following equation for F, 


(He af *) F (vr) + 2, Ve vy) + 2: | Hr Wee, YF’) dr’ =0 (31.48) 
where 


V, 


ax 


= [ He, aca 


ie 


Way = (By 3 , VO) (31.50) 


The opcrator W’,,- describes the exchange interaction between the 
incident elcctron and the valency electron. The cxchange interaction with 
the inner electrons is described by the operator Win //,. 
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If we let F have the asymptotic form 
gikar 


Pio deal es ae (31.51) 


then by analogy with the case of the hydrogen atom, the scattering cross- 
section assumes the form 

da* = |qi,|* dQ (31.52) 

Separation of the angular variables in (31.24) or (31.48) is carried out 

exactly as in the case of a collision of an electron with a hydrogen atom 


($ 11). 


The equations for the partial waves are of the form 


dq? sts 1) 
pad Act + ki. = ~~ | Saal =2 d V Aa’ ‘all’ Ta’ ‘a’l’ 
r r? 
+2 Dy W garaawl s")faav(r adr’ (31.53) 
0 
where ¥ 4 4a and W 4 goat are given in terms Of Vigyrag ANd Wy grag by 
V aaa = | YOO Waaralt’)Yry(0'9') dO (31.54) 


W aaraa'u(r sr) = rr ve (Op) Wa araa(t FY Yr y(O'') dQ’ dQ" (31.55) 


Owing to the presence of the term |r’ +r"|~* in W,4,9(r', r”), the expression 
given by (31.55) will be different for r <r’ andr > vr’. The two cases will 
be represented by W“ and W?, respectively. The exchange term on the 
right-hand side of (31.53) may be written in the form 


j Woy Sydr' = | (WY—W2) fy dr’ + f W2f.dr' (31.56) 


where vy represents the group of subscripts Aa/. Equation (31.53) cannot 
readily be solved numerically because of the presence of the integral 


: 2 , . : : 
| we) J, dr’, and successive iterations are necessary. However, if we use 


Q 
approximate atomic functions which are expressed in terms of the one- 


electron functions, we can write 
WY? = d, Ry (T)Qiw(’) (31.57) 


This may be used as in § 14 in order to transform equation (31.53) into a 
form which can be solved without iteration. 


Chapter 9 
COLLISION OF AN ELECTRON WITH A SODIUM ATOM 


32. ELASTIC SCATTERING OF ELECTRONS BY SODIUM ATOMS 


As an example of the application of the general theory developed in the 
preceding chapter, let us consider one of the simplest problems, the collision 
of an electron with a sodium atom. In accordance with the discussion in 
§ 31 this problem may be looked upon as involving the motion of two 
electrons, i.e. the incident and the valency electrons in the field of the nucleus 
and of the inner electrons. The latter comprise ten electrons in the ls, 2s 
and 2p shells. 

Among existing calculations of the inelastic scattering of electrons by 
sodium atoms, the most accurate is that reported in [94]. This work is 
concerned with the scattering length, making allowance for exchange and 
3s—3p coupling. It turns out that this coupling plays a very important 
role. Using experimental data on oscillator strengths, the authors of this 
paper conclude that 99-4°% of the polarizability of the atom is due to the 
3s—3p coupling. Therefore, when the single virtual transition 3s—3p is 
included, this accounts almost completely for the polarization of the atom 
by the incident electron. 

The equations discussed in [94] are, in fact, identical with (31.48), 
although they are derived in a different way. The corresponding equations 
for the partial waves on the total angular momentum representation are 
of the form 


(HEP —1kKe fy + > Wwe Svr a 2» | Woyt dvr dr’ — 0 (32.1) 
0 
ta? ve 1) 
aa were ta + U(r) + Ws (32.2) 


and y= (a, th, fy) 


The function U(r) and the exchange operator }f’f can be expressed in 
terms of the densities p(r, r) and p(r, +’) constructed from wave functions 
for the closed shells using (31.33) and (31.34). In the present case 


prs) = Ras") + RoR) 


a y Rap Rap) Vin, 9) ¥11(0 9) (32.3) 


pa 
where R,,, R., and K,, are the radial wave functions for the atomic electrons 


in the closed shells. 
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Using (32.3), we have the following expressions for U and W 


11 
U=—-—+4+2 [ (RE, + R34 3RE Vor r')r’? dr’ (32.4) 
r a 


0 


cO 


lea Te | yrs r)r’'Ry(r’)F(r’) dr’ Ry,(r) 
0 


ee) 


; | y(t re Rod ’)F(r') dr’ Ra,dr) 
0 
31, ! t , t 
bg | tel r Ral FO) dr Rag) 
2 


0 


(3/,+3) 


(21,43) Meir rR (r)F(r') dr’ R2,(r) p (32.5) 
0 


where 


CS 


The quantities ¥,,,, and %,,-, which enter into (32.1) can be expressed 
in terms of the wave functions for the valency electron in the ground 3s 
state or the excited states. 

As was pointed out above, calculations of the scattering length can, 
in practice, be confined to the single 3p excited state. 

The wave functions for the electrons in the closed shells R,,, R,, and R3, 
and the function for the valency electron in the 3s and 3p states were taken 
from [95, 96]. 

The system of equations given by (32.1) can, on the present approxi- 
mation, be reduced to two equations for f399,9 and /311,9 (the subscript A 
which characterizes the state of the atom is now replaced by the triad of 
quantum numbers n, /, 4). In order to find the scattering length we must 


set k3g = 0. We then have k?, = —x? = —0-137. The asymptotic behaviour 
Of f3o0,0 and f3144.9 aS CO is 
fso0,0(1) ~ r—a*; fsitorwe (32.6) 


In [91] the equations for fg99,9 and /3;;,.9 were solved by the method 
of successive approximations. At first the exchange operator W was neglected 
11—2 
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in the operator H, so that H, became the energy operator 7, on the central 


fo @) 


field approximation. The cxchange term | Whar’ was also neglected. 
0 


The zero-order approximation, f°, to the functions f was then substituted 
into the exchange term. 


The solution of the equations with the exchange term | Wf dr gave 


0 
the first approximation f“ which was again substituted for fin the exchange 
term, and so on. The values of a* and a™ given in Table 16 were obtained 
after five successive approximations (# represents the order of the 
approximation). 


TABLE 16 
nN 0 l 2 3 4 
at —2°4 10-0 14-0 13-9 — 
aé —2°4 —9-6 — 5-3 9-5 --6°6 


Next, the effect of the exchange interaction with the ls and 2s shell on 
the function {399.9 was determined in a stmilar way by substituting the 
function /$o.o into the operator Wf. 

The final values of a= were found to be a* = 12,a~ =—9. The cross- 
scction averaged over the spins could be determincd from known a* anda’: 


<a) = n[3(a7)? +(a*)*] 


It was found to be 380za%. Moreover, it is possible to use these values of 
a* and a™ to detcrmine the amplitude when the scattering process is accom- 
panicd by the spin flip of the valency clectron. 


Na( t )+e({) = NaC} )+e(t) 


Let ,, and ji,. be the spin z-components of the valency and incident clectrons 
respcctivcly. Spin flip involves the transition 


(fs, Its) mo (fls1, — fs) (32.7) 


Since fig; +f,. must be conserved, this transition can only occur when 
Hs ths. = 0. 

Let us cxprcss the scattcring Icngth a(piy. feo. —hg,. —f,2) for the 
process defined by (32.7) in terms of a* and a~. Instcad of the symbols 
+ and — we shall write (0) and (1) in accordance with the magnitude of 
the total spin. We then have 


ACLs, —~ Hest, — Hsis fst) = is ee ar aL Can Oe (32.8) 
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Substituting the values of the Clebsch-Gordan coefficients we obtain 
(Het. —Hsts — sto st) = 3(a? — a) (32.9) 


The corresponding effective cross-section is o = 440nag, which is con- 
sistent with the experimental result ¢ > 260za6 [97]. 


33, EXCITATION OF THE SODIUM ATOM 


The excitation of the 3p level of sodium was considered in [94] on the 
modified Born approximation (§ 18). In the equations for the partial waves, 
only the matrix element V, 4, tysn’1)’.4.7 Was taken into account, i.e. the 
3s—3p interaction. In this matrix element only the dipole part was retained. 
The latter 1s proportional to 


1 
(Yi)3s.3p = > } R3,R3,r° dr + | R3,R3,rdr (33.1) 

r 

0 r 
This quantity was then replaced by its asymptotic form 
3 
(y)35, 3p = ne R3,R3,P dr (33.2) 
0 


It would appear that only the long range part of the intcraction was 
taken into account in this particular calculation. However, this long range 
part plays the dominating role in the case under considcration. 


80 , 
60 
2 
’ 40 
20 
0 2 4 6 8 Te 


E (eV) 


Fic. 46. Cross-section for the excitation of the 3p-state of the sodium atom (in units 
Of mao7): 1—calculated, 2—experimental 


Figure 46 shows the numerical results of the calculation and the expcri- 
mental data as given in [94]. 


Chapter 10 


EXCITATION OF MULTIPLY-CHARGED IONS 


34, EXCITATION OF THE 25 AND 3s STATES OF HYDROGEN-LIKE JONS 


WHEN an electron collides with a multiply-charged ion, the unscreened 
Coulomb field has a considerable effect on the electron. It is then convenient 
to use an approximation in which the Coulomb field is taken into account 
exactly, whereas the effect of the electrons in the electron shell is regarded 
as a small perturbation. This approximation involves the replacement of 
Fy in (31.37) by gy. If, in addition, the exchange interaction is also 
neglected, then the excitation amplitude becomes 


1 i 
GAoA = a7 | Da Vado DP Ao dr (34.1) 
1 


Calculations based on this formula were carried out in [98] for transitions 
to the 2s and 3s states in hydrogen-like ions and for optically allowed 
transitions on the dipole approximation. 

Let us consider the excitation of hydrogen-like ions. The matrix element 
V 4a, 1S then given by 


Vado = | WaoWa dr’ (34.2) 
rr" 


where w,, and wy, are the wave functions for the initial] and final states of 
the ion. We separate out in (34.1) and (34.2) the dependence on the nuclear 
charge Z. The functions w,, and Ww, contain the quantity 


p=Zr (34.3) 


as an argument, while the normalizing factor is proportional to Z?. There- 
fore, V,,, may be written in the form 


Vado = ZU 44,(P) (34.4) 


The function m, contains as its arguments the quantities Zy/A, and kyr 
where Z, is the charge on the ion, Le. 


aa 


0am 04(cakyt) 34.5) 


° 


tA 


If we now express r in terms of p in accordance with (34.3) and substitute 


ae 34.6 
irs (34.6) 
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we have 
Zo 
Pa = Pal 5 Pa-P (34.7) 
PA 
For hydrogen-like ions Zy = Z—1, so that 
J 
a 
Z 
Pa= Pa :Pa-P (34.8) 
A 
Substituting (34.4) and (34.8) into (34.1) we have 
1 Fe 
Lod(Paos PA) =~ 573 | aU ara Pao dp (34.9) 
The cross-section is then given by 
1 p, l 
I 4oa = 73 7, 2 M|7dQ (34.10) 
where 
J | 
1—— 1-—— 
z Z Z 
M = | 6. >» Pa-P U 4c PAo\ > Pap: P | dp (34.11) 
‘ A AAo Pao 
| O}- 
sett 
N 
31S 0-5 
Coa ae 
O 90 180 
6 (deg) 


Fic. 47. Angular distribution of electrons in the case of the excitation of the 2s-state 
of the helium ion 


When Z > | we can neglect Z~' in the arguments of ~, and o,,, and M 
is then independent of A. According to (34.10) the cross-section is then 
inversely proportional to Z*. The integration in (34.11) is carried out numeri- 
cally. The angular distributions of the scattcred electrons are anisotropic, 
and when the encrgy of the incident clectron is close to threshold, the 
maximum of the angular distribution lies near 180°. As an example, Fig. 47 
shows the angular distribution for Z = 2 and the Is—2s transition at 
different encrgies. Figure 48 shows the total excitation cross-scction (obtained 
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by numerical integration of the differential cross-section) for the Is—2s 
transition for different Z. The results of calculations using the Born approxi- 
mation and the distorted wave approximation [99] are given for the five-fold 
ionized carbon atom. In the units chosen in Fig. 47 (p = k/Z along the 
abscissa axis) all the thresholds for the excitation of the 2s state lie at the 
single point p* = 0:75. We note that the excitation cross-section at the 
threshold is finite. The reason for this was discussed in § 23. 


k 


Zz 


Fic. 48. Excitation cross-sections for the 2s-state in hydrogenic ions of different nuclear 
charge: points—distorted wave approximation (in units of za9?); broken curve—Born 
approximation 


35. OPTICALLY ALLOWED TRANSITIONS ON THE DIPOLE APPROXIMATION 


Let us consider now the excitation of optically allowed transitions. 
If we replace the matrix element V,,, by its asymptotic expression for 
large r, we have 


; 
Vado ~ 73 Daao (35.1) 


where Daag, 1s the matrix clement for the dipole moment. The quantity 
Daa, is in its turn related to the oscillator strength for the particular transition. 
The excitation cross-section may be expressed in terms of the oscillator 
strength saa, as follows [98]: 
Shek AE E 
Cig Se Ob a ee 33:2 
me Cheers & AE a) 


where AE is the transition energy, F is the energy of the incident electron; 
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Jo IS given by 


& ) ee “1 ¢in, in’ 7 (35.3) 
O9o( =, — | = 22° — SS -——- 0X1) —— | F Cin, in’, 1, x ; 
°\'2?’ AE E (e2"—1(e2""—1) dng) P80 
where 7 
SS 
/2AE 
Z 
y= - 
\/2(E— AE) 
4nn’ 
S22 35.4) 


XG 
(1y' —11)° 
and Fis the hypergeometric function. 


4 


wie 
\ x 
512 
4 
O 
1 2 3 
£ 
AE 
Fic. 49. A plot of a) as a function of E/AE (points represent the distorted wave 
approximation) 


Figure 49 shows go in units of za) as a function of E/AE for different 
AE/Z*. This figure also shows (AE)*o,,4/54,4 calculated from the data in 
{99}, which gives the calculated cross-section for the transition ls—2p 
in five-fold tonized carbon ions using the distorted wave approximation. 
A detailed table of oo for different values of the arguments is given in [98]. 


CONCLUSION 


THE aceount presented in this book is in some respeets incomplete. This is 
unavoidable beeause the theory of eleetron-atom eollisions is in the proeess 
of rapid development, and in many eases only preliminary results are avail- 
able. Nevertheless, it is possible to discern general trends in the development 
of the subject, and this enables us to indieate eertain problems which will 
remain of eurrent interest. 

Firstly, we must mention the general analysis of cross-seetions near 
threshold with allowance for eoupling between channels and for polari- 
zation of the atom. Studies of resonanees in the excitation eross-seetion for 
a particular atomic state near the excitation threshold for the next state (with 
higher energy) are of particular interest. The M-matrix formalism is con- 
venient in such investigations (Section 12). Seeondly, we must reeall studies 
of the asymptotic behaviour of the wave function for the p + e + e system 
when both eleetrons are in the eontinuous spectrum. This is neccessary for 
the elucidation of the ionization eross-section near threshold. Caleula- 
tions of eross-seetions for the elastie and inelastic scattering of electrons by 
hydrogen atoms show that partial waves with large L play an appreciable 
part even relatively near threshold. As indicated in [61], the develop- 
ment of approximate methods of ealeulation which do not require the 
partial wave expansions will be of major importanee in this eonneetion. 

In calculations of eross-sections for the scattering of electrons by com- 
pheated atoms, one eneounters the unavoidable problem of the rational 
ehoiee of approximate atomic funetions, and also the problem of the general 
strueture of expressions describing the electron-atom interaetion when 
approximate atomie functions are employed (this was considered in Section 
26 in connection with the form of the exchange operator). 

Among the current problems in collision theory, there is the dcvelop- 
ment of approximate formulae for estimating the cross-seetions (possibly 
with the aid of some cxperimentally obtained parameters), and also approxi- 
mate formulae which are convenient in the integration of the velocitics of 
the incident electrons over the spectrum. Simple formulae of this nature 
may be obtained on the basis of elassieal mechanics. In this connection we 
must mention the work reported in [100] and [101]. 

Finally, the sueeessful applieation of the methods of quantum ficld 
theory in a number of branehes of thcoretieal physies suggests that the 
introduetion of these methods into the theory of atomie eollisions will turn 
out to be fruitful. 


158 


REFERENCES 


-ONmrpEEZese 


. A. Lippmann and J. Schwinger, Phys. Rev. 79, 469 (1950). 
. N. Hack, Phys. Rev. 96, 196 (1954). 
. Gell-Mann and M. L. Goldberger, Phys. Rev. 91, 398 (1953). 
. E. Moses, Nuovo Cim. 1, 103 (1955). 
. Eckstein, Phys. Rev. 101, 880 (1956). 
.G. Newton, J. Math. Pigs: 1, 319 (1960). 
. B. Robinson, Phys. Rev. 117, 1281 (1960). 
. A. Bethe, Phys. Rey. 76, 38 (1949). 
. F. O’Malley, L. Spruch and L. Rosenberg, J. Math. Phys. 2, 491 (1961). 
. F. Drukarev, ZhETF, 19, 247 (1949). 
. D. Faddeev, ZNETF, 35, 433 (1958). 


; A, Klein and C. Zemach, Ann. Phys. 7, 440 (1959). 

. W. Kohn, Rev. Mod. Phys. 26, 292 (1954). 

. G. F. Drukarev, Vestn. LGU, No. 22 (1958). 

. G. F. Drukarev, Zh. vychislit. matematiki 1 matem. fiziki, 1, 11, 33 (1961). 

. Yu. I. Demkov, DAN SSSR, 89, 249, 1953. 

. Yu. I. Demkov, ‘‘Variatsionnye piintsipy v teorii stolknovenii Fizmatgiz” 


(1958). ‘*On the Variational Principle’’, Pergamon Press, London, 1963. 

Schwartz, Ann. Phys. 16, 36 (1961). 

. spruch and L. Rosenberg, Piiys. Rev. 116, 1034 (1959). 

. D. Kamenetsktl, Opt. i spektr. 9, 111 (1960). 

.M. Yavoiskii and V. D. Kamenetskit, /zv. vnzov MV i SSO SSSR, Fizika, 
No. 4, 66 and 74; No. 5, 26 (1960). 

. A. Lippmann, Phys. Rev. 102, 264 (1956). 


© 


. T. Epstein, Phys. Rev. 106, 598 (1957). 

. D. Faddeev, ZHETF, 39, 1459 (1960). 

. Gerjuoy, Phys. Rev, 109, 1806 (1958). 

. Castillejo, I. C. Percival and M. J. Seaton, Proc. Roy. Soc. A254, 259 
(1960). 


L 
V 
B 
B 
. L. L. Foldy and W. Tobocman, Phys. Rev. 105, 1099 (1957). 
S 
L 
E 


. I. C. Percival and M. J. Seaton, Proc. Camb. Phil. Soc. 53, 654 (1957). 
. M. H. Ross and G. L. Shaw, Ann. Phys. 13, 147 (1961). 

. R.H. Dalitz and S. F. Tuan, Ann. Phys, 10, 307 (1960). 

. A. I. Baz’, ZhETF, 33, 923 (1957). 

. R. G. Newton, Phys. Rev. 114, 1611 (1959). 

. G. F. Drukarev, ZhETF, 32, 601 (1957). 

. R. Damburg and R. Peterkop, Proc. Phys. Soc. 80, 1073, 1962. 

35. 
36. 
37. 
38. 
39, 


M. Gailitis and R. Damburg, ZHETF, 44, 1643 (1963). 

S. Chandrasekhar and F. Breen, Astrophys. J. 103, 41 (1946). 
M. J. Seaton, Proc. Roy. Soc. A 241, 522 (1957). 

P. Morse and W. Allis, Z. Piiys. 70, 567 (1931). 

T. L. John, Proc. Phys. Soc. 76, 532 (1960). 


40. G. F. Drukarev, ZhETF, 25, 139 (1953). 


159 


83. 
84. 


rer FPF BAZMONMe<yr rane 


REFERENCES 


. R. Marriott, Proc. Phys. Soc. 72, 121 (1958). 
. M. H. Mittleman and K. M. Watson, Phys. Rev. 113, 198 (1959); 116, 920 


(1959). 

A. Lippmann and H. M. Schey, Piiys. Rev. 121, 1112 (1961). 

. Smith, R. P. McEachran and P. A. Fraser, P/ys. Rev. 125, 553 (1962). 

G. Burke and H. M. Schey, Piys. Rev. 126, 147 (1962). 

. Temkin and J. C. Lamkin, Phys. Rev. 121, 788 (1961). 

. D. Ob’edkov, ZHETF, 43, 649 (1962). 

. Temkin, Phys. Rev. Letters, 6, 354 (1961). 

Rosenberg, L. Spruch and T. F. O'Malley, Phys. Rev. 119, 164 (1960). 

. schwartz, Phys. Rev. 124, 1468 (1961). 

Gerjuoy and N. A. Krall, Péys. Rev. 119, 705 (1960). 

. A. Krall and E. Gerjuoy, Phys. Rev. 120, 143 (1960). 

Bederson, J. Hammer and H. Malamud, New York University Technical 
Report. 

. T. Brackmann, W. L. Fite and R. H. Neynaber, Phys. Rev. 112, 1157 
(1958). 

. B. Gilbody, R. F. Stebbings and W. L. Fite, Phys. Rev. 121, 794 (1961). 
A. Vainshtein and G. G. Dolgov, Opt. i spektr. 7, 3 (1959). 

D. Landau and E. M. Lifshits, Kvantovaya mechanika, Gostekhizdat 
(1948). 


. J. Leech, D. Bates, A. Fundaminsky and H. S. W. Massey, Pii/. Trans. Ro). 


Soc. A 243, 93 (1950). 


. M. J. Seaton, Proc. Phys. Soc. 77, 174 (1961). 

. V. M. Burke and M. J. Seaton, Proc. Phys. Soc. 77, 199 (1961). 

. V. M. Burke and K. Smith, Rev. Afod. Phys. (1962). 

. R. J. Huck, Proc. Phys. Soc. A 70, 369 (1957). 

. G. A. Erskine and H. S. W. Massey, Proc. Roy. Soc. A 212, 521 (1952). 

. R. Marriott, Proc. Phys. Soc. 72, 121 (1958). 

. R. Peterkop and R. Damburg, ZhETF, 43, 1763 (1962). 

. W. Lichten and S. Schultz, P/ys. Rev. 116, 1132 (1959). 

. R. F. Stebbings, W. L. Fite, D. G. Hummer and R. T. Brackmann, Pins. 


Rev. 124, 2051 (1961). 


. S. Khashaba and H. S. W. Massey, Proc. Phys. Soc. 71, 574 (1958). 

. I. C. Percival and M. J. Seaton, Phil. Trans. Roy. Soc. A 251, 113 (1958). 
. R. McCarroll, Proc. Phys. Soc. A 70, 460 (1957). 

. T. J. M. Boyd, Proc. Phys. Soc. 72, 523 (1958). 

. D. McCrea and T. V. M. McKirgan, Proc. Phys. Soc. 75, 235 (1960). 

. G. C. McCoyd, S. N. Milford and J. J. Wahl, Phys. Rev. 119, 149 (1960). 
. L. Fisher, S. N. Milford and F. R. Pomilla, P/iy's. Rev. 119, 153 (1960). 

. 8S. N. Milford, J. J. Morrissey and J. H. Scanlon, Piiy's. Rev. 120, 1715 (1960). 
. B. H. Bransden and A. Dalgarno, Proc. Phys. Soc. A 66, 268 (1953). 

. P. Burke, D. D. McVicar and K. Smith, Proc. Phys. Soc. 83, 397 (1964). 
. R. McCarroll, Proc. Phys. Soc. 83, 409 (1964). 

. E. D. Trifonov, ZhETF, 34, 1643 (1958). 

. V. A. Fock, ZhETF, 10, 961 (1940). 

. Yu. N. Demkov, ZAETF, 34, 714 (1958). 

. B. L. Moitsetwitsch, Proc. Roy. Soc. A 219, 102 (1953). 

. B. L. Moiseiwitsch, Proc. Phys. Soc. 77, 721 (1961). 


P. M. Morse and W. P. Allis, Péys. Rev. 44, 269 (1933). 
I. M. Matora, Opt. i speAtr. 9, 707 (1960). 


<P ZOrAMOV.L 


REFERENCES 161 


.S. W. Massey and B. L. Moiseiwitsch, Proc. Roy. Soc. A 227, 38 (1954). 
. M. Morse, L. A. Young and E. S. Haurwitz, Phys. Rev. 48, 948 (1935). 
. J. Schulz and R. E. Fox, Phys. Rev. 106, 1179 (1957). 

. Baranger and E. Gerjuoy, P/iys. Rev. 106, 1182 (1957). 

. Marriott, Proc. Phys. Soc. A 70, 288 (1957). 

. A. Vainshtein and G. G. Dolgov, Opt. i spektr. 11, 301 (1961). 

. F. Drukarev, ZAHETF, 31, 287 (1956). 

. J. Seaton, Phil. Trans. Roy. Soc. A 245, 469 (1953). 

. A. Vainshtein and I. I. Sobel’man, ZAETF, 39, 767 (1960). 

. Salmona and M. J. Seaton, Proc. Phys. Soc. 77, 617 (1961). 

. A. Fock and M. IJ. Petrashen’, P/ys. Zs. Sowjetunion, 6, 368 (1934). 


D. R. Hartree and W. Hartree, Proc. Roy. Soc. A 193, 299 (1948). 
. H. G. Dehmelt, Phys. Rev. 109, 381 (1958). 

. M. Gailitis, Opt. i spektr. 14, 465 (1963). 

. L. A. Vainshtein, Opt. i spektr. 11, 301 (1961). 

. M. Gryzinski, Phys. Rev. 115, 374 (1959). 

. V. L. Ochkur and A. M. Petrun’kin, Opt. i spektr. 14, 457 (1963). 
. L. Vriens, Physics Letters, 9, 295 (1964); 10, 170 (1964). 

. R. Stabler, Phys. Rev. A 133, 1268 (1964). 

. R. Peterkop, Izv. A. N. Latv. S.S.R. 9, 79 (1960). 

. V. Ochkur, ZHETF, 45, 743 (1963). 


